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The vibration-rotation energies of the axially symmetric X Y3Z molecular model are investi- 
gated to second order of approximation, including cubic and quartic anharmonic terms in the 
potential energy, all Coriolis interactions between rotation and vibration, centrifugal stretching 
and dependence of the rotational constants on the vibrational state. The energies are expressed 
in term value form, E=hc(Gy+ Fr— Fc), where Gy, Frand Fe are, respectively, the vibrational, 
rotational, and Coriolis terms which are given explicitly as functions of the molecular constants 
appearing in the Hamiltonian. 


I. INTRODUCTION 


HE axially symmetric XY3Z molecular model, of which examples are the methyl halides, is of 
considerable interest to chemists and physicists and the infra-red and Raman spectra of many 
molecules of this type have been subjects for much experimental study. Various aspects of the theory 
underlying the interpretation of these spectra have been considered by many workers. The normal 
vibrations have been studied by several workers'~ using different types of potential functions and 
other writers ® have studied the symmetry properties of the normal vibrations from the standpoint 
of group theory. Several authors*’? have attempted to set up simplified quadratic potential functions 
to fit the observed vibrational frequencies. Dennison! has given the rudimentary theory of the 
rotational fine structure of the infra-red spectra of XY3Z molecules. Other writers have treated 
certain aspects of Coriolis interactions,’ centrifugal stretching" and accidental degeneracies.'” 

In view of the fact that the theory of the rotation-vibration spectra of XY;Z molecules exists in 
the literature in fragmentary form and it has not been carried to the high order of approximation 
necessary for the interpretation of certain features of experimental data obtainable with improved 

1D. M. Dennison, Rev. Mod. Phys. 3, 280 (1931). 

2G. B. B. M. Sutherland and D. M. Dennison, Proc. Roy. Soc. A148, 250 (1935). 

3 J. E. Rosenthal and H. H. Voge, J. Chem. Phys. 4, 134 (1936). 

‘E. B. Wilson, Jr., J. Chem. Phys. 2, 432 (1934). 

5 J. E. Rosenthal and G. M. Murphy, Rev. Mod. Phys. 8, 317 (1936). 

°H. H. Voge and J. E. Rosenthal, J. Chem. Phys. 4, 137 (1936). 

7Z. 1. Slawsky and D. M. Dennison, J. Chem. Phys. 7, 522 (1939); J. W. Linnett, J. Chem. Phys. 8, 91 (1940). 

SE. Teller and L. Tisza, Zeits. f. Physik 73, 791 (1932). 

°E. Teller, Hand- und Jahrbuch der Chem. Physik (1934), Vol. 9. 

10M. Johnston and D. M. Dennison, Phys. Rev. 48, 868 (1935). 
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techniques, it has seemed desirable to carry out a unified treatment of the rotation-vibration 
energies of such molecules and to extend the calculations to second order of approximation in such 
a manner as to include, in addition to the energies associated with the harmonic oscillator and rigid 
top approximations, contributions arising from cubic and quartic anharmonic terms in the potential 
energy, all types of Coriolis interactions between rotation and vibration, centrifugal stretching and 
changes in the rotational constants arising from vibration. In subsequent papers the expressions 
developed in this paper will be applied to the interpretation of some new experimental observations 
on the infra-red spectra of the methyl] halides. 


Il. THE NORMAL COORDINATES 


The model adopted for the axially symmetric XY;Z molecule consists, in its equilibrium con- 
figuration, of a symmetrical XY; pyramid, similar to that assumed by the writer" for the pyramidal 
XY; model, with a Z particle situated outside the pyramid on its symmetry axis and adjacent to 
the X particle. For describing the positions of the particles a body-fixed rectangular coordinate 
system, xyz, is adopted whose origin lies at the center of gravity of the model, whose positive 2 axis 
coincides with the symmetry axis of the model and points from the X particle toward the Z particle, 
and whose xy plane is parallel to the X; base of the pyramid. The particles, which are assigned 


— — fad : 0 0 0 
indices 7 (i=1, ---5), for reference, have the masses m; and equilibrium coordinates (x;, yi, 2i) 
shown in Table I. The equilibrium coordinates are illustrated in Fig. 1. 

The instantaneous position of the ith particle during oscillation is given relative to the body-fixed 


, 0 0 0 . jus 
coordinate system by (x;, yi, 2:) = (xix, yity's, z:+z,;) where the primed quantities denote com- 
ponents of displacement from the equilibrium position. From a consideration of the six equations 
satisfying the Eckart conditions": ® and the expressions for the relative displacements of the particles 
during small oscillations it is evident that the following is a suitable choice of intermediate generalized 
coordinates} for discussing the rotating and vibrating model: 
’ ’ -} 
x=x4—[(xi+ao+x3)/3]; €={3 (2x1 —x2—2x3)+(ys—vy2)} /2 
’ ’ , 3 ’ 
y=ya—L(nityet+ys)/3]; 2={-3 (291-92 —ys) + (13-42) } /2 
’ ’ , P =} 
a=%—[(ateete)/3]; 6={—3° (2x1—x2—45)+(y3—y2)} /2 
; p=Xs—%4; Q=Ms— Ms = — Bi. 
13 W. H. Shaffer, J. Chem. Phys. 9, 607 (1941). 
14 C, Eckart, Phys. Rev. 47, 552 (1932). 
18 E. B. Wilson and J. B. Howard, J. Chem. Phys. 4, 260 (1936). 
+ These intermediate coordinates and the similar ones of Slawsky and Dennison, reference 11, satisfy all the Eckart 
conditions and are therefore suitable for discussion of the rotating and vibrating model. The coordinates of Rosenthal 


and Voge, reference 3, do not satisfy the Eckart condition requiring conservation of internal angular momentum in zero 
approximation and could not readily be used for study of interactions between rotation and vibration. 
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0 0 0 
TABLE I. Masses m; and equilibrium coordinates x;, yi, 2i. 














i mi xi® yi? si? 
Y 1 my, —do 0 —= C4 
Y 2 my, +d o/2 —3iao/2 —C1 
Y 3 my, +ad/2 +34a9/2 mm G5 
X 4 Me 0 0 —C4 
n Z 5 


Ms 0 0 +6; 








The kinetic energy of vibration, Ty, is given in terms of the above intermediate coordinates by 
the expression, 


Ty = (1/2) {mi(pP+ @) + u2(@ +9) + 2us(pe+ Gy) +im(E+ P+) +s? + us2 + 2usrZ}, (1) 


where , 
ws=[ms5(3my+m4)/M]; wr =4s t+ Oma" 

n us=[3mi(mat+ms)/M];  w2=ys +68? 
- Me = [3mym;/M ]); M3 = ue + Om a8 
il in which M=(3m,+ms+ms); a= —(ms5¢5/3m,a0); B= —(c1/ao). It can readily be shown that the 
d generalized quadratic potential energy function, Us, is 
1S Up = (1/2) {ki( + 0?) + ho(x? + y?) +hs(p? +9?) + 2ka(pxt+qy) 
+2ks(éx+ny) + 2ko(EP+ 0g) + mike? + m2" + ngr? + 2nyrz+2ns6zt+2nefr}, (2) 


where the constants hi, ke, etc., are generalized force constants which may be expressed as functions 
of the constants occurring in any special type of potential function such as valence or central force 
n- types.* 





al The axial model XY;Z belongs to the symmetry point group C;, and has three non-degenerate 
to parallel normal modes of class A; and three degenerate perpendicular modes of class E as can be 
te shown by the group theory treatment of Wilson‘ or Rosenthal and Murphy.’ Approximate pictures 
is of the normal modes are given by Sutherland'® and Wu” where the single modes are arbitrarily 
e, called w1, w3, ws, and the double modes we, ws, ws, following the notation of Sutherland and Dennison.’ 
od From an inspection of the Lagrange’s secular determinant for this case, which breaks up into 
4 three threefold steps, it is obvious that the intermediate coordinates ¢, r, and z are associated with 
; the single modes and the sets (£, p, x) and (n, g, y) are each associated with one set of the double 
ad modes. The transformations to normal coordinates can be somewhat simplified by the following 


intermediate transformations: 


3 ’ P a ’ 4 , 
ns meE=Q1; mn=R, me = Si, 
les 5 -3 ’ , 5 4 ’ , Hy 3 v v 
ed box =2 (o102+0203), boy =2 (o1R2+o2R3),  s2=2 (03S2-+045s), 
i 3 ’ ’ Hy 3 , ’ H -} e , 
Mp= 2 (o1Q2 — a203), Mig >= 2 (0:Re+o2R3), pyr =2 (o3.S2— 0453), 
where 
o1= {1+[ms/(mime)* }}!; o3= {1+[ue/(uams)? ]}}; 
o2={1— Cus/ (ume)? 4s og={1- Cus/ (mus)? ] ie 
* Computations of the constants in simplified potential functions for the methyl halides have been made by Voge and 
Rosenthal, reference 6, and by Slawsky and Dennison, reference 7, (SD), all basing their work on the generalized potential 
functions given by Rosenthal and Voge, reference 3; they obtained complex values for some of the force constants. The 
{ writer has made computations, basing his work on the generalized expression (2) above, on the constants in a potential 
art function similar to that of SD; in the results, which will be published later, the difficulty involving complex roots does 
hal not seem to arise. 
ero 16 G. B. B. M. Sutherland, Infrared and Raman Spectra (Methuen, London, 1935). 


17 Ta-You Wu, Vibrational Spectra and Structure of Polyatomic Molecules (Prentice-Hall, New York, 1940). 
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As a result of this transformation 7'y becomes 


Ty =(1/2) x (QF +R +38!) (3) 
and Up» becomes : 
Ur= (1/2) & CK(Qe + RIAN S! +2.E (Ki Qe + RR + Nw SiSeIN, (4) 
where : - 
K,=(k:/m), N= (m/m,), 


K,= (01/2) (ke/ pe) + (es/ 1) + (2ks/(urme)?) J, No= (03/2) (me / Ms) + (13/ us) +(214/(uaus)*) |, 


Ks=(02/2)[ (e/a) + (ks/4s) — (2ks/(urme)!)], Ns = (05/2) (rt2/us) + (1t5/ ma) — (204/ (uaas)*) J, 


Ky= o1(2m;) (ua Rstun e), Nye= ox(2ms) (us Ms-+ua Ns), 
Kis=02(2m) (uo ks—m ke), Nis=o4(2m) (us Ms— a M6), 
Ky3= (o102/2)[(Re/ m2) = (Rs/p1) |, N23 = (o304/2)[ (n2/ Ms) _ (n3/ ma) |. 


The generalized coordinates Q; and R; (¢=1, 2,3) are associated with the twofold degenerate 


modes (we, ws, ws) and the S’, with the single modes (w, ws, #5). The normal frequencies of oscillation 
w; (in cm=) are related to the roots \; of Lagrange’s determinantal equation by the expression, 


—13, : P . : ; ; 
w;=(2mc) \j, in which c denotes the velocity of light, and where the \; satisfy the following relations: 


Needs Ap= Ki-+Ke+Ko, 
iat Dale Nake Kile 4 Kes + KK. —Ka—Ko~-Ko. 
bade Si Malls ~ReKs~Kele—KnK, 
hibAs+As= Mi + Ne+Ns, 
Arda-EAgde-t Ag = Ny No-+ N2Ns+ NgNi— Ni2— Nis— Nos, 


Adsds = Ni NoN3— NieNs— NisN2— NosM1. 


The transformation to normal coordinates, Q;,"R;, S;, is accomplished by means of the following 
substitutions: 
- Q=LiaijQ;; Ri=DXj aijR; (i=1, 2, 3; j=2, 4, 6) 
Si= Le BueSr (G=1, 2, 3; R=1, 3, 5) 

where a;;=A;A;; and By = B,Bi if 

Ay;=(K2—\i)(Ks—\;)— Ki, Bu=(N2—Ms)(Na—s) — Nu 

Ao;=KiKs— Ke(Ks—\j), By. = NsNs— No(N3— x), 

A3;= KiKe—Ks5(K2—\)), Bs. = NiNe— Ns(N2— dx), 


A;= (hited hed B.=(But+But+By). 
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wm 


The transformation coefficients, a;; and B;,, satisfy the relations 
Liai=Liaj=1; Li aijoiye = Lj aijav;=0, 
2 2 
Di B= De Bx=1; Li BuBu = Le Bu Birr =O. 


, ‘ . , . Soot 
If dimensionless coordinates qj, 7;, s:, are introduced by the relations, Q;=(h /d,;)q;, etc., where 
h=(h/2m) and h is Planck’s constant, the vibrational energy of the system can be written, to zero 
approximation, in the Hamiltonian form as 


Hov=hel Dj wil (ba +Pr)/B)EG AN) IF Le wil (be /b) +5], (5) 


where ps, = (dT y/ds,) and is the momentum conjugate to s;, etc. The pairs of coordinates, g; and rj, 


associated with the doubly degenerate modes can be conveniently replaced by polar coordinates 
according to the substitution 
i=pj COS $j; 7j;=p; SIN ¢;. 


It can be shown from consideration of the symmetry properties of the normal modes that the 
cubic and quartic parts of the anharmonic potential energy function contain the following terms: 


Ur=he De sel Bewesi +X; Beivi + X Bex-nSe-, (6) 


“Se 


2 2 2 
U,= he{ ar SCY KickeSie A ie (Vice Ske HV cane SESK?) | 


k'>k 


2 2 2 2 
+205 oi viii DL viiriveey + De veniiSet+ De Lo viinwsise J}, (7) 


">i aaa 


where the constants Bxx, Yer, EtC., are expressed in cm. 


III. THE INTERNAL ANGULAR MOMENTUM 


The Coriolis interaction between rotation and vibration gives rise to the term, —[(p,P:/Joz) 
+ (pyPy/Toy) +(p-P-/Io:) |, in the quantum-mechanical Hamiltonian, where p= (pz, py, p.) =angular 
momentum associated with vibration, P=(P,, P,, P:)=total angular momentum of the molecular 
system, and the equilibrium moments of inertia are 


Tos = Toy = 3m: c1-+ (ao/2) ]-+macs+ mcs, (8) 


Ioz= 3mydp. (9) 
The components of the internal angular momentum p have the following forms: 


bo= Qi Dee Eel (wx/0;)'r ipo, — (wi/x)*sepr -] 

+ Lin Esl (wy/05) (Gib, +1 iP a,,) — (wil) (ria +4rbr dM}, 
Pu=— Qj (Lee Siel (wx /e5) gids, — (wi/wr)'SeP a] 

—L ir L(wir/oi)(ridr ~ iba )— (wi/ay ripe — 9i'Pa ) J} 


P2= Li (Ei(Qibr —TiPa ALG Fir L@r/wi)(Qibr —1iP a) — (w;/wr) (tba —Girhr ) ]}, 
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where the indices j and 7’ assume values 2, 4 and 6; the index & runs over 1, 3, and 5; and the double 
subscripts on é;;, and £;; occur cyclically as 24, 46 and 62. The constants éj, &);, ¢; and ¢;; which 
occur in the above expressions are defined in the following manner: 

E ix = 01; 820x022 + 043 83.033 + O12 83423 + 03 {824.032 + (a2;012 + 23 ;013) Bix, 

fj = (aj; ae 01 02; )\y2-+ (a1 ;@3; = 1 jQ13; 113, 


9° 


9 ° 
6) = aajb02 + «3 jb33 + Zee j013 bos — a1 ;, 


C57" _ Ot j0¢2 j'D 22 + a1 043 ;'D33-+ (a2 j0%3 ;7 + a2;013;) 23 — &1jQ1;', 
where 


A222 = (0103/2) (52+ o%+012+051), a23 = (o104/2)(o52—o41 —os2+051), 


33 = (0203/2) (os2+01—o12—051), A32= (0203/2) (052 — 041 +042 —951), 
2 2 2 2 2 
bee = (01/2) (o52+0411+ 2042041), bo3 = b32 = (o102/2) (052 — 041), 


2 2 2 
b33 = (02/2) (os2+ou = 20420141), 
if o5.= (us/me2)?*; on = (us/p1)?; O42 = Me (Maule); 05: = we ( Msp). 
The parameters ¢; satisfy the relation, 


2 2 
Sot Sst 56 = (uaus — me) /(uime— ws) = (Lo:/2Loz), 
which has been pointed out by Johnston and Dennison.” Other identities involving the &j, &)., 
¢;;, and ¢; may be derived with the aid of the definitions given above. 
Certain constants are defined below which occur in the rotational part of the Hamiltonian and 
are used in Section V in connection with the eigenvalues of the Hamiltonian operator: 


€;= —2(a2;05+03;06), fu=(1— (1/2)81.— di En], 
f= (Bur— 2Bx.07 — 2808), giv = [1 —(3/2)a1;—6;- Le Bie, 


fig=(l- (1/2)a1;-Ei ta-Zp i hjj= ftr¥Seau~t;-Er tek 


where 
o5= o1(3m/2) (ar +Bu2 ), o7= o3(3m/2) (Bus +a ), 


— ~f 7 il 
o¢=02(3/2) (Buz —am ), o3=o4(301/2) (Bus —aps ). 
IV. THE QUANTUM-MECHANICAL HAMILTONIAN 


The quantum-mechanical Hamiltonian operator appropriate for studying the rotation-vibration 
energy states of a polyatomic molecule can be set up according to the method outlined by Wilson 
and Howard® in the form 


H=(1/2) D w(Pa—Pa)base *(Pa—pa)+(1/2)ut Xe pew pet Uot Uit U2t+---, (10) 
a, 8 
where the quantities P, and p, are the operators associated with the x, y, and z components of total 
and internal angular momentum, respectively; the quantities u and wag are functions of the normal 
coordinates and are related to the instantaneous moments and products of inertia; the p,; are oper- 
ators associated with the momenta conjugate to the vibrational coordinates; and the quantities Uo, 
U, and U; are the zero-, first- and second-order parts of the potential function. For the purpose of 
computing the energies to second order of approximation the expansions of the quantities u and pas 
are broken off with terms quadratic in the vibrational coordinates. The details of the procedure for 
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setting up the Hamiltonian operator H for the XY3;Z model are similar to those in papers on the 
XY, model by Shaffer, Nielsen and Thomas'* (SNT) and the bent XY2 model by Shaffer and 
Nielsen. 

A contact transformation of the type suggested by SNT may be applied to the Hamiltonian 
operator HT to yield a transformed operator H’ given by the expression, H’=THT~-', where 
T=exp [(—1)*«S] and S is a function of the operators occurring in /7. The transformed Hamil- 
tonian, /7’, has the form 


W ath path bete-->, (11) 


where ¢ is the parameter of smallness according to which the Hamiltonian is separated into its 
various orders of magnitude, 


Hy=Mp, H,=H,—(— 1)}UI)S— SH), H,=H.+[(— 1)3/2)S(, +H) — (+h) S]. 
By means of such a contact transformation it is possible to reduce very much the number of terms 


, , 
in H, by transferring them to H2; this procedure simplifies the derivation of the second-order con- 


tributions to the energy expression. It can be shown that all terms can thus be eliminated from Hy 
except the Coriolis interaction terms arising from degenerate oscillations and certain terms in U, 
in the case of accidental degeneracies. Details of the application of the contact transformation have 
been given by SNT and by Silver and Shaffer.” 

The basic wave functions for the perturbation calculation of the energies are similar to those 
given by Silver and Shaffer.” The quantum numbers defining the rotational states are denoted by 
J, K and M following the notation of Dennison;! the vibrational quantum number associated with 
normal mode w, is denoted by v, and, in the case of the twofold degenerate frequencies, the quantum 
number associated with the internal angular momentum is denoted by ?,. 


V. THE ENERGY EXPRESSION 


The zero-order Hamiltonian contains the terms, 


Ho=Hov+(1/2)[(Pz/Toz) + (Po/Ioy) + (P./Ioz)] (12) 


and has eigenvalues, 


6 
Eo=he{ © wnlon+ (gn/2)J+BeJ(J+1)+(C.—B.)K*}, (13) 


n=1 


where B,=(h/82°cIoz), Ce=(h/82°cIo.) and g, is the degree of degeneracy of normal frequency 
Wn (Q1=23=25=1; go=21=26=2). 


, 
The first-order Hamiltonian, 1;, can be written, for the case of no accidental degeneracies, in the 


form 


My = — Di (Fi/Toz) PoiP:, (14) 
where p4;=(¢ibr —rjPq ); the first-order contribution to the energy is 
j i 


E,= —heFo, (15) 
where Fco= 2C. > €,1;K. 


, 
The second-order Hamiltonian, H2, for the XY3Z model is a very complicated expression and will 
not be given here since it is similar to the corresponding expression given by the writer” for the 
18 W. H. Shaffer, H. H. Nielsen, and L. H. Thomas, Phys. Rev. 56, 895 (1939). 


19 W. H. Shaffer and H. H. Nielsen, Phys. Rev. 56, 188 (1939). 
20S. Silver and W. H. Shaffer, J. Chem. Phys. 9, 599 (1941). 
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pyramidal XY3 model. The second-order energy, E:, can be obtained by the usual perturbation cal- 


culation in a manner similar to that employed earlier by the writer.” 


When the contributions Eo, E, and E; are added together the total rotation-vibration energy E 
can be written, to second-order of approximation, in the simple form 


E=he{Gv+Fr-—Fc}, (16) 


where F¢ is the quantity defined in (15), Gy and Fp are, respectively, the vibrational and rotational 
term values. The quantities Gy and Fp are given by the following expressions: 


Gy=Got+ y { onl n+ (gn/2) ]+XnnLOnt (gn/2) P4- ¥ xnrnLone + (gn'/2) Lent (gn/2) J} 


n=1 n'<n 
+25 (Gi(G-ID+ LD Gyiljlj}, (17) 
i'<i 
Fr=ByJ(J+1)+(Cyv—By)K?-—DyJ?(J+1)?-—DixJ(J+1)K?—DxK', (18) 
in which the rotational constants By and Cy effective in the various vibrational states have the forms, 
6 
By=B,{1-am+> nL Unt (gn/2) ]}, (19) 
n=l 
6 
Cv=C.{1—Bo+X Bulent (gn/2) J}. (20) 
n=1 


The constants occurring in (17) are 


Go= (1/16) Xe {Oveeee—7(Brre/on) +3 X [or Beree/(4er—ww) J} 
. 


+23 0.-F SiC 42GB.) -Ze En} 


7’ 


xn = (1/4) | Ov enee — 15(Binx/cx) +3> [ (Sco - 3con-) Brae /con” (Ae, _ wy) ]} ; 
7 


wi7= (1/4) {Orsini — Le [ (Bw; — 30x) Brs/eoe(40j— ex) J}, 
Xen =Verk ke — 3L(BrteBerks/ wr) + (Berern Bunn /we) ]— 24 [onBern/ (400, — wy) + [oom Bese / (Acar? —w)}} ; 
2305 =Viririst (Ev Bet Sir Co) (wi/wy) + (wr /) 1— Cee (BeiiBeirir/on); 
Xjn= — 3 (Bier Buji/wr) — py (Be siBr-nx/tox’) — 2[00)81js/ (Ato; — coe) J+-Yenss + EBL (co;/cor) + (con/w;) J, 


Gis= —(1/2)vitii— Le [onBej;/4(40;— on) ]+6;Ce, 
Gig = 26 FC 


The centrifugal expansion constants Dz, Dyx and Dk, of (18) are defined as 
Dy =(1/4)(3tez2z+ 27 2y2yt Tez); 
Dyx=(1/2)(—3tesse— 2Teysy— Teayy t+ BT es2e 47 2x22); 
De=—Dy—DyKttezes 
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where 
- sm 2,2 2,2 
‘ome (B./C.) Li (a1;/w;) + Le (fi-/ wi) }, 
2 2 
Tzzzz— Cc wr (4B1;./wx), 
Txyxy (B./C.) Di (a;/«;), 
Taexyy — Trrrz 2 Tryzrys 
Tzzrz = B.C, oe (e; ‘ «;) , 
Trrzz— B.C, dD - (2Bu-fi fem). 
The constants ap, Bo, a, and 6,, occurring in (19) are 
a= (1/2) ( — Trzzz 2Tryzyt Tea T y | om Be 
Bo = (1 /4) (3 Traxzet OTasey —3 Trxryy 1 6Tzsz2) ’ 


ox, = (h/2acopToe) {[(Toe/ Toe) fu — fer +23 CE ix(3eo¢-+0)) / (wx — 0) ]} 
j j 
+ (h/2nc)} { (Sf Brict/ wr ) + De (fer Burn / wn) } ’ 


‘ i i 
Bi = (h/2rcwrTo2)(3Bix) + (h/2mc)}{ (3hiBere/wn) +X Ne Berice/wr)}, 
k’ 


a= (h/2mcwo;Toz) {fi(Toe/ Toe) + (1/2)(€;— fii Bi) FL ir L (30; +0;-)/(wji— 0) ] 
2 2 2 2 2 ; 
+(1/2) Dee Ejel(3@j;+ wx) /(wj— wx) ]} +(h/2mc)! Se (fiBric/wr), 


. i 
B;=(h/2mcwToz) {e;(Ioz/Toz) —hjit+ E hirL(30;+e)/(wj—w5") J} +(h/2mc)* ee Beii/wn). 


Discussions of the selection rules governing transitions between the energy states have been 
given by Dennison,! Wilson,‘ Sutherland,'* Teller,? and others. 
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Infra-Red Spectra of Hydrocarbons 


I. Some Investigations of the Temperature Dependence of Absorption Bands 


W. H. Avery anp C. F. ELtis 
Houston Refinery Research Laboratory of Shell Oil Company, Incorporated, Houston, Texas 


(Received October 14, 1941) 


Theory indicates that the separation of the maxima of the P and R branches of an infra-red 
band should be proportional to the square root of the absolute temperature. This relation has 
been tested by experiments with ethane, propylene and 2, 2-dimethyl butane in which the 
temperature was varied between — 195°C and room temperature. It has been found that the 
relation is approximately valid, the band widths being somewhat less than theory predicts in the 
liquid and considerably less in the solid. The positions of the bands are independent of the 
temperature in 2, 2-dimethyl butane and only slightly affected in the other compounds. Some 
bands of propylene display a pronounced variation in relative intensity with temperature. The 
usefulness of spectra obtained at low temperatures in the theoretical assignment of funda- 
mentals, as a tool for the study of the liquid state and for analytical purposes, is pointed out. 





INTRODUCTION 


JERRUM first showed in 1912 that the 
infra-red absorption bands of diatomic 
molecules possess P and R branches symmetri- 
cally disposed with respect to the center of the 
band and that the separation in wave numbers 
of the maxima of the two branches is given by 


the relation, 
1 sky * 
weary 
m\ I 


K being the Boltzmann constant, T the absolute 
temperature, J the moment of inertia and c, the 
velocity of light. Von Bahr and others have 
demonstrated that increasing the temperature 
causes a separation of the two maxima, in accord 
with the theory. She also pointed out the 
significance of variation of the temperature as a 
means of distinguishing double maxima produced 
by two bands and those due to P and R branches 
of a single band. However, because of the 
extreme temperatures required to produce a 
significant change in band width and the experi- 
mental difficulties involved, these studies at high 
temperatures do not appear to have been used 
in the interpretation of molecular spectra. 
Serious experimental obstacles are not encoun- 
tered in the use of low temperatures, and in view 
of the problems presented by the theoretical 


* See footnote on page 14. 


interpretation of spectra where the fundamental 
bands overlap, it is surprising that low tempera- 
ture studies have not heretofore been attempted 
as a means of obviating the difficulty.t 

The use of infra-red spectra for analytical 
purposes is limited by the amount of overlapping 
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Fic. 1. Low temperature cell. 


+ A few investigations of infra-red spectra of gases and 
liquids have been carried out at low temperatures. How- 
ever, these studies were undertaken in investigations of 
rotation in the solid state. See Ruheman, Low Temperature 
Physics (Cambridge University Press, 1937). G. B. B. M. 
Sutherland and C. K. Wu in an article published after this 
work was begun (Proc. Roy. Soc. 176, 485, 1940) mention 
the possibility of using low temperatures in the interpre- 
tation of molecular spectra. 
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Fic. 2. Spectrum of ethane at room temperature. 


of the bands in the various components of the 
mixture to be analyzed. Since this overlapping 
would be only half as great at —195°C as at 
room temperature, use of low temperatures 
might also be of value in this field. 


EXPERIMENTAL 


The spectrometer has been described.'! The 
absorption cell and optical arrangement used for 
the low temperature experiments are shown in 
Fig. 1. Parallel light was reflected by the front- 
surface aluminized mirror M, to the silvered 
bottom of the cell 1M; then to Mz which returned 
the beam to its original direction. The cell was 
constructed by the American Instrument Com- 
pany, and the surface M; was flat to 6 wave- 
lengths of sodium light. To facilitate filling, the 
cell was mounted on a tripod set in V-ways 
which permitted the cell to be moved and 
returned accurately to its original position. A 
similar arrangement allowed the mirrors M, and 
M; to be removed and replaced without change 
of the focus. Leveling screws were provided for 
both tripod and mirror mounting, and the angle 
of the two mirrors with respect to each other or 
to the mirror M; could be adjusted. The cell 
was 7 cm in diameter and 15 cm in depth. 

Absorption spectra were obtained as follows: 
After the cell had been filled to a known gas 


1W. H. Avery, J. Opt. Soc. Am. 31, 633 (1941). 


pressure, the gas was condensed in the side arm 
S with liquid nitrogen and a given region of the 
spectrum, usually about 2 mu, recorded on the 
chart. The side arm was then warmed, the gas 
condensed on the cooled plane surface 1/3; and 
the same spectral region was again recorded. 
Percentage transmissions were computed in the 
usual manner. This procedure, in addition to 
being convenient and rapid, eliminated errors 
which inevitably arise in regions of small 
absorption from the use of null cells. 

The fraction of light reflected at a liquid 
surface depends indirectly on the absorption 
spectrum of the compound being investigated 
and could become appreciable (~10 percent) in 
regions of strong absorption.? A correction for 
this reflection would be desirable but this is not 
possible in the absence of dispersion data. 

In many of the experiments at the lowest 
temperatures a polycrystalline layer was formed 
which scattered an appreciable fraction of the 
incident radiation. This fraction varied from 
experiment to experiment, the nature of the solid 
formed evidently being sensitive to small differ- 
ences in the rate of cooling. Therefore, in 
preparing the absorption curves, it was assumed 
that regions showing no absorption in the gas 


2 See, for example, the dispersion data for three methyl 
halides given by A. I. Mahan, J. Opt. Soc. Am. 31, 248 
(1941) and dispersion data on water by M. Centeno, J. 
Opt. Soc. Am. 31, 244 (1941). 
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Fic. 3. Low temperature spectra of ethane. 


represented 100 percent transmission in the 
solid; the assumption seemed justified by the 
fact that the solid displayed wide regions in 
which the transmission remained at a high 
(60-70 percent) constant value. Since the cor- 
rection is large, numerical values of the trans- 
mission listed for the lowest temperature may 
be somewhat in error, but the positions of the 
absorption bands and the relative intensities of 
the bands are presumably not seriously altered 
by the correction. 

A temperature of —195°C was obtained by 
immersing the lower part of the cell in liquid 
nitrogen contained in the Dewar flask, D, Fig. 1. 
Other temperatures could be obtained by proper 


TABLE I. Observed frequencies of ethane, with a possible assignment. 


adjustment of the distance between the bottom 
surface of the cell and the upper surface of the 
liquid nitrogen. The KBr window remained at 
room temperature. Temperatures were estimated 
by measurement of the difference in potential 
of two opposed iron constantan thermocouples, 
one immersed in liquid nitrogen and the other 
attached to the outside bottom surface of the 
cell. The temperatures are uncertain by +10° 
with the exception of —195°C. To assure a 
uniform layer the gases were condensed on the 
surface M3; at a temperature where the liquid 
viscosity was not great, then cooled to the 
desired temperature. Direct condensation at 
very low temperatures yielded viscous liquids 














Observed Observed Assigned 
wave-length frequency frequency 
mu cm7} Intensity* cm"! 
12.09 827 Ss Fund 
8.22 1216 vw 1460 —275 
7.25 1380 m Fund 
7.14 1400 —' 
6.80 1470 Ss Fund 
6.50 1538 vw — 
5.84 1712 w 1460+275 
5.46 1831 w 827+993 
4.97 2012 vw 1170+827 
4.50 2222 w 7d 
1 +99: 
4.25 2353 mes { ater 
3.74 2674 w 2960 — 275 
3.64 2750 m 1380+1375 
3.36 2980 s Fund 











Calculated 

trequency 
em*! Dx Disa 

827 E’ Ey 

1185 E' Ey 
1380 A's Av, 
1470 E’ Bie 
1735 E’ pe 
1820 E’ Px 
1997 A"; +A” e+ EF” A jutAout Ex 
2202 E’ | 
2377 A me A Qu 
2340 A';\+A'o+F’ a 
2685 E’ z 
2755 A'’s Aox 
2980 E' Ey 











* Meaning of symbols: s =strong, m =medium, w =weak, vw =very weak. 
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which collected in drops and showed little 
tendency to spread into a thin film.* 


MATERIALS 


Ethane was prepared by Mr. E. Manning of 
this laboratory by hydrogenation of Ohio Chem- 
ical Company’s ethylene at room temperature 
using a platinum catalyst. The ethylene was 
stated to be 99.5 percent pure, the hydrogen 
99.9 percent. Analysis of the ethane by Mr. E. 
Gelus of this laboratory using a method of 
isothermal distillation showed the purity to be 
99.8 percent or better. 

Propylene obtained from the Ohio Chemical 
Company was stated to be 99.95 percent pure. 
The purity was confirmed within the limits of 
error of the isothermal distillation, and it was 
used without further purification. 

2, 2-dimethyl butane was prepared from a 
petroleum fraction by Mr. W. H. Marsh of the 
Rice Institute by repeated azeotropic distillation 
with a 100 plate Hicks-Brunn column. The 
melting point was —102.0°C and the boiling 
point 49.80°C. 


ETHANE 


Because of its significance in the determination 
of the potential restricting rotations of the 
methyl groups, ethane has been the subject of 
several recent investigations.*~> The existence 
of a barrier of about 3000 calories is now well 
established but spectroscopic evidence still is 
conflicting with regard to the orientation of the 
methyl groups. The low temperature results 
were expected to be of interest for their bearing 
on this question. ; 

The spectrum of ethane has been obtained at 
room temperature, at —170°C and at —195°C, 
the pressures used ranging from 50 to 500 mm. 

The spectrum found at room temperature is 
presented in Fig. 2. The wave-lengths and 
frequencies of the bands are presented in Table I 





* In a few experiments at low temperatures the gas was 
first condensed in the side arm, the reflecting surface next 
cooled to — 195°C and the gas then allowed to vaporize to 
the surface. 

3 F, Stitt, J. Chem. Phys. 7, 297 (1939). 

‘ B. L. Crawford, Jr., W. H. Avery, and J. W. Linnett, 
J. Chem. Phys. 6, 682 (1938). 

®§ J. Karweil and C. Schafer, Zeits. f. physik. Chemie B40, 
376 (1938). 
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(columns 1 and 2). The bands at 1216, 1538 and 
2012 cm have not been found by previous 
investigators. However, all these bands are 
weaker than the weakest ethane bands previously 
reported. They were found on every record at 
500 mm and their positions do not correspond 
to bands of any likely impurities. Column 3 lists 
the relative intensities of the bands, 4 presents 
the probable assignment, 5 the frequency calcu- 
lated for the combination or overtone bands 


TABLE IT. Low temperature spectrum of ethane. 











Frequency 
Wave-length Frequency difference 
mu cm" cm7 Intensity* 
13.6 735 87 vw 
12.22 818 . s 
12.11 826 
11.10 901 w 
10.15 985 vw 
8.3 1205 _ vw 
7.80 1282 88 w 
7.30 1370 s 
6.84 1462 76 s 
6.50 1538 , m 
5.67 1764 w 
5.50 1818 vw 
5.00 2000 vw 
4.52 2212 w 
4.25 2353 m 
3.85 2597 -- m 
3.74 2674 : 6 m-s 
3.64 2750 , m-s 
3.45 2900 s 
3.20 3125 m 
” come of symbols: s =strong, m =medium, w =weak, 7w =very 
weak. 


listed in 4. Columns 5 and 6 give the symmetry 
for the point groups D3, and D3a. 

Figure 3 presents the spectra of ethane 
observed at —170 and —195°C. In Table II 
the wave-lengths and frequencies are listed. 


DISCUSSION 


Stitt® has recently discussed in full the 
assignment of frequencies in ethane. He con- 
cluded that the available spectroscopic evidence 
provided no clue to the orientation of the 
methyl groups. Although three new ethane bands 
have been discovered in this work, no evidence 
has been adduced which alters Stitt’s conclusions. 
As is indicated in Table I, the bands at 1200 and 
2012 cm— may be assigned to either of the point 
groups D3, or Da. No assignment is apparent 
for the minimum observed in the gas at 1538 
cm~!, The 1400 cm band was found at 1414 by 
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Crawford, Avery, and Linnett using a spec- 
trometer of less resolving power. The value 
1400 cm-! agrees with the calculated frequency 
of the P branch of the 1380 cm band and 
leaves little doubt as to the origin of this mini- 
mum. Stitt suggested this assignment of the 
band. It is interesting to note that three assign- 
ments involve the torsional frequency. 





FREQUENCY IN CM 


PROPYLENE 


195 °C 


Fic. 4. Room temperature and Jow temperature spectra of propylene. 


Before discussing the data obtained at liquid 
nitrogen temperature it is well to consider 
briefly how the spectrum of a gas might be 
affected by a lowering of the temperature. 

(1) In a gas the band ‘“‘width’’** should vary 


** In parallel type bands the band width may be defined 
as the wave number separation of the P and R branches. In 
a perpendicular band of a symmetrical top molecule, 





PS RG TI EBL ER EGER BT ss Sa ie aia 


ESNet eter SRE 


tas Sea 


outer 


TY PRET ONS 


SO ane. a ot ob me Ute 








er 


be 


ry 


1ed 





ahi. me 


tied an ein im ea te AG 


{PRA EES PR BE aa Fat 








Af FPR TOS Cee Tee LIN SET tat 


as the square root of the absolute temperature 
as the temperature is lowered. When condensa- 
tion to a liquid takes place, we should expect the 
positions and shapes of the bands to be altered 
through intermolecular interactions and restric- 
tions on the rotation. In the solid, rotations 
should be further restricted if not prevented, 
producing a further change in the band shape. 

(2) In general, three translational and three 
rotational degrees of freedom of each molecule 
would be converted upon solidification to lattice 
vibrations. Since intermolecular forces in hvdro- 
carbon molecules would be expected to be 
considerably smaller than intramolecular forces, 
these vibrations would be of too low frequency 
to appear as fundamentals in the near infra-red 
even in the absence of restrictions imposed by 
the selection rules. However, combination bands 
would be possible. 

(3) As a rule we should not expect the 
symmetry of the solid lattice to conform to that 
of the gas, so that vibrations which are de- 
generate in the gas might not be so in the solid. 
Gas phase selection rules might also be altered 
by liquid or solid interactions. 

(4) Since overtones are possible only for 
anharmonic types of vibration and since inter- 
molecular interactions might affect the character 
of the potential functions, overtone bands 
obtained at low temperatures might be expected 
to show a change in intensity or position relative 
to the room temperature results. 

Comparison of Figs. 2 and 3 shows that most 
of the above anticipations are verified. The 
reduction in the width of the bands and the 
attendant increase in resolution are striking. The 
perpendicular 827 cm band is suitable for a 
semi-quantitative estimate of the temperature 


dependence. Considering the width at 60 percent 


transmission, this band is 141 cm~! wide at 
room temperature, 56 cm~! at —170°C (meas- 
ured as twice the wave numbers separation 
between the center and low frequency edge) and 
27 cm™ at —195°C. The predicted widths for 


maxima are observed only if the moment of inertia about 
the symmetry axis is not less than one-third of the other 
moment. However, the separation of points on the band 
where the intensity has a given value should show the same 
dependence. See S. L. Gerhard and D. M. Dennison, Phys. 
Rev. 43, 197 (1933). Similar considerations apply to bands 
of asymmetric top molecules. R. M. Badger and L. R. 
Zumwalt, J. Chem. Phys. 6, 711 (1938). 
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the gas would be 141, 83, and 80 cm~'. These 
measurements are, of course, very rough but 
the fact that the band widths are less than the 
predicted values suggests that there is some 
hindrance to molecular rotation in the liquid 
and considerable hindrance in the solid. The 
doublet structure of the band at —195°C 
suggests that the vibration is no longer degener- 
ate in the crystal lattice. This interpretation is 
strengthened by the lack of structure of the 
band in the liquid at —170°C. However, the 
doubling could also be ascribed to the appearance 
of a new band at —195°C. 

Comparing the spectra obtained at room 
temperature and liquid nitrogen temperature we 
note that all the major bands found at room 
temperature also appear in the liquid and solid 
spectra. The lowest degenerate fundamental is 
found at 822 cm~', 5 cm~ lower than the gas 
value. The A» fundamental appears at 1370 cm~' 
and the second degenerate fundamental at 1462 
cm~', the frequencies being lower than the gas 
phase values by 10 and 8 cm~!. The resolution 
of the spectrometer at 3000 cm™ is too small to 
permit conclusions to be drawn regarding any 


TABLE III. Partial vibrational spectrum of propylene. 








Room temperature —195°C 
Fre- Wave- Fre- Wave- Wells ; 
quency length quency length and Raman Calcu- 


cm7 mu cm”! mu Wilson = _ Effect lated 





Notinvestigated Notinvestigated 417 432 428 
Notinvestigated Notinvestigated 578 580 574 
Notinvestigated Notinvestigated 755 — — 

915 10.93 912 10.96 919 920 917 


940 10.64 936 10.68 — — 932 
os — 963 10.38 — — 1015 
991 10.09 994 10.06 ca. 996 —_ 976 
— — 1043 9.59 — — 1022 
1170 8.55 1172 8.53 1166 — 1170 
1206 8.30 _— — — — — 
— _— 1224 8.17 1224 
1236 8.09 — ~— — - —_ 
1250 8.00 1246 8.03 1244 - - 
1288 7.76 — — 1287 - 
1304 7.67 1297 7.71 — 1297) 1287 
1320 7.58 — — 1317 — — 
— — 1370 7.30 — — — 
1400 7.14 — <— 1399 — 1386 
— — — — 1416 1415 1429 
1444 6.93 1445 6.92 1448 1448 1471 
1466 6.82 —_ — 1472 _ 1474 
1489 
1516 6.60 1506 6.64 1508 — —_ 
1520 


1650 6.06 1648 6.07 1647 1648 1681 
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Fic. 5. Room temperature and low temperature spectra of 2-2 dimethyl butane. 


shifts in this region. The weaker bands all 
appear at —170°C and —195°C with the 
exception of the minimum found at 1400 cm~. 
This is to be expected from the assignment of 
this minimum to the P branch of the 1380 cm™ 
band, the lowering of the temperature reducing 
the band width so that the branches of the band 
are no longer resolved. 

Several bands appear at —170°C and —195°C 
which are not found at room temperature. In 
every case the bands are separated from strong 
bands by multiples of roughly 80 cm~!. Such a 
low frequency cannot correspond to a funda- 
mental vibration of the ethane molecule and 
must, therefore, be ascribed to a lattice vibration, 
possibly a restricted rotation of the ethane 
molecule in the crystal or liquid. 

Some differences in the intensities and posi- 
tions of the overtone bands have been observed, 
but in view of the large scattering correction 
mentioned above we reserve comment until 


more data have been obtained. 





PROPYLENE 


The study of the spectrum of ethane disclosed 
that no significant shifts in the positions of the 
bands occurred when the temperature was 
lowered to —195°C. Since it is well known that 
the spectra of strongly polar compounds are 
affected by the change from gaseous to liquid 
state it was interesting to investigate the 
temperature shift in a more polar hydrocarbon. 
Propylene was chosen for study because, in 
addition to being somewhat polar it has been 
the subject of a recent theoretical analysis by 
Wells and Wilson.* They found that the over- 
lapping of the bands at room temperature 
prevented a complete assignment of the fre- 
quencies. Thus, a low temperature measurement 
of the compound might be of interest from two 
standpoints. 

The region 6 to 12.5 mu was chosen for 
examination; for a ready comparison with the 


6 A. J. Wells and E. Bright Wilson, Jr., J. Chem. Phys. 
9, 319 (1941). 
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data of Wells and Wilson, transmissions have 
been plotted as a function of frequency rather 
than of wave-length. The curves are presented 
in Fig. 4. Numbers on the curves indicate 
pressures in mm of Hg. The low temperature 
cell was used for all measurements. 

In Table III are presented the frequencies and 
wave-lengths of bands observed in propylene at 
room temperature and —195°C (columns 1, 2, 
3 and 4). In column 5 the results of Wells and 
Wilson are presented. Column 6 lists the fre- 
quencies observed in the Raman effect and 


TABLE IV. Wave-lengths and frequencies of significant 
bands in 2, 2-dimethyi butane. 





Room temperature —195°C Raman effect* 





Wave-length Frequency Wave-length Frequency Frequency 
mu cm"! mu em" em! 
259 (1) 
335 (0) 
361 (1) 
24.4 410 
23.9 418 
23.7 422 
20.64 484 485 (1) 
14.3 700 
14.09 710 14.07 711 711 (10-b) 
13.91 719 
12.9 775 
12.79 782 12.81 781 
12.68 789 
11.48 871 11.50 870 869 (3) 
10.75 930 10.74 931 927 (3-b) 
10.0 1000 10.14 986 
9.82 1018 10.03 997 1017 (2-b) 
9.78 1022 
9.32 1073 9.31 1074 1075 (2) 
8.22 1216 8.23 1215 1217 (3-d) 
8.0 1250 7.99 1251 1251 (2-b) 
7.65 1307 7.65 1307 1304 (1-d) 
7.45 1340 
7.27 1376 7.33 1364 
7.17 1395 
6.95 1439 1447 (6) 
6.78 1475 6.82 1466 1468 (6) 
6.75 1480 
6.42 1558 6.25 1600 
6.32 1582 
5.99 1669 6.00 1667 
5.77 1733 5.86 1706 
5.22 1916 5.22 1916 
4.88 2049 4.80 2083 
4.30 2326 4.32 2315 
37 2700 2720 (1) 
2760 (1) 
2860 (3) 
3.45 2899 3.47 2882 2874 (10) 
2903 (3-d) 
2936 (3-d) 
2965 (6-d) 
2980 (3-d) 
3.19 3135 
2.37 4219 2.37 4219 








_ *E. J. Rosenbaum, A. V. Grosse and H. F. Jacobson, J. Am. Chem. 
Soc. 61, 689 (1939). Note: Numbers in parentheses refer to intensity; 
b means broad, d diffuse. 
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column 7 the calculated positions of the absorp- 
tion bands. A comparison of column 3 with 
column 5 shows that three bands were found in 
this research which were not reported by Wilson 
and Wells, i.e., 936, 963 and 1043 cm-. It is 
interesting to observe that the _ theoretical 
analysis of the spectrum of propylene predicts 
four bands, which Wilson and Wells, on the 
basis of their data, were unable to correlate 
with observed frequencies. The unknown fre- 
quencies were computed by them to be at 
about 920, 1015, 1020 and 1100 cm-. There 
seems to be little doubt that the bands brought 
out by the low temperature experiments corre- 
spond to three of the predicted frequencies. 

The frequency shifts associated with lowering 
of the temperature are seen to be small with the 
exception of the band at 1370 cm~'. The explana- 
tion of this displacement is not clear but it may 
be related to another phenomenon which is 
produced by a lowering of the temperature: a 
variation in the intensity of certain bands. Thus, 
a preliminary series of curves obtained at 
temperature intervals of about 5°C from —195 
to —170°C in the spectral region 900 to 1050 
cm! show a pronounced temperature dependence 
of the relative intensities of the bands. The 
disappearance of the band at 1466 cm~ and the 
apparent shift in the position of the 1400 cm~! 
band may well be other examples of this temper- 
ature dependence, but a definite explanation 
requires a more thorough investigation than was 
possible in this survey. Perhaps the minima at 
1466 and 1444 cm are P and R branches of a 
single band whose center is shifted to 1445 cm! 
at —195°C. 


2, 2-DIMETHYL BUTANE 


This compound was examined because it is 
fairly representative of the higher hydrocarbons 
and because the liquid may be supercooled to 


TABLE V. Comparison of widths of bands at room tempera- 
ture and —195°C. 








Frequency of 





band center Room 
(cm~!) temperature —195°C Ratio Transmission 
782 25 13 0.52 60 
930 24 12 0.50 80 
1073 24 12 0.50 60 
1216 56 28 0.50 60 

















18 CLEVELAND, MURRAY, 
—195°C without crystallization. The spectrum 
is presented in Fig. 5, the wave-lengths and 
frequencies of the bands at the two temperatures 
are given in Table IV in columns 1 through 4. 
Column 5 presents the Raman spectrum obtained 
by Rosenbaum, Grosse, and Jacobson. 

The increase in resolution obtained by the use 
of low temperatures is particularly evident in 
this spectrum. It is interesting to note that the 
group of bands of 10 mu appear in the Raman 
effect only as a single broad line. The positions 
of the bands are affected only slightly if at all 
by the temperature change, except in the 
overtone region, and the agreement between the 
Raman and infra-red frequencies is excellent. 
As is shown in Table V, the reduction in the 
width of the bands at — 195°C is in close agree- 
ment with the theoretical value. Column 4 gives 
the ratio between the width at room temperature 
and that at —195°C. The ratio predicted for a 
gas is 0.51. Column 5 gives the value of the 
percentage transmission at which the width of 
the band was determined. 

The close agreement between calculated and 
observed band widths and gas and _ liquid 
frequencies suggests that the vibrational and 
rotational partition functions for 2, 2-dimethyl 


COLEY, 





AND KOMAREWSKY 

butane are unaffected by the change from the 
gaseous to the liquid state and that the entropy 
change on vaporization could be accurately 
correlated with the change in free volume. 


CONCLUSION 


Although the experiments described above are 
intended only to serve as an introduction to low 
temperature spectroscopy, they indicate that 
the method promises to be very valuable in 
several types of problems. In the spectroscopic 
analysis of molecular structure, it provides a 
better basis for comparison of Raman and 
infra-red data and a means of distinguishing 
bands unresolved by conventional methods. It 
provides a very useful tool for the study of the 
liquid state. The range of analytical applications 
of infra-red spectroscopy should be greatly 
enlarged by the use of low temperature technique. 
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Displacements, estimated intensities and depolarization factors have been measured for 
acetone, 2,6-dimethyl-4-heptanone, 5-nonanone, 7-tridecanone, 8-pentadecanone and 5- 
tetradecanone. The results for acetone are compared with the previous data: Eight lines 
previously reported by one or more investigators are believed to be non-existent. A new weak 
line was observed at 1676 cm~. The depolarization factors obtained agree well with those 
reported by Simons, except that the 1356 line was found to be considerably less polarized and 


data for five additional lines were obtained. 


INTRODUCTION 
HE dialkyl ketones, 2,6-dimethyl-4-hepta- 
none, 5-nonanone, 7-tridecanone, 8-penta- 
decanone and 5-tetradecanone were prepared in 
another investigation by a new and convenient 





method of ketone production and it seemed 





desirable to obtain their Raman spectra while d 
they were available. In searching the literature, P 
it was noted that there were some discrepancies 

S| 


in the Raman data reported by many investi- 
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TABLE I. Raman spectrum of acetone as obtained by several investigators. 


RAMAN SPECTRA OF 


DIALKYL KETONES 














Mean values Kohlrausch K6ppl Sirkar Present results Probable values 
N Av I p Av I v Av I p Av I p 
1 360 0 
8 389 2 0.87 391 1 392 1.5 390 2 0.9 391 2 0.87 
7 489 1 486 i 487 3 492 1 0.9 488 1 0.9 
12 528 2 0.70 532 3b 527 2 531 4 0.7 530 3 0.70 
1 544 1 
1 594 00 594 00 580 0 587 0 
1 727 735 0 731 0 
15 788 7 0.20 787 8 789 5 + 786 8 0.1 787 7 0.20 
1 821 
5 903 5 902 00 903 0 899 } 901 } 
3 930 1 
12 1068 2 0.64 1066 2 1065 1 1066 4 P 1066 2 0.64 
1 1093 0 
2 1192 3 0.88 
12 4 3 0.88 1220 2 1221 1 1221 5 0.8 1221 3 0.88 
1 124 
9 1350 1 0.23 1354 } 1358 0 1355 2 0.6 1356 1 0.6 
1423 3 o 1423 3 ) 
11 1429 3h 0.90 1438 4 1423 3h 1427 6b 0.8 1438 af 0.90 
1 1489 
1676 0 1676 0 
12 1711 4 0.38 1706 5) 1705 5 1710 7 0.5 1710 6 0.38 
1 1748] 1 1751 ! 1750 i 
3 2692 1 2689 00 2689 0 2697 1 P 2692 1 P 
5 2849 1 2838 1 2857 Ob 2843 3 2840 3 ) 
» 2p ms 
1 2879 2865 2) 2871 2} 
14 2924 7 0.22 2922 12vb 2923 10 2921 10 0.1 2922 10 0.22 
6 2967 4 2965 5 2964 4 2963 4 0.9 2965 4 0.9 
8 3005 4 3005 4d 3006 4 3004 5 0.8 3005 4 0.8 











Av=displacement in cm~, J=estimated intensity, p=depolarization factor, + =Stokes and anti-Stokes, brackets 
enclose data in regard to which there is some uncertainty, the brace joins lines that were unresolved on the polarization 
spectrogram, b=broader than usual, v) = very much broader than usual, P= somewhat polarized but a quantitative esti- 
mate of the depolarization factor was not possible, N=number of times the line has been reported by independent 


investigators. 


gators for acetone, the simplest dialkyl ketone, 
and that the depolarization factors for several 
strong lines had not been reported. Consequently, 
the Raman spectrum of this compound was 
obtained from carefully purified samples. 


EXPERIMENTAL 


The apparatus and experimental technique 
were essentially the same as described in previous 


papers.! A few changes have been made: The 


mercury arc now has tungsten-Pyrex seals at 
the electrodes and was baked out for two hours 
at 450°C, under high vacuum, before the mercury 
was distilled into it. The side tube, which now 
has a somewhat larger diameter than previously 
in order to increase its effectiveness as a con- 
denser was then sealed off while still on the 
pumps. The arc thus prepared has been in 

1F. F. Cleveland, M. J. Murray, H. H. Haney, and J. 


Shackelford, J. Chem. Phys. 8, 153 (1940); Forrest F. 
Cleveland and M. J. Murray, J. Chem. Phys. 7, 396 (1939). 


practically continuous operation for six months 
and shows no sign of leakage or deterioration. 
The arc is now started without changing the 
series resistance by placing the tip of a high 
frequency coil (commonly used for testing 
vacuum assemblages for leakage) near the 
mercury surface at the negative terminal.* 

In the depolarization measurements, the 
Polaroid disk has been replaced by a Nicol 
prism, on account of the rather strong absorption 
of light by the Polaroid. 

The ketones which were solids at room temper- 
ature were kept in the liquid state during expo- 
sure by means of an electric heating clement 
wrapped loosely about the glass tube containing 
the sample. The width of the slit of the spectro- 
graph was set at 0.08 mm for the regular expo- 


2We are indebted to Professor Paul L. Copeland, 
Department of Physics, for assistance in preparing the 
mercury arc and for suggesting this convenient method of 
starting it. 
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TABLE IT. Raman spectra of some dialkyl ketones. 
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Desig- 2,6-dimethyl-4-heptanone 5-nonanone 7-tridecanone 8-pentadecanone 5-tetradecanone Desig- 
nation v I p Av I p Av I p Av Av I p nation 
224 2b 
+262 4 0.7 2572.) 
- 0.6 
295 2 271, — 2 J 
a +349 3b 394 1 
b +419 4 0.6 
533 2b 
715 1 715 2b 
c 779 3b 0.3 772 1 776 1 
815 2 823 2b 0.6 819 3b ) c 
d 840 60.3 8431 8421 835 2b) 
870 1 874 3 | 863 2 872 4 872 o Ff d 
- 0.6 
8922) 889 3 0.6 893 2 893 3b P ° 
929 1 
e 962 5 0.8 970 1 0.9 962 2 
994 b 
f 10391 1020 1b 0.9 f 
1059 1 1056 5 ) 1062 5 0.6 
> 0.7 § 
1074 1 1065 3 } 1074 4b 0.7 1074 5d 1083 5 0.6 
h 1133 5b 0.7 1121 4 0.5 1119 4 0.6 1120 5 1123 5 0.4 h 
Z 1171 4 0.9 1168 1 1167 1 i 
1215 1 
1301 4. 0.7 1308 5 0.7 1303 6 1301 7 «0.6 k 
l 1333 4 0.8 
1411 3 0.7 1413 3. 0.7 1411 2 G7 1411 1 1405 3 0.9 
1449 8 | 1452 7 ) 1441 7 | 1438 8 1441 8 | 
m , 0.8 ¢ 0.7 > 0.7 > 0.7 m 
1465 8, 14624] 1456 5 | 1455 6 1455 8) 
1714 6 0.6 1717 5 0.6 1717 3. 0.4 1713 3 1718 4 0.7 
n 2719 , + 
2736 2 2731 1 0.6 2731 2 2732 2 @7 n 
2762 2 
2858 10 0.1 2852 10 2853 10 0.1 7) 
7) 2870 10 0.3 2869 10d 0.2 2884 10 | [2865 6] 2870 =10) 
- 0.4 0.6 
p 2004 9 06 2911 9 06 2907 10) 2004 «28 = 2897-10 p 
2930 5 P 2936 10 0.2 2933 8 0.1 2931 9 2933 10 = 0.1 q 
r 2966 10 0.9 2965 5 0.9 2965 4 0.9 2960 6 2967 5 0.9 r 
2983 3 
C.M.* C.W. C.W. C.W. C.M. 








sures and 0.12 mm for the polarization exposures. 
The time of the exposures was 7-10 hours. 
Excitation was by Hg 4358A. 

The ketones used in this investigation, with 
the exception of acetone and 2,6-dimethyl-4- 
heptanone which were both Carbide and Carbon 
Chemical Corporation products, were prepared 
by a catalytic method,’ involving the simultane- 
ous dehydrogenation and condensation of pri- 
mary alcohols. The liquid ketones were purified 
by distillation in a column equivalent to twelve 
theoretical plates. The solid ketones were purified 
by recrystallization and finally were melted and 


3V. I. Komarewsky and J. R. Coley, J. Am. Chem. Soc. 
63, 700 (1941). 


« C.M.=continuum medium, C.W.=continuum weak; bold type indicates previously observed lines. 





distilled in vacuum directly into the Raman 
tube. The boiling or melting points were as 
follows: Acetone, b. p. 56.1°C, 2,6-dimethyl-4- 
heptanone, b. p. 168°, 5-nonanone, b. p. 187- 
188°, 7-tridecanone, m. p. 32.5-33.0°, 8-penta- 
decanone, m. p. 39-40°, 5-tetradecanone, m. p. 
25.5-26.0°. Two samples of acetone were used, 
one of which was purified by use of the sodium 
bisulfite addition compound. 


RAMAN SPECTRUM OF ACETONE 


The Raman data for acetone, as obtained by 
many investigators‘ are summarized in the first 


4 See bibliography in Der Smekal-Raman Effekt, Ergans- 
ungsband (Julius Springer, Berlin, 1938), p. 220. 
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columns of Table I. The depolarization factors 
are those measured by Simons.*® In several of the 
earlier investigations unfiltered light was used. 
This made it difficult to decide which of the 
mercury lines was the exciting line and fre- 
quently non-existent lines were reported. With 
this in mind, the later work of Kohlrausch and 
KX6ppl® and of Sirkar’ are listed separately in 
Table I. Finally, the present results and the 
values which seem to the authors to be the most 
probable ones are listed in the columns at the 
right. Frequencies that have been observed by 
only one of the investigators must be regarded 
as questionable and have for the most part not 
been included in the list of probable values. 

Kohlrausch and Pongratz® have pointed out 
that the line 1192 is not real, but corresponds to 
the 2965 cm line excited by Hg 4047A. This 
conclusion is supported by the present results 
in which, not only was the 1192 line absent, but 
also the depolarization factor of the 2965 line 
(previously unknown) was found to be the same, 
within experimental error, as that of the line 
measured by Simons as 1192. The line 930, 
observed by three investigators, probably corre- 
sponds to the 2692 line excited by Hg 4047A. 

The polarization data agree quite well with 
those of Simons, except that the 1356 line was 
found to be considerably less polarized and data 
for five additional lines were obtained. 


RAMAN SPECTRUM OF 2,6-DIMETHYL- 
4~-HEPTANONE 


To facilitate comparison, the spectra of the 
other five dialkyl ketones are listed side by side 
in Table II. Of these, the spectrum was most 
complete for the 2,6-dimethyl-4-heptanone. Dis- 
placements and estimated intensities have been 


reported for this compound by Kohlrausch and 


K6ppl.* Displacements which agree within 
experimental error with their results appear in 
bold type in Table II. The weak lines at 1276 


5 L. Simons, Soc. Scient. Fenn. Com. Phys. Math. 6, No. 
13 (1932). 

6K. W. F. Kohlrausch and F. Képpl, Zeits. f. physik. 
Chemie B24, 370 (1934). 

7S. C. Sirkar, Ind. J. Phys. 7, 61 (1932). 

8 K. W. F. Kohlrausch and A. Pongratz, Zeits. f. physik. 
Chemie B27, 180 (1934). 





and 1619 cm™ were not observed in the present 
investigation, but nine new lines were found. 
The doublet 1449, 1465 corresponds to the broad 
line at 1456 cm™ in the previous report. 

The lines designated by the letters a, 0b, etc. 
are those listed by Kohlrausch and Képpl® as 
characteristic of isobutyl compounds. They were 
obtained by examination of the spectra of eight 
compounds containing the (CH3;)2»CHCH.-group. 
The surprisingly large intensities of the 2870 
and 2966 lines may result from coincidence of 
CH; and CHg frequencies, since acetone, con- 
taining only CH; groups, had frequencies at 
2871 and 2965 cm™. 


KETONES OF THE TYPE 
CH;(CH:),CO(CH,),,CH; 


In general, the spectra of the four ketones of 
this type are not as complete as for the one just 
discussed. The frequencies given by Kohlrausch 
and Képpl® as characteristic of C,Hen+1 are 
indicated by the letters in the right-hand 
column of Table II. A new frequency, not 
observed for acetone or 2,6-dimethyl-4-hepta- 
none, appears near 1303 cm. It is somewhat 
polarized and seems to become more intense as 
the length of the chain increases. The possibility 
that this may be a chain frequency is suggested 
by the fact that Whitcomb, Nielsen and Thomas’ 
have calculated values of 385, 1106, 1354, 1463, 
2320 and 2778 cm for the chain frequencies of 
straight chain paraffins. Values observed by 
them in the infra-red spectrum of undecane were 
1066, 1308, 1463, 2320 and 2778 cm™, the 
exact coincidence of the last three observed 
values with the calculated values being explained 
by the fact that these frequencies were used to 
adjust the force constants. Frequencies observed 
here for the normal ketones which may corre- 
spond to these are 394 (observed for only one), 
1075, 1303 and 2733 cm™. 

In the 3000 cm region an exceedingly 
intense, polarized line appears at 2855 cm for 
the last three ketones and the 2983 line, found 
for the 2,6-dimethyl-4-heptanone, was not ob- 
served for any of the four normal ketones. 


9S. E. Whitcomb, H. H. Nielsen, and L. H. Thomas, J. 
Chem. Phys. 8, 143 (1940). 
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An investigation has been made of the reactions of ethylene photosensitized by 5'P; cadmium 
atoms, i.e., using the resonance line at 2288A. It is found that reaction occurs readily, in 
contrast to the action of 5°P,; atoms, which have been previously shown to be ineffective. The 
products of the reaction are butene, hexene, and higher olefins, together with smaller amounts 
of acetylene and hydrogen. No pressure rise is detectable in the early stages of the reaction, and 
acetylene is a much less important product than is the case with mercury photosensitization. It 
is concluded that the primary step may be either 


Cd (P)) +C2H, = Cd (So) +C.H 3+ H 


or 


Cd(#P1)+C2Hs=Cd('So) +C2Hs*. 


Subsequent polymerization probably occurs by a free radical mechanism. 





N recent years we have made a number of 
investigations of reactions photosensitized by 
53P; cadmium atoms, the atoms being excited 
by absorption of the resonance line at 3261A.!~4 
These results are of interest by comparison with 
mercury photosensitization, since the magnitude 
of the excitation energy of the photosensitizer 
is much smaller, viz., 87 kcal. as compared with 
112 kcal. for Hg(6*P;). In view of this it seemed 
worth while to extend the investigations to a 
number of excitation energies of different magni- 
tudes, and experiments have accordingly been 
made with Cd(5'P;), Zn(4*°P1), and Zn(4'P;) 
atoms. The results with zinc will appear shortly.® 
The present paper deals with photosensitization 
by Cd(5'P;), the atoms being excited by absorp- 
tion of the lower resonance line at 2288A. The 
excitation energies, etc., corresponding to the 
use of both resonance lines for zinc, cadmium, 
and mercury are given in Table I.° 
In view of the interest attaching to the 
photosensitized polymerization of ethylene, we 
have chosen this reaction as the initial investi- 
gation of photosensitization with Cd(5!P;) atoms. 


1E. W. R. Steacie and R. Potvin, Can. J. Research B16, 
337 (1938). . 
039) W. R. Steacie and R. Potvin, J. Chem. Phys. 7, 782 

1 A 

3 E. W. R. Steacie and R. Potvin, Can. J. Research B18, 
47 (1940). 

4E. W. R. Steacie, D. J. LeRoy and R. Potvin, J. Chem. 
Phys. 9, 306 (1941). 

5 E. W. R. Steacie and H. Habeeb, unpublished. 

6 E. W. R. Steacie, Ann. N. Y. Acad. Sci. 41, 187 (1941). 
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A number of investigations of the mercury 
photosensitized reactions of ethylene have been 
made, largely by Taylor and his collaborators. 
Virtually all investigators adopt for the primary 
step the suggestion of Taylor and Bates? * 


C.H4y+Heg(*P3) = C.H2+H2+Heg('So). 


It is then assumed that the subsequent reactions 
are the photosensitized polymerization of acety- 
lene, together with reactions of hydrogen atoms 
and of radicals with ethylene, viz., 


H +C.Hy= C.H; 
Co.H;+CoHy= C4Hg, etc. 


The main evidence for this mechanism is the 
detection of large amounts of acetylene in the 
products of the reaction, and the fact that the 
pressure rises considerably in the early stages 
of the reaction, passes through a maximum, and 
then decreases. 

We have recently? made a more detailed 
investigation of the reaction with special empha- 
sis on the early stages. It was found that in the 
initial stages the pressure rise and acetylene 
production were exactly equivalent, but that the 
initial rate of acetylene production (or of 
pressure increase) decreased with increasing 


7H. S. Taylor and J. R. Bates, Proc. Nat. Acad. Sci. 
(Washington) 12, 714 (1926). 

8 J. R. Bates and H. S. Taylor, J. Am. Chem. Soc. 49, 
2438 (1927). 

*D. J. LeRoy and E. W. R. Steacie, J. Chem. Phys. 9, 
829 (1941). 
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TABLE I. Excitation energies, etc. 








‘Maximum energy 
available”’ to 
split off H atom, 





Energy of Heat of assuming inter- 
excited dissociation mediate hydride 
Resonance atom, of hydride, formation, 

Substance line, A Transition kcal. kcal. keal. 
Mercury 1849 61S )—6'P, 153.9 8.5 162.4 
2537 6'So— 6P 112.2 8.5 120.7 
Cadmium 2288 5'1S9—5'P1 124.4 15.5 139.9 
3261 5'1S9— 5°P 87.3 3.5 102.8 
Zinc 2139 41S9—4'P, 133.4 23.1 156.5 
3076 41S,—43P, 92.5 23.1 115.6 


pressure. This points definitely to an excited 
molecule type of mechanism, viz., 


CoH4+Heg(@P1) = Hg (So) +C2H4* 
CoH4*=CoH2+He 
CoH4+Ce2H a* = 2C 3H. 


One of the most striking features of the 
investigation was the explanation afforded for 
the subsequent pressure decrease. It has usually 
been assumed that the decrease was caused by 
the photosensitized polymerization of acetylene 
formed in the early stages of the reaction. It 
was found that this was definitely contrary to 
the facts, since the acetylene concentration 
increased without interruption after the pressure 
had begun to fall. It was found that hydrogen 
was responsible for the pressure decrease, 
presumably by an atom and radical sensitized 
polymerization of ethylene, viz., 


H.+He(*P;) = Hg('So) +2H 
H +C.H, = C.H; 
2C.H;= C,H 10 
C.H4+CoH; = C,4Hg, etc. 


Acetylene proved to be surprisingly stable to 
attack by excited mercury atoms in the presence 
of ethylene, its concentration in some cases 
rising to 30 percent of the total gas mixture. 

Bates and Taylor” made the first investigation 
of the cadmium photosensitized reactions of 
ethylene-hydrogen mixtures, using the resonance 
line at 32614, i.e., 53P; cadmium atoms. A more 
detailed investigation of the reactions photo- 
sensitized by Cd(5*P,) atoms has recently been 
made by one of us.'! The results are striking in 


107. R. Bates and H. S. Taylor, J. Am. Chem. Soc. 50, 
771 (1928). 

NE. W. R. Steacie and R. Potvin, Can. J. Research B16, 
337 (1938). 


that they indicate a quantum efficiency of about 
0.5 for a 1 : 1 ethylene-hydrogen mixture, but a 
yield of only about 0.015 for pure ethylene. 
Furthermore, this does not seem to be due to a 
lack of quenching by ethylene, since qualitative 
indications are that ethylene quenches reason- 
ably strongly. 

Two possible explanations suggest themselves 
for the difference between the results with 
Cd(?P;) and Hg(*P;) atoms. The reaction 


Co.H,=C2H2+He (1) 


is endothermic to the extent of 42.9 kcal.” The 
activation energy of the reaction has not been 
determined, but Taylor and van Hook" have 
found that of the reverse reaction to be roughly 
42 kcal. Combining these results we obtain a 
value of approximately 85 kcal. for the activation 
energy of reaction (1). It is therefore possible 
that the 112 kcal. associated with a *P; mercury 
atom can cause reaction, while the 87-kcal. 
excitation energy of a *P; cadmium atom cannot 
be used efficiently enough to do so." 

On the other hand, there is a great deal of 
evidence that the polymerization of ethylene 
can be induced by free radicals, and this suggests 
the possibility of the thermal reaction proceeding 
by a free radical mechanism.'® It also suggests 
the possibility of the photosensitized reaction 


12 F.R. Bichowsky and F. D. Rossini, Thermochemistry of 
Chemical Substances (New York, 1936). 

13H.A. Taylor and A. van Hook, J. Phys. Chem. 39, 811 
(1935). 

4 For a discussion of this point and of the mechanism of 
quenching from the standpoint of potential energy surfaces 
and the theory of absolute reaction rates see K. J. Laidler, 
J. Chem. Phys. 10, 34 (1942). 

16 See, e.g., H. D. Burnham and R. N. Pease, J. Am. 
Chem. Soc. 62, 453 (1940). 
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involving a primary split into radicals, i.e., 


C,H,;+M*=C.H3;+H+M 
or =C2H3+MH. 


If this is the case, the Cd(*P:1) photosensitized 
reaction would be ruled out if the strength of 
the C—H bond in ethylene is greater than 103 
kcal. It is, however, impossible to reconcile this 
with the facts as far as the mercury photo- 
sensitized reaction is concerned. However, it 
seemed possible that such a mechanism might 
hold for other sensitizers, and accordingly we 
have investigated the polymerization of ethylene 
photosensitized by Cd(5!P,) atoms. 


EXPERIMENTAL 
The Light Source 


The cadmium resonance lamp was similar to 
one previously described’? except that the 
emitting portion of the lamp was made of quartz 
rather than of Corex D glass, in order to transmit 
\2288. The lamp consisted of a 6-inch length of 
11 mm outside diameter quartz tubing connected 
by graded seals to Pyrex tubing of the same 
diameter. To the ends of this were sealed 
electrode chambers of Pyrex, 1 inch in diameter 
and 3 inches long. These contained standard 


Fic. 1. Spectro- 
gram of cadmium 
resonance lamp. 
Successive expo- 
sures 60, 10, 2, 1, 
3, and 3 sec. 
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neon sign ‘‘coated”’ sheet iron electrodes in the 
form of hollow cylinders } inch in diameter and 
13 inches long. The ends of the electrodes 
contained annular porcelain insulators to prevent 
excessive sputtering. The lamp contained a piece 
of cadmium, weighing about 0.1 gram, and was 
filled with neon at a pressure of 3 mm at room 
temperature. Mercury vapor was excluded during 
the evacuation and filling. The filling of the 
tube, outgassing, etc., was done by the usual 
neon sign technique. 

The lamp was operated at 120 milliamperes 
from a 6000-volt sign transformer. 

A set of spectrograms of the lamp, taken with 
a Hilger intermediate quartz spectrograph, is 
shown in Fig. 1. Considering the diminished 
plate sensitivity in the short ultraviolet it is 
apparent that the two resonance lines, 2288A, 
and 3261A, are much more intense than any 
other lines. The zinc resonance line 2139A, and 
the mercury line at 2537A appear, but their 
intensities are far too low to have any appreciable 
effect. The zinc line at 3076A is probably 
present also, but is masked by the fairly strong 
cadmium line, 3082A. It is evident that the 
lamp is highly efficient for photosensitization 
purposes, since the intensity is mainly in the two 
resonance lines. In this respect it compares 
very favorably with commercial high voltage 
mercury lamps. 

When desired the lower resonance line could 
be completely filtered out, and 3261 almost 
completely transmitted by surrounding the lamp 
with a filter consisting of an 8-inch length of 
Corex D tubing with a wall thickness of about 
2 mm. 

The output of the lamp was determined only 
very roughly, since it is intended to make 
quantum yield determinations of this and other 
cadmium photosensitized reactions in the near 
future. There appears to be no suitable acti- 
nometer for 2288A which could be used at 300°C. 
The intensity was therefore merely estimated 
geometrically from the amount emerging from a 
quartz window on the side of the furnace. This 
indicated that the amount of \2288 entering the 
reaction vessel was of the order of 10~* to 10-7 
einstein per sec. 

The output was very constant, and a series of 
measurements with a thermopile, and with a 








wi 


th 


Wi 


ac 
re 
m. 


co 
‘oe 
th 


n 


nD = 


Cn O kr OMS OO mes 


preteens 





hs brat 


2 8 as 


PHOTOSENSITIZED 


} CIRCULATING 








=o 


! 
uL 


POLYMERIZATION OF C.H, 25 


PUMP ; 





AND 
MERCURY 





VALVES 


























Fic. 2. Apparatus. 


General Electric F.J. 405 sodium photo-tube 
indicated that the output was constant to within 
5 percent over long periods of time. 

More accurate measurements were made of 
the relative intensities of \2288 and \3261. In 
these a beam of radiation from the quartz 
window was defined by stops (but without 
lenses, to avoid focal isolation) and then used in 
various ways: 

(a) A uranyl oxalate actinometer was used 
with and without a chlorine filter. Allowing for 
the transmission of the filter this gave approxi- 
mately \3261/A2288 = 1.5. 

(b) A sodium photo-tube was used together 
with various combinations of chlorine and Corex 
D filters. Allowing for the transmission, and for 
the variation of the sensitivity of the photo-tube 
with wave-length, this gave in 4 determinations 


3261/2288 = 0.75, 0.97, 0.97, 1.04. 


This value of the ratio will be somewhat low on 
account of the presence of other lines in the 
region of 2800A where the photo-tube has 
maximum sensitivity. 

(c) A thermopile was used with various filter 
combinations. This gave a value for the ratio of 
2.3, which will be somewhat high since some of 
the long ultraviolet will appear as \3261. 

(d) Photographic determinations were made 
using Eastman 40 U. V. sensitized plates. With 
rough corrections for the falling off in plate 
sensitivity at the shorter wave-lengths these 
indicated a value of about 2 for the ratio. 


~e 


It thus appears that the ratio of the intensities 
of \3261 to \2288 is not far from 1.5. 

For the interpretation of the results of photo- 
sensitization experiments, some knowledge of 
quenching efficiencies is required. An investiga- 
tion of the quenching of cadmium resonance 
radiation by a number of gases is in progress, 
and pending the completion of this work only 
rough “absorption’’ measurements were made. 
For these cadmium turnings were placed directly 
in the reaction cell, and the intensity of the 
resonance radiation transmitted at 225°C was 
measured photographically, and also in some 
cases with the photo-tube or thermopile and 
suitable filter combinations. These approximate 
measurements indicated a rather low absorption 
of \3261 by ethylene, i.e., witha 4-cm path 
there was about 30 percent absorption at 35-mm 
pressure, while under the same circumstances 
over 90 percent of 42288 was absorbed. 


The Reaction System 


The reaction system is shown diagrammati- 
cally in Fig. 2. The lamp was contained in a 
furnace A, which was cylindrical in shape and 
was heated electrically. The lamp was surrounded 
by an annular fused quartz reaction vessel, 
10 cm long, 8 cm in diameter, with a 4-cm 
diameter hole through the middle. The tube 
leading from the right-hand end of the reaction 
vessel passed into a second furnace B, and was 
connected through a graded quartz to Pyrex seal 
to a Pyrex tube 1 inch in diameter and 1 foot 
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Fic. 3. Pressure decrease—time curves for ethylene alone. 
Ethylene pressures in mm as shown. The lowest curve 
shows the rate with the filter on, i.e., with 43261 only. 


long which contained cadmium turnings. This 
acted as a saturator to furnish cadmium vapor 
to the reaction system. The temperatures of A 
and B could be independently varied. The 
furnace A contained a fused quartz window 
opposite the reaction vessel through which 
spectroscopic observations of the lamp could be 
made. 

The temperatures of the two furnaces were 
measured by means of copper-constantan thermo- 
couples. The temperatures were regulated by 
Cenco bimetallic regulators. We have found that 
these can be lengthened to over three feet 
without affecting their operation. This is done 
by merely sawing the outer tube and inner rod 
in half, and brazing on an additional length of 
brass tube and of drill rod with silver solder. 

The experimental arrangement could be used 
for static, circulatory, or single-pass flow runs. 
For circulatory runs a circulating pump was 
used which consisted essentially of two Toepler 
pumps with mercury reservoirs on a_ rocker 
which alternately raised and lowered them. This 
arrangement avoided volume, and _ pressure, 
changes in the system during the cycle of 
operation of the pumps. The pumps had a speed 
of about 500 cc per minute. A system of mercury 
bubblers was used as valves, so that the gas 
always circulated in the same direction, i.e., 
through the saturator and then the furnace. 
After leaving the furnace the gases passed 
through a trap C. This could be maintained at 
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any desired low temperature by an arrangement 
which has been previously described.’ After 
passing the trap the gases returned to the 
circulating pump. 

A manometer was connected to the system at 
M. By a suitable manipulation of stopcocks the 
trap C could be isolated or by-passed when 
desired. The usual oil pump, mercury vapor 
pump, McLeod gauge, gas reservoirs, etc., were 
connected at D. Gaseous or liquid products 
could be removed through E. The total apparent 
volume of the circulating system, when the 
furnaces were at 300°C, was 945 cc. 

Static runs could be made when desired by 
sealing off the tubing between the reaction 
vessel and the saturator at the point F, putting 
cadmium turnings in the reaction vessel, and 
shutting the stopcocks leading to the circulating 
pump. The apparent volume of the system when 
used in this way was 512 cc. 

To make single-pass runs stopcocks G and [I 
were closed. A cylinder of ethylene with a 
regulating valve set at about 5 lb. was connected 
at I, stopcock J was opened, and gas flowed 
into the apparatus at a rate determined by the 
size of the capillary K. The gas passed through 
the saturator, reaction vessel, and the trap C, 
and then through the file-scratched stopcock L 
to a liquid-air trap N. The filed stopcock was 
adjusted until the desired pressure was indicated 
on the manometer. In the absence of permanent 
gases all the gas flowing through the apparatus 
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Fic. 4. Pressure decrease—time curves for ethylene- 
hydrogen mixtures. Ethylene pressure approximately 200 
mm. Hydrogen pressures in mm as shown, 
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was collected in N. If, however, hydrogen were 
formed in appreciable quantity it was necessary 
to pump out WN occasionally during a run. The 
apparent volume of the single-pass system was 
738 cc. 


Analysis of Products 


The gaseous products of the reaction were 
analyzed in a small-capacity low temperature 
distillation apparatus of the Podbielniak type. 
The liquid products were fractionated in a 
micro-distillation apparatus, and the fractions 
were analyzed by more or less conventional 
micro-analytical methods. Acetylene was deter- 
mined as previously described.® 


TABLE II. Acetylene production in static runs at 280°C. 








Time, min. 
(corrected Pres- 





Initial for change sure 
ethylene in trans- de- 
pressure, mission crease, Acetylene formed —A4P 

mm of cell) mm cc mm CeH2 

438.5 9 12.5 0.28 0.41 30.0 

440 23 29.5 0.58 0.86 34.4 

201.5 10.5 8.5 0.40 0.59 14.4 

199.5 16.5 14.0 0.48 0.71 19.8 

201 23 17.5 0.69 1.02 17.2 

199.5 23.5 20.0 0.67 0.99 20.2 

201 46 33.0 0.95 1.41 23.4 

201 88 51.5 100 148 34.8 

100.5 10 55 0.39 0.58 9.5 
99.5 38 19.0 0.70 1.04 18.3 

101 10 5.0 0.33 0.48 — (A) 
ES 10.5 3.0 0.43 0.64 4.7 
50.5 18.5 6.5 0.56 0.83 7.8 
50.5 20 7.0 0.56 0.83 8.5 
51.0 31 11.0 0.61 0.91 12.1 
51.0 58 16.5 0.59 0.88 18.8 
23.0 105 25 040 059 4.2 
24.2 21.5 4.4 0.52 0.77 5.7 
11.0 5.0 0.65 0.20 0.30 Be 
10.1 10.5 1.3 0.30 0.45 2.9 
10.0 23 2.0 0.33 0.49 4.1 
94.7 

+4.8 mm of He 10 9.0 0.26 0.39 — (B) 
93.7 
+4.8mmofH. 43 25.5 0.48 0.71 — (B) 








(A) In this run the products were allowed to stand, after turning 
off the lamp, for 15 hours in the reaction system at 280°C before 
withdrawing the sample for analysis. The results show that any loss 
ot acetylene due to thermal reactions is negligible in the length of time 
required for an ordinary run. 

(B) These runs show that the addition of an amount of hydrogen 
much larger than that which might be formed in the reaction cuts 
down acetylene production only by about 30 percent, probably merely 
owing to quenching by hydrogen as well as ethylene. 
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Fic. 5. Static runs—initial rates. 


Materials 


The cadmium used was the purest commercial 
grade. It contained a small amount of zinc. The 
metal was cast into thick rods and turned on a 
lathe into thin coherent turnings. 

Ethylene was obtained in cylinders from the 
Ohio Chemical and Manufacturing Company, 
and was purified by simple bulb to bulb distilla- 
tion. For single-pass runs it was taken direct 
from the cylinder. 

Hydrogen was of commercial grade. It was 
passed over hot copper and through a liquid- 
air trap. 


RESULTS 
(1) Static Runs 


A series of static runs was first made at 280°C 
with various pressures of ethylene. The results 
are given in Fig. 3, and show that the rate of 
pressure decrease increases strongly with in- 
creasing pressure. The lowest curve represents 
runs made with the Corex D filter in place, so 
that only 43261 was transmitted. It is evident 
that 3261 is relatively inefficient compared to 
2288, and that when both lines are present little 
error will be introduced if we disregard the 
effect of 43261, and attribute any effects observed 
to 2288 alone. 

It will be noted that the curves show no 
indications of a pressure rise in the initial stages 
of the reaction, as occurs in the mercury photo- 
sensitized reaction. However, the total rise in 
the case of mercury is small, amounting to from 
1 to 2 millimeters at the pressures employed 
here. In view of the difficulties in thermostating 
at 300°C, such small effects would be difficult to 
establish. However the absence of any signs of 
an induction period in the pressure-time curves 
is proof that if a rise does occur it is considerably 
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Fic. 6. Acetylene production in static runs. Ethylene 
pressures in mm as shown. 


smaller than in the case of mercury. An initial 
rise of 0.5 mm or less at a total pressure of 
200 mm might, however, have escaped detection. 

Figure 4 shows the results of a series of runs 
made with 200 mm of ethylene and varying 
amounts of hydrogen at 280°C. The rate is 
obviously strongly increased by hydrogen, pre- 
sumably largely because, as shown previously, 
the photosensitized reaction with 3261 occurs 
readily in the presence of hydrogen, but not in 
its absence. Since the result is a reaction in 
which both resonance lines participate the 
situation is complex, and no further work was 
done with ethylene-hydrogen mixtures. 

The results given in Fig. 3 show the general 
trend of the effect of pressure on the rate of 
pressure decrease. They are not very accurate 
since a layer of brown polymer very slowly 
forms on the walls of the reaction vessel and 
cuts down the intensity of the radiation entering 
the cell. The effect is, however, quite slow and 
only becomes troublesome when a series of long 
runs are being made. To correct for this a series 
of runs were made with the vessel freshly 
cleaned. The duration of these was as short as 
possible, and only initial rates of pressure 
decrease were determined. These are given in 
Fig. 5. The total time of illumination for the 
whole series was short enough so that no appre- 
ciable polymer deposition occurred. Figure 5 
was then used as a calibration curve for future 
runs to correct for the loss of intensity due to 
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Fic. 7. Ratio of pressure increase to acetylene production in 
static runs. 


polymer formation. In all cases this correction 
was only applied until the intensity had fallen 
off about 30 percent, after which the cell was 
cleaned by taking the apparatus apart, washing 
with fuming nitric acid, and flaming in a stream 
of oxygen. Check runs were made from time to 
time during all phases of the investigation to 
correct for the variation in transmission of the 
cell due to polymer formation. 

The static runs indicated that liquid products 
boiling in the butene-hexene range were formed, 
but on account of the small quantity of material 
available they were not suitable for an investi- 
gation of the liquid products of the reaction. 
However the quantity of gas available was 
ample for acetylene determinations, and in view 
of the importance of acetylene as a product of 
the mercury photosensitized reaction an extended 
series of runs was made to determine the relative 
magnitudes of acetylene formation and of the 
pressure decrease under various conditions. The 
results of these runs are given in Table II. 
The fourth column of the table gives the amount 
of acetylene found, and the fifth column the 
partial pressure this would exert in the system. 
The last column gives the ratio of the pressure 
decrease to the acetylene formed, both quantities 
being expressed in mm. 

The results on the production of acetylene are 
shown in Fig. 6. This shows that, within the 
experimental error, the initial rate of acetylene 
production is independent of the ethylene 
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pressure. As a run proceeds and ethylene is used 
up, the acetylene production falls off and as 
would be expected this effect occurs sooner with 
low pressure runs, where the total amount of 
ethylene in the system is relatively small. The 
initial rate of acetylene production is about 
0.056 mm per min., or 2.7X10-® mole per sec. 
This corresponds to a quantum yield of acetylene 
production of the order of 0.05 to 0.1. 

Figure 7 shows the ratio of the pressure 
decrease to acetylene production, both expressed 
in mm. It will be seen (a) that the ratio increases 
as a run proceeds, and (b) that the initial ratio 
is approximately constant up to 50-mm pressure 
and equal to about 2, and that at higher pressures 
it increases rapidly to nearly 30 at 440-mm 
pressure. Assuming that the pressure decrease 
in the early stages of the reaction is due solely 
to the formation of butene, this means that at 
pressures below 50 mm in the early stages of the 
reaction about 6 ethylene molecules disappear 
by forming butene for every one which disap- 
pears by forming acetylene. At 440-mm pressure 
on the same basis about 58 molecules of ethylene 
disappear by polymerization for each one 
disappearing by forming acetylene. 

It may be noted from Fig. 6 that the acetylene 
concentration in any run eventually reaches a 
constant value. Thus, for example, at 200-mm 
initial ethylene pressure acetylene builds up to 
about 1.5 mm in 60 min. and then remains 
sensibly constant. Calculation shows that at 
60 min. slightly less than half the ethylene has 
disappeared, mostly by polymerization. Since 
the initial rate of acetylene formation is inde- 
pendent of the pressure, it follows that the drop 
in acetylene production is not due merely to the 
decreasing concentration of ethylene. It must 
therefore be due to (a) competing processes 
which destroy acetylene, (b) to a reaction of 
hydrogen, or (c) to a diminished formation of 
acetylene owing to quenching by butene rather 
than ethylene. It seems unlikely that the effect 
is due to a loss of acetylene by direct photo- 
sensitized polymerization, since the quantity of 
acetylene present is always very small compared 
with that of ethylene and butene, and acetylene 
cannot have a much larger quenching cross 
section than the olefins. The indications are that 
hydrogen atoms have little action on acety- 
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lene,'®- although no experiments have been 
made at high temperatures. Acetylene may, 
however, disappear indirectly through a reaction 
of hydrogen atoms, i.e., 


H +CoH, = CoH; 
C.H;+CsH»2= polymer. 


It may also disappear by reaction with higher 
radicals produced from butene. 

It appeared possible that the small amount of 
acetylene formed might arise by a reaction 
involving \3261 rather than \2288. To check 
this runs were made in the static system at 
200-mm pressure with the filter on. These 
showed a production of acetylene of 0.03 cc of 
acetylene per hour. This is negligible as far as 
the present investigation is concerned. 

As mentioned previously, acetylene is sur- 
prisingly inert to the action of excited mercury 
atoms in the presence of ethylene. This suggests 
an excited molecule mechanism for the acetylene 
polymerization, with efficient deactivation by 
ethylene. 

TABLE III. Production of liquid products in circulating 
runs. Temperature = 290°C; trap temperature = — 130°C; 
saturator temperature=300°C; ethylene pressure =ap- 
proximately 20 cm. 














Liquid 
products, 
ce per 
Run —AP Time, min. 
No. mm per hr min. x10 Remarks 
1 0.57 280 8.7 
2 0.43 275 7.0 
3 0.42 245 7.4 
4 0.38 272 5.8 
5 0.33 330 5.3 
6 0.29 305 4.7 
7 0.64 250 10.0 System cleaned out and 
8 0.55 275 — reaction vessel flamed 
9 0.46 275 7.5 in oxygen after run 6. 
10 0.40 280 — 
11 0.39 275 — 
12 0.32 280 — 
13 <0.02 350 0.0 Filter on. Only 3261 


transmitted. 


14 <0.02 280 0.0 Blank run. Lamp not 


on. 








16K. F. Bonhoeffer and P. Harteck, Zeits. f. physik. 
Chemie 139, 64 (1928). 

17H. v. Wartenberg and G. Schultze, Zeits. f. physik. 
Chemie B2, 1 (1929). 

18K. H. Geib and E. W. R. Steacie, Zeits. f. physik. 
Chemie B29, 215 (1935). 

19 J. R. Bates and H. S. Taylor, J. Am. Chem. Soc. 49, 
2438 (1927). 
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To test this point further runs were made in 
the circulating system with the 3261A line alone 
(i.e., with a Corex D filter). This has previously 
been shown to have very little effect on pure 
ethylene. With acetylene alone at 200-mm 
pressure, the acetylene consumption was 90 cc 
per hour in a two-hour run. This will be a 
minimum rate since the transmission of the cell 
was greatly diminished by the deposition of 
polymer on the walls, and the initial rate must 
have been much higher. After cleaning the 
reaction vessel an experiment was made with 
200-mm total pressure of a mixture consisting of 
77 percent ethylene and 23 percent acetylene. 
In one hour 1.9 cc of acetylene was consumed, 
together with a trace of ethylene. It follows that 
in a 3 to 1 mixture of ethylene and acetylene 
the rate of polymerization of acetylene has 
fallen to less than 1/47 of the rate for pure 
acetylene. It seems highly unlikely that the 
quenching cross section of ethylene for *P, 
cadmium atoms can be 10 times greater than 
that of acetylene, and it must therefore be 
concluded that excited acetylene molecules are 
easily deactivated by collision with ethylene. 


(2) Circulating and Single-Pass Runs 


The runs in the circulating and single-pass 
systems were mainly of long duration, and 
were designed to accumulate products for 
analysis, especially liquid products. During such 
long runs a very thin film of highly polymerized 
material built up on the walls of the vessel and 
cut down the light transmission. Examples of 
the over-all rates of successive runs of about 4 
hours duration are given in Table III. These 
runs were made primarily for the purpose of 
collecting liquid products, and the trap was 
maintained at —130°C to diminish secondary 
reactions of the products of higher molecular 
weight. 

The results on the production of acetylene in 
circulating and single-pass runs are given in 
Fig. 8. The times given in the figure have been 
corrected for the drop in transmission of the 
reaction cell in successive runs, using the data 
of Table III. In these runs the temperature of 
the reaction vessel was 290°C, and that of the 
saturator was 300°C. The pressure was 20 cm 
in both cases. In the circulating runs the trap 
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Fic. 8. Acetylene production in all types of run; open 
circles =circulating runs; filled circles=single-pass runs; 
half-filled circles = static runs. 


temperature was —130°C. In the single-pass 
runs the flow rate was 20 cc per min. at N.T.P., 
corresponding to a contact time of about 2.0 min. 
Some single-pass runs were also made at higher 
flow rates, viz., 136 and 394 cc per min. In these 
runs, however, considerable difficulty was en- 
countered in obtaining satisfactory saturation 
with cadmium vapor. The results were therefore 
rather unreliable and are not reported, but as 
far as they went they indicated that the change 
in flow rate made no appreciable change in the 
rate of acetylene production. In the single-pass 
runs all the gas flowing was condensed in a 
liquid-air trap and an aliquot was analyzed for 
acetylene. 

Some explanation is required of the fact that 
the maximum acetylene production in static 
runs of long duration is approximately 1 cc, 
while in circulating runs it rises to above 4 cc. 
(Obviously it will increase linearly with time 
indefinitely in single-pass runs.) There are 
several factors which together will account for 
the larger acetylene production in circulating 
runs of long duration. 

(a) The volume of the circulating system is 
larger by a factor of nearly 2. Hence if the 
maximum occurred at the same partial pressure 
it would be nearly 2 cc in the circulating system. 

(b) The vapor pressure of acetylene is about 
10 mm at —130°C, which corresponds to 12 cc 
of acetylene in the gas phase. However, it is 
probable that some solubility of acetylene in 
the butene, etc., in the trap will occur, and thus 
acetylene may be removed from the system, and 
its apparent steady-state concentration increased. 

(c) The destruction of acetylene may involve 
reaction with butene or higher hydrocarbons, as 
pointed out previously. Since these are trapped 
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out in the circulating system, acetylene can 
accumulate to a greater extent. 

(d) In the circulating system the total pressure 
was kept constant throughout a run at 20 cm 
by changing the level of the mercury reservoirs 
of the circulating pump. Since higher products 
were trapped out this 20 cm consisted only of 
hydrogen, ethylene, acetylene, and a trace of 
butene. Hence the partial pressure of ethylene 
was maintained much higher than in the later 
stages of a static run. It follows that acetylene 
should build up to a greater extent if the relative 
amount of quenching by ethylene and acetylene 
is the deciding factor. 

However, in spite of the differences discussed 
above, and of the lower accuracy of circulating 
and especially of single-pass runs, it appears 
from Fig. 8 that the initial rate of acetylene 
production is approximately the same in static, 
circulating, and single-pass runs. It may be 
concluded, therefore, that the initial rate ob- 
served in static runs is not influenced by second- 
ary reactions involving butene, higher olefins, 
or hydrogen. 


(3) The Products of the Reaction 
(a) Acetylene 


The production of acetylene has already been 
discussed in detail. 


(b) Hydrogen 


The detection of an amount of hydrogen of the 
same order of magnitude as that of acetylene is 
very difficult, and hence it was determined 
merely as ‘‘non-condensable” gas. Thus in a 
typical 33-hour circulating run at 290°C and 
20-cm pressure the gas non-condensable at 
—180°C was determined. This, including the 
vapor pressure of ethylene, was equivalent to 
2.0 cc of gas at N.T.P. In a corresponding 
blank run the non-condensable gas, including 
the vapor pressure of ethylene, was 0.4 cc. 
Hence the non-condensable products amounted 
to 1.6 cc or 0.46 cc per hour. The acetylene 
formed was 1.0 cc per hour under the same 
conditions. The hydrogen formed therefore 
appears to be equivalent to approximately half 
the acetylene. It would appear from this result, 
which was typical of a number of determinations, 
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that some of the hydrogen formed with acetylene 
is used up. However, in view of the very small 
amount of hydrogen the possibility of its loss by 
adsorption, solubility, etc., cannot be excluded. 


(c) Propylene 


In low temperature distillations of the products 
some indications of the presence of propylene 
were obtained. Thus in 2 nine-hour single-pass 
runs at 20-cm pressure distillation’ analyses 
indicated 74 and 91 cc of propylene at N.T.P. 
A nine-hour blank under the same conditions, 
however, gave 86 cc, and hence the propylene 
appears to be merely an impurity in the ethylene. 
This is not a large amount of impurity, since in 
such a single-pass run over 11 liters of ethylene 
is used. As might be expected, little propylene 
was found in circulating runs where the amount 
of ethylene used was much smaller. Here blank 
runs and ordinary runs both indicated the 
presence of about 0.5 percent propylene as an 
impurity in the ethylene. It may therefore be 
concluded that propylene is not produced in 
the reaction in appreciable amounts. 


(d) Butene and Higher Olefins 


A typical distillation analysis of the accumu- 
lated products from 7 circulating runs each of 
approximately 4 hours duration, with the trap 
at —125°C and an initial pressure of 200 mm 
gave 


Sample1 CsHs b.p. —53°to —44°C 19 ccat N.T.P.=0.03 g 
Sample2 CsHs_ b.p. — 8° to + 4°C 266 cc =0.61 g 
Sample 3 ‘“‘CeHi2”" b.p. 17° to 70°C 1.01 g 
Sample 4 b.p. >70°C 0.40 g 


These results suggest that the products other 
than hydrogen and acetylene are about } butene 
and the remainder higher olefins. Analytical 
data are given later to substantiate this. 

A number of circulating runs were made 
specifically to determine the rate of production 
of butene. In these a distillation analysis of the 
products of each run was made in the low 
temperature still. Only butene could thus be 
determined: the remaining liquid products were 
merely weighed. The results are given in Table IV. 

The data of Table IV are in general agreement 
with those reported above, i.c., from 25 to 30 
percent of the olefinic products appear to be 
butenes. From these results we can make an 





TABLE IV. Distillation analyses of products of circulating 
runs. Temperature = 290°C; trap temperature = — 125°C; 
initial pressure-200 mm. 








Time, CsH«4, CsHs, CsHs, Liquids, 
min. cc cc 4 g 


250 252 26.5 0.063 0.16 

235 — 20 #£0.047 0.14 

220 195 31 0.072 0.19 System cleaned after 
245 248 35 0.082 0.21 second run, hence 
Blank 294 2.5 — _ 0.00 higher transmission in 
Blank 278 00 —_ 0.00 following runs. 


Remarks 











approximate check of the relative rates of 
acetylene production and of pressure decrease. 
From Table IV for a four-hour run in a clean 
reaction vessel the rate of butene production is 
about 9 cc per hour. The liquid products are 
about 2.7 times this by weight. Assuming that 
they are hexene, this corresponds to 16 cc. 
Now the pressure decrease, expressed as cc of 
gas disappearing is given by 


—AP= 1 XC4Hs+2 XCeHie 
=41 cc per hour. 


From Fig. 8 the production of acetylene is 
approximately 1.0 cc per hour. Hence for the 
circulating runs 


—AP/C:H2=41. 


Actually, from Fig. 7 the ratio in the static runs 
is 35 at 90 min. and is rising with time. The 
agreement between the static and circulating 
runs is thus excellent. 

As a check on the composition of the fraction 
of the products boiling at about 0°C, which has 
been reported as butene above, a molecular 
weight determination was made by combusting 
a known volume of gas and determining the 
weights of carbon dioxide and water formed.” 
Two determinations gave 53.58 and 53.79, as 
compared with the theoretical value for butene 
of 54.1. For the composition of the gas the 
combustion analyses gave 


C= 84.44, 84.96 percent 
H = 14.78, 14.56 
C/H= 5.75, 5.87. 


(Theoretical C/H for CsHs=6, for CsHe=8, for 
C4Hy)=4.8.) 


20... Marion and A. E. Ledingham, Ind. Eng. Chem. 
Anal. Ed. 13, 260 (1941). 
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As a further check on the composition of the 
0° fraction, a ‘‘butene” cut of —8° to +4° was 
further separated by distillation into two cuts, 
(a) 70 percent of the total, boiling in the range 
—8° to 0°, and (b) 30 percent boiling in the 
range 0° to +4°. Samples (a) and (b) were then 
quantitatively brominated. 

The results were: 

Percent of double bonds 


assuming a molecular 
weight of 56 


Sample (a) 92.4 
(b) 90.5 
It thus appears probable that a small amount of 
butane is present with the butene. Pentene or 
hexene with the butene would give a similar 
result, but this seems unlikely in view of the 
agreement between the two fractions. 


(e) Higher Liquids 


As pointed out above the higher liquids 
amounted to about twice the butene. Most of 
these were in the 17° to 70° cut. The small 
amount of material made it difficult to fractionate 
this cut further, but the indications were that 
it was largely hexene, with some pentene. 
There was also a small amount of liquid product 
which boiled in the range 70° to 120°C, and a 
still smaller amount which boiled above this. 

The liquid products from a number of circu- 
lating runs were roughly fractionated into cuts 
boiling in the ranges 30° to 40°C, 60° to 75°C, 
100° to 120°C, and over 120°C. These were 
analyzed by combustion with the following 
results: 


30° to 40° cut C/H=5.97, 6.07 
60° to 75° cut C/H=5.97, 5.92, 6.14, 5.97 
100° to 120° cut C/H=6.03, 6.14 
over 120° cut C/H=6.10, 6.09. 


Since all olefins give C/H=6, these results 
strongly confirm the previous conclusion that 
the liquid products are olefinic. 

As a further check on the olefinic nature of 
the products, a number of brominations were 
made followed by fractional distillation of the 
bromides and micro-analyses of the fractions for 
bromine. The products used were a composite 
sample from a number of circulating runs. 
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(a) The analysis of the bromides from the 0° 
fraction was 
Found Br=76.15 percent 
calc. for CsHsBr2= 74.06 
calc. for CsH¢Brs= 85.57. 


Hence if any diolefins are present they are not 
in large amount. 

(b) The 60° to 75° fraction was brominated in 
about 25 cc of carbon tetrachloride and the 
product, after removal of excess bromine and 
solvent under reduced pressure at 50°C, was 
fractionated. Fraction I was a colorless oil with 
b.p. 75°/2 mm. This was the bulk of the product. 
On redistillation it again gave b.p. 75°/2 mm. 
Fraction II was of relatively small amount, with 
b.p. 80° to 120°/2 mm. The bromides from 
fraction I were analyzed with the results shown 
in Table V. Fraction I was then redistilled twice 
and a fraction of about } the amount, coming 
over at about 68° to 70°/2 mm, was again 
analyzed. This gave C=28.48, 28.48, H=4.75, 
4.79, Br=68.07. The results suggest that the 
bulk of the product is a single Cs, compound, 
i.e., hexene, but that some diolefins or pentene 
may be present. The latter is probably the more 
likely alternative. 


DISCUSSION 
The main facts established above are: 


(a) The products of the reaction are olefins, 
together with smaller amounts of hydrogen and 
acetylene. 

(b) Since circulating runs with the trap at 
— 80° and —125°C and single-pass runs all give 
much hexene, it appears that hexene is not 
merely formed by secondary. reactions of butene, 
but arises in reactions which are an integral part 
of the mechanism. 

(c) The concentration of acetylene rises to a 
maximum and then remains sensibly constant. 
The initial rate of acetylene formation is approxi- 
mately independent of pressure. 

(d) The ratio of pressure decrease to acetylene 
formation is high at high pressures and falls to 
about 2 at lower pressures, i.e., the rate of 
pressure decrease rises rapidly with increasing 
pressure. 

(e) There is no measurable pressure rise at the 
beginning of the reaction as occurs with mercury. 


There appear to be 4 possible mechanisms for 
the primary step: 


CoHi+Cd(*P,) =Cd(4So) +2CH2 (2) 
=Cd(1So)+CsH2+He (3) 
= Cd(So) +CoH4* (4) 
gel 5) 

or =CdH+C;2H; 


Reaction (2) appears to be unlikely on 
energetic grounds. In any case the predominance 
of C4 and C, olefins in the products and the 
absence of propylene and methane seem to rule 
it out. 

Reaction (3) can certainly be ruled out as the 
primary step since acetylene is not the major 
product of the reaction, and evidence has been 
cited above to show that acetylene is surprisingly 
inert in the presence of ethylene. 

Reaction (4) is favored by analogy with the 
mercury photosensitized reaction, and reaction 
(5) is favored by analogy with other cadmium 
photosensitized reactions. 

As we have seen, in the case of the mercury 
photosensitized reaction the main steps are 


C.H,4* =C:H2+Hae, (6) 
C.H,+Co2H4* =2C2H, (7) 


and the subsequent polymerization arises by a 
free radical mechanism based on H atoms formed 
by 

Ho+Hg(*P1) = Hg("So) +2H. (8) 


In view of the large amount of polymerization 
compared with acetylene formation in the 
present investigation, it does not seem possible 
to explain all the polymerization on this basis. 
However, it is possible that the higher excitation 
energy in the case of cadmium causes (7) to be 
largely replaced by 


C.H4+C2H4* = CaHs. (7a) 


TABLE V. 








Calculated for 








Found CeHi2Bre CeHioBra CsHioBre 
C 28.48, 28.75 29.52 22.55 26.10 
H 4.77, 4.91 4.92 2.35 4.35 
Br 67.71, 67.64 65.55 75.11 69.55 
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It seems more likely, however, that by virtue of 
the higher excitation energy reaction (5) is also 
an important primary step. If so the subsequent 
polymerization can proceed simply by a free 
radical mechanism involving H and C;Hs, e.g., 


H+C.H,=CoH; 
CsH4+CoHs=CaHy 
CoH4+Cy4Hy = CoH, etc. 


As a matter of fact a mechanism could be 
devised based on reaction (5) as the only 
primary step, provided that we assume acetylene 


KEITH J. 
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to be formed by disproportionation, 
C2H3+C2Hs=CoH2+CoHs. 


However, all this is entirely speculative. The 
main features of the work are that there is 
practically no polymerization of ethylene by 
Cd(5*P1), while polymerization proceeds effici- 
ently with the higher energy associated with 
Hg(6*P1) and Cd(5'P,). 

We are indebted to Dr. L. Marion and 
Mr. A. E. Ledingham of this laboratory for 
making the micro-analyses. 
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The potential-energy surfaces relevant to a number of simple processes involving excited 
(?P) sodium have been constructed, and the mechanism of the reactions discussed in terms of 
them and of the theory of absolute reaction rates. It is shown that quenching by atoms is 
normally a very inefficient process, the effect of additional atoms in a quencher being to stabilize 
a quenched complex, commonly a polar compound, which finally decomposes to give the 
deactivated atom and a vibrationally excited product. The activated state may be either at the 
top of a rotational-energy barrier or at a crossing point to a polar surface, and the theory is 
shown to be consistent with the experimental values of the quenching cross sections. In the 
cases examined the quenching is of a physical nature, the energy passing chiefly into vibrational 
energy, with a certain amount going into translational and rotational degrees of freedom. 
Saturated hydrocarbons appear to quench like hydrogen except that there is probably a small 
energy of activation, but the unsaturated ones quench by means of an interaction between the 
sodium atom and the unsaturated group, and the efficiency is considerably higher. 


INTRODUCTION 


FTER an atom in the gas phase has ab- 
sorbed energy and become excited elec- 
tronically, three distinct processes may ensue: 

(1) The electronic energy may be re-emitted in 
the form of radiation (fluorescence) ; 

(2) The atom may collide with an atom or 
molecule and become deactivated, the excitation 
energy passing into translational, rotational or 
vibrational degrees of freedom, but not inducing 
reaction (quenching) ; or 


* Contribution No. 1031 from the National Research 
Laboratories, Ottawa, Canada. 





(3) The atom may collide with another atom 
or molecule and initiate chemical reaction in the 
latter (photosensitization). All of these processes 
may, and frequently do, occur simultaneously in 
a given reaction system, and the individual rate 
constants may be separated by investigating the 
effects of changing the concentrations of the 
reacting species on the intensity of the emitted 
fluorescent radiation or on the rate of formation 
of the reaction products. Thus if A denotes the 
atom which is excited electronically to A* by the 
absorption of radiation, and M is a foreign gas 
molecule, possible reactions, with their corre- 
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sponding rate coefficients, are 


A+hy=A* (kx), (1) 
A*+A=A’+A’ (ka), (2) 
A*+M=A’+M’ (ks), (3) 

A*=A+hy (ky). (4) 


A’ denotes an A atom which has an abnormal 
amount of translational energy, while M’ may be 
excited translationally, vibrationally or rota- 
tionally, and may subsequently undergo reaction. 
The assumption of a steady state for the excited 
atom leads at once to the equatien! 








I ky 
= (5) 
Io RatRocathkseu 
1 
= (6) 
1+hkorcatksreu 


where 7(=1/k,) is the mean life of A*; J is the 
intensity of radiation emitted in the absence of 
foreign gas and extrapolated to zero concen- 
tration of A, and J is that in the presence of A and 
M at concentrations ca and cy. Since 7 can 
usually be determined readily, ke and k3 can be 
found from the effect on the ratio I/Jo of varying 
the concentrations of A and M. 

The rates of such processes as (2) and (3) have 
previously been discussed almost entirely in 
terms of classical collision theory according to 
which the number of collisions per second per 
cubic centimeter, for unit concentration (in 
molecules per cc) of the two colliding molecules, 


is given by 
mx+mp\ 7} 
z=o} 8kT(“~—) , (7) 
MamB 


where o is the distance between centers on 
collision, and ma and mg are the actual masses of 
the colliding species. This is related to the 
velocity coefficient k by the equation 


k=Ze-=!RT cc molecule sec.—', (8) 


and if the activation energy E can be assumed to 
be zero, as is probably the case in most of these 


10. Stern and M. Volmer, Physik. Zeits. 20, 183 (1919). 
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reactions, k and a are related by 


ma+mp\ 7 

bo brir(—“—"*)| : (9) 

MamMmp 
The experimental results for quenching and 
primary photosensitization processes are fre- 
quently expressed in terms of o”, known as the 
“quenching cross section,”’ and derived from k by 
means of this equation. It will later be shown 
that this equation is also obtained when the 
theory of absolute reaction rates is applied to the 
processes, and that values of the quenching cross 
section in satisfactory agreement with experi- 
ment can be derived when the potential-energy 

surfaces have been constructed. 

Most of the experimental work on quenching 
and photosensitization has been done with ex- 
cited sodium, mercury and cadmium. The 
simplest reactions to interpret theoretically are 
those of sodium, since the potential-energy sur- 
faces are fewest in this case, and because the 
quenching of a ?P sodium atom must involve the 
formation of sodium in its lowest, *S, state since 
this is the only lower state; with mercury and 
cadmium, on the other hand, quenching often 
occurs by the transfer of the atoms from their 
3P; levels to the metastable *P» levels (see Part 
II, following paper). A further simplification in 
the case of sodium arises from the fact that the 
amount of energy, 48.3 kcal. per mole, liberated 
when an excited (?P) sodium atom becomes 
deactivated to the normal (2S) state is not 
sufficient to dissociate the average quenching 
molecule. 

The quenching of excited (?P) sodium by a 
number of gas molecules has: recently been 
studied by Norrish and Smith.? Some of the data 
are given in Table I, which also includes, for 
comparison, the data for mercury and cadmium. 
The values given for sodium are those of Norrish 
and Smith, and of Terenin and Prileshajewa (the 
halogens). 

Of the previous theoretical treatments of re- 
actions of this type, most have been of a general 

2 R. G. W. Norrish and W. MacF. Smith, Proc. Roy. Soc. 
A176, 295 (1940); see also R. Mannkopff, Zeits. f. Physik 
36, 315 (1926); J. G. Winans, ibid. 60, 631 (1930); A. 
Terenin and N. Prileshajewa, Zeits. f. physik. Chemie 
B13, 72 (1931), Physik. Zeits. Sowjetunion 2, 337 (1932); 


and B. Kisilibasch, V. N. Kondratjew, and A. Leipunsky, 
ibid. 2, 201 (1932). 
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Fic. 1. Po- 
tential-energy 
curves for the 
quenching of an 
excited atom A* 
by an atom M 
(Case I; cf. Eq. 
(10)). 











A-M distance 


nature. Teller* has discussed the general principles 
involved in the dissipation of electronic into 
other types of energy, while Kallmann and 
London,‘ and Morse and Stuechelberg,® have 
made calculations based upon assumed potential- 
energy surfaces. Recently, a number of papers 
have treated special aspects of the application of 
the theory of absolute reaction rates to reactions 
involving excited electronic states,* and in par- 
ticular Magee and Ri’ have discussed the quench- 
ing of sodium by hydrogen. In the remainder of 
the present paper the general forms of the 
potential-energy surfaces upon which the simpler 
reactions involving excited sodium atoms take 
place will be deduced, and the mechanisms and 
rates of the reactions discussed in the light of 
them. 


QUENCHING BY ATOMS 


The set of potential-energy curves for the two- 
atom system formed from the interaction of the 
excited and the quenching atom belongs to one of 
at least three classes. In the first place the 
interactions between the atoms may be of such a 
nature that the lowest curve for the interaction of 
A* with M is intersected by one of the curves for 
A-—M;; this case is illustrated in Fig. 1 in which 
the rounding-off of the curves due to resonance 
between them is shown. If the probability of 
crossing from the upper to the lower surface is 
equal to p, the probability that atoms A* and M 


3E. Teller, J. Phys. Chem. 41, 109 (1937); Ann. N. Y. 
Acad. Sci. 41, 173 (1941); see also J. v. Neumann and E. P. 
Wigner, Physik. Zeits. 30, 467 (1929). 

4H. Kallmann and F. London, Zeits. f. physik. Chemie 
B2, 207 (1929). 

5P. M. Morse and E. C. G. Stuechelberg, Ann. d. 
Physik 9, 579 (1931). 

6 For a summary, with references, see S. Glasstone, K. J. 
Laidler, and H. Eyring The Theory of Rate Processes 
(McGraw-Hill Book Co., 1941), Chapter VI. 

7J. L. Magee and T. Ri, J. Chem. Phys. 9, 638 (1941). 
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approaching one another along the upper curve 
will separate along the lower can readily be 
shown to be® 


k=2p(1—p). (10) 


Since 1>p>0, x has a maximum value of 4 
when p=3, while, when p~1 or p~0, x is very 
small. If there is no difference of multiplicity 
between the two states the interaction energy 
will usually be large owing to a considerable 
overlapping of orbitals; p will therefore be ~0 
and the rate of quenching low. The reaction is 
also slow for the same reason in the rather 
common case of a change of multiplicity, when 
p~i. Only in especially favorable cases, there- 
fore, will the quenching efficiency be appreciable. 

In the second case the upper and lower curves 
do not intersect, but both are intersected by a 
third curve, so that quenching can occur by 
means of two transitions; this case is illustrated 
in Fig. 2. The curve HBFK is usually a polar 
curve, for the ionic molecule AtM-. In order to 
derive an explicit expression for the transmission 
coefficient of such a reaction, an expression is first 
obtained for the fraction of systems which, 
beginning at A, find themselves traveling towards 
F on the polar curve BF, and then for that 
fraction of the latter which finally arrive at E. 
The first fraction is found to be 2(1—p)/(2—p)), 
where p is the probability of crossing at B; the 
second fraction is 


20(1—<c) 
1—p*{o*+(1—0)?}) 





where o is the probability of crossing at F. The 


TABLE I. Quenching cross sections (A*). 











Quenching gas Na(?P) Hg(*P;) Cd(3P1) 
Inert gases ~0 ~0 ~0 
He 7.4 6.01 0.67 
Ne 14.5 0.192 0.021 
CO 28.0 4.07 0.14 
Ie 40.0 — oo 
Bre 100.0 “——- — 

I 40.0 — — 
CH, 0.11 0.059 0.012 
C:He 0.17 0.421 — 
CH, 44.0 48.0 — 
C3He 52.0 oa ae 








8 W. A. Noyes and F. C. Henriques, J. Chem. Phys. 7, 
767 (1939). 
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net transmission coefficient is therefore 


2(1—p) 20(1—<¢) 
c= . . (11) 
2—-p 1—p*{o?+(1—c)?} 





Other things being equal this function is 
greatest when p=0; its maximum value is } when 
p=0and o=3, while it is zero when ¢=0 or c= 1. 
When there is no change of multiplicity the 
resonance energies are usually large so that p~0 
and o~0; «x is therefore small in this case. 
Otherwise, if there is a change of multiplicity, 
o~1 so that x~0. 

The third possible case is when the polar curve 
does not lie sufficiently low to cut both curves; 
quenching is extremely inefficient in this case. In 
general, therefore, atoms are poor quenchers, 
owing to the high probability of the dissociation 
of the complex formed by the two atoms to give 
the original excited atom; the efficiency with 
which electronic energy passes into translational 
is thus low. That atoms are poor quenchers seems 
to be confirmed by experiment; thus the inert 
gases quench very weakly.® The case of the 
quenching of excited sodium and thallium by 
iodine atoms appears, however, to be in contra- 
diction to this conclusion, since Terenin: and 
Prileshajewa’® give values of 40 and 30A2, re- 
spectively, for o?. It is conceivable that this 
arises from the case of s-~} and p=0, but in view 
of the indirect method by which the results were 
obtained it would seem to be more likely that the 
results are in error, the quenching actually being 
effected by iodine molecules. 


QUENCHING BY HYDROGEN 
Quenching by Hydrogen Atoms 


Before considering the quenching of excited 
sodium by hydrogen molecules, it is necessary to 
construct the potential-energy curves for the 
system Na—H,;; these are shown in Fig. 3.4% A ?P 
sodium atom interacts with a 2S hydrogen atom 
to give 12+, *Z+, 4II and “II states, of which the 
first is the lowest and attractive. The curve for 
this state can be plotted making use of the 


®R. Mannkopff, Zeits. f. Physik 36, 315 (1926); L. V. 
Hamos, ibid. 74, 379 (1932). 

10 A, Terenin and N. Prileshajewa, Zeits. f. physik. Chemie 
B13, 72 (1931); Physik. Zeits. Sowjetunion 2, 337 (1932). 
1 Cf. R. S. Mulliken, Rev. Mod. Phys. 4, 7 (1932). 





spectroscopic data;” the value of Do, the dis- 
sociation energy plus the zero-point energy, is 
34.15 kcal. The form of the curve is remarkable in 
that in spite of being very shallow at its minimum 
it rises sharply in the region of 4-6A;;" this arises 
from resonance with the ionic state, its effect 
being to give negative values for the band- 
spectrum constants a, and X.Wwe. 

A normal (2S) sodium atom interacts with a 
hydrogen atom to give an attractive '=*+ state and 


H Atm" 





Fic. 2. Potential-en- 
ergy curves for quench- E 
ing by an atom M, the 
process involving the 
intermediate formation 
of a polar complex 
(case IT; cf. Eq. (11)). 








A-M distance 


a repulsive *2+ state, of which the former can be 
constructed from the spectroscopic data. The 
only other curve for the system which falls 
sufficiently low to be of any interest is that for the 
ionic Na+—H state. For an infinite Na+—H- 
separation this is higher than the curve for 
Na(?S)—H(2S) by an amount J(Na)—E(H) 
where J(Na) is the ionization potential of sodium 
and E(H) is the electron affinity of the hydrogen 
atom ; this quantity has the value 101.5 kcal. The 
form of the curve for the ionic state is given 
approximately by the expression 


=—é/r+br-", (12) 


where V is the energy relative to the separated 
ions, ¢€ is the electronic charge, r the distance 
between the Na+ and H- ions, and 6 and 7 are 
constants ; the first term represents the Coulombic 
attraction of the ions, while the second is the 
repulsive contribution due to exchange forces. 
The values of 6 and can be evaluated from the 
crystal data collected by Sherman ;" at values of r 
greater than 4A the influence of the second term 
becomes negligible, and the energy in kcal. per 


% Values from H. Sponer, Molekiilspektren (J. Springer, 
Berlin, 1935), Vol. I. 

13 Cf, R. S. Mulliken, Phys. Rev. 50, 1028 (1936), who 
treats LiH similarly. 
14 J. Sherman, Chem. Rev. 11, 93 (1932). 
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Fic. 3. Potential-energy curves for the system Na—H. 


mole may be written as V= —330/r, where r is 
expressed in angstroms. 

Examination of the curves in Fig. 3 shows that 
for quenching to take place the system must 
proceed along the Na(?P)+H(2S) curve and 
either make a direct transition to the repulsive 
3>+ curve, with the ultimate formation of 
Na(??S)+H(?2S) which will move away from one 
another with an extra relative kinetic energy 
equal to the excitation energy of sodium, or make 
two transitions, first to the polar curve and then 
to the repulsive curve. The first possibility is 
case I of the previous section (Fig. 1 and Eq. 
(10)) ; as there is a change of multiplicity p~1, so 
that x~0 and the rate is low. The second possi- 
bility is case II (Fig. 2 and Eq. (11)) with p~0 
and ¢~1, so that again «~0. The normal course 
of events is therefore that the colliding atoms 
readily form a polar compound Na*H-, but that 
this, being unable to get rid of its excess of 
vibrational energy, redissociates at once. 


Quenching by Hydrogen Molecules 


The potential-energy surfaces for the three- 
body system Na—H—H may be represented in a 
three-dimensional diagram in which the coordi- 
nates are potential energy, Na—H distance and 
H —H distance. The two profiles of the surfaces 
which correspond to infinite separation of Na and 
H and of H and H are shown on two sides of the 
cube in Fig. 4. The process of interest, the 
quenching of Na(?P) to Na(?S) is represented by 
a transition from B to A in this figure. The 
manner in which the curves are connected by 
surfaces may now be discussed. 

When a hydrogen molecule in a !2+, state 
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interacts with a 2S sodium atom to give a linear 
complex, the resulting state is a 22+ state, which 
is therefore the lowest for the system, the surface 
describing it connecting the two lowest curves. 
This surface corresponds to that upon which the 
reaction H+H, takes place and is of the same 
general form; its energy may be expressed by the 
London equation" taking the negative sign: 


E=Q-[3{(a—8)?+(8—vy)?+(y—a@)?} J}, 


where Q is the sum of the three Coulombic 
energies and a, 6 and vy are the exchange energies. 
The form of the surface has been calculated using 
the semi-empirical method, the potential-energy 
curves being estimated from the data given by 
Sponer.” The proportion of Coulombic energy 
was taken as 12 percent for the H—H bond and 
22 percent for the Na—H bond; this latter figure 
is the average of 12 percent for H—H and 32 
percent for Na—Na, these being the most 
satisfactory figures for bonds involving electrons 
having principal quantum numbers of 1 and 3, 
respectively.!® The interactions of the higher 
states are neglected. 

The interaction of a 2S sodium atom with a 
3>+,, hydrogen molecule gives rise to two states, 
4y+ and 22+; similarly H(?S)+NaH(®=*) gives 
4>+ and *2+. The two repulsive curves in Fig. 4 
for Hz and NaH are therefore connected by two 


(13) 
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Fic. 4. Potential-energy surfaces for the system 
Na—H-—H (linear configuration). The curves on the left- 
and right-hand faces are for Na—H and H—H, respec- 
tively, and the manner in which they are connected by 
— the cube is described in the text. AB= MN 
= 48.3 kcal. 


18 F, London, Probleme der modernen Physik (Sommerfeld 
Festschrift) (1928), p. 104; Zeits. f. Elektrochemie 35, 552 
(1929). 

1% J. H. Bartlett and W. H. Furry, Phys. Rev. 38, 1615 
(1931); N. Rosen and S. Ikehara, ibid. 43, 5 (1933). 
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surfaces. The same conclusions have been arrived 
at by the method of bond eigenfunctions;!’ 
neglecting double-exchange integrals, and making 
other approximations, the energy of the *2* state 
is found to be 


E=Q+([3{(a—8)?+(B—y)?+(y—a@)*} J, (14) 
while that of the ‘2+ state is 
E=Q-a-8-y. (15) 


As the exchange integrals are of negative sign 
these surfaces represent highly repulsive states; 
this being so, and since the energy of Na(2S)+2H 
(M in Fig. 4) is considerably higher than that of 
Na(?P)+H2(!2+,) (B in Fig. 4) it is clear that a 
quenching process utilizing one of these states 
would have a very high activation energy. As 
such energies are not involved the surfaces are 
not used, but it will be seen later that the 
repulsive surfaces are of significance in the 
reactions of excited mercury and cadmium. 

Of the surfaces involving excited sodium, only 
the lowest, connecting the curves for Na(?P) 
+H,(!2+,) and NaH('Z+)+H(?S), and of desig- 
nation *=+, is of any interest. As Magee has 
shown,'® the triangular configuration of atoms is 
probably the most stable in this case, but the 
error in assuming linearity is not large. A plot of 
the surface for the linear configuration has been 
made in the same manner as for the lowest 
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Fic. 5. Potential-energy profiles in the plane ABCDE 
of Fig. 4, corresponding to the normal H—H separation. 
AB=48.3 kcal., BP=c. 200 kcal. 


17J. C. Slater, Phys. Rev. 38, 1109 (1931). 
18 J. L. Magee, J. Chem. Phys. 8, 677 (1940); J. L. 
Magee and T. Ri, ibid. 9, 638 (1941). 
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Fic. 6. Potential-energy profiles in the plane FGHIJKL 
“ Fig. 4, corresponding to an extension of the H—H 
ond. 


surface, of sufficient accuracy to establish two 
facts. The first is that the lowest point in the 
valley representing the approach of Na(?P) to He 
runs parallel to the Na — H axis up to small (c. 1A) 
values of the Na—H distance; on this surface the 
excited sodium atom therefore approaches close 
to the hydrogen molecule without appreciable 
extension of the bond in the latter. The second 
point is that the surface is at all points well above 
the lowest surface, a conclusion that can also be 
accepted for the triangular configuration. The 
possibility was originally envisaged that the 
interaction between excited sodium and hydrogen 
might be so different from that between normal 
sodium and hydrogen that the two surfaces might 
cut, with the result that quenching would involve 
a simple transition from the upper to the lower 
surface.!* Such a possibility is, however, excluded. 

The only remaining possibility is that the 
surface representing the ionic state Nat —-H-—H 
cuts both upper and lower surfaces, so that the 
reaction proceeds via a transition first from the 
upper to the polar surface, and then from the 
polar to the lower surface. The curves repre- 
senting NatH~-(!2+) and H~-—H(?*) are shown 
in Fig. 4. The construction of the former has 
already been discussed. The molecule H2- 
has been treated theoretically by Eyring, 
Hirschfelder, and Taylor®® who find a dissociation 
energy of 37.4 kcal. and an equilibrium separa- 
tion of 1.8A. At the separation of 0.74A, equal to 

12 Cf. M. G. Evans, Trans. Faraday Soc. 35, 59 (1939). 


20H. Eyring, J. O. Hirschfelder, and H. S. Taylor, J. 
Chem. Phys. 4, 484 (1936). 
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Fic. 7. Potential energy surfaces for the linear system 
Na—X-—X, where X is a halogen atom. 


the separation in He, the energy is considerably 
(by more than 50 kcal.) higher than corresponds 
to dissociation into H+H~-. The sodium ion 
Na*('S) interacts with the H- ion (22+) to give a 
single 22+ state in which it is assumed for the 
purpose of calculation that the charge is located 
at the center of the hydrogen molecule: The 
surface is best constructed by assuming that the 
interaction between Nat and H:~ obeys the 
same law of force as that between Nat and H-, 
which is justifiable as the exchange forces in the 
region of interest are negligible in comparison 
with the coulombic ones. Owing to the high 
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repulsion in H2~ at the separation of 0.74A the 
polar surface does not cut the upper surface in the 
profile that corresponds to the normal H—H 
distance; this profile is shown schematically in 
Fig. 5, in which the curve for the upper surface is 
drawn with a slight minimum which Magee and 
Ri have shown to exist for the triangular con- 
figuration. For crossing to take place there must 
be an extension of the H—H bond until the polar 
surface cuts the upper one. This extension is 
clearly not as great as to 1.8A, the normal 
distance in H2"; if it were, the effect would be to 
bring the polar surface down as low as possible, 
but the activation energy would be high (Fig. 6). 
The value of r, at the crossing point in this profile 
is about 8A. The crossing point through which 
reaction actually proceeds will clearly correspond 
to a smaller extension of the bond, so that the 
value of r, will be only slightly greater than 2.4A, 
the value of r at the minimum of the polar 
surface. 

Before considering the nature of the activated 
complex of the reaction it is necessary to consider 
the rate equation. If the transmission coefficient 
is unity and the activation energy zero, the 
theory of absolute reaction rates*! gives for the 
rate expression for a non-linear complex,” 
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(16) 
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From the fact that the activation energy is zero, 
it follows that the hydrogen molecule retains its 
two rotational and one vibrational degrees of 
freedom in the activated state, and that one of 
the vibrational degrees of freedom in the 
initial state becomes replaced by a vibrational 
degree for the complex as a whole. With these 
assumptions the rate equation becomes 


4 
=r] seer (~~—"*) | (17) 


Mams 


the same reduction being here obtained as in the 
case of a reaction between atoms; 7; is the 
separation in the activated state. This equation 
is equivalent to Eq. (9), so that the calculated 


2h? 





value of r; may be compared directly with the 
experimental value of o. 

Magee and Ri take the activated complex to be 
at the top of the rotational-energy barrier, the 
value for r; being 7.64A. However, in view of the 
fact that a stretching of the H—H bond has to 
take place before there can be a transition to the 
polar state, it also appears possible that the 
configuration of highest free energy may be at the 
crossing point, particularly as the entropy is 
lower here. If this is so, r; for the activated com- 


21 For reviews, see H. Eyring, Chem. Rev. 17, 65 (1935); 
Trans. Faraday Soc. 34, 41 (1938); S. Glasstone, K. J. 
Laidler, and H. Eyring, reference 6, p. 184 ff. 

2 Cf. S. Glasstone, K. J. Laidler, and H. Eyring, 
reference 6, p. 309. 
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plex should be only slightly greater than 2.4A. 
The experimental value of 2.7A for o suggests 
that this interpretation is to be preferred. 

The probability of transition to the polar 
surface will be of the order of unity, since the 
interatomic separation is small and the over- 
lapping of orbitals therefore considerable. In 
order that reaction shall proceed to completion 
further extension may be necessary in order that 
the polar surface shall cut the lower surface, and 
from the form of the polar surface it is seen that 
this can take place without expenditure of energy. 
The value of the transmission coefficient depends 
upon the ease with which the system passes from 
the polar surface to the lower surface as compared 
with the probability that it returns to the upper 
surface. In the case of such a reaction as the one 
under consideration the mass point representing 
the system will perform a circular movement 
within the basin representing the polar surface, so 
that the probability of a return to the upper 
state is small; the transmission coefficient of the 
over-all reaction may therefore be taken as unity. 
Physically this means that owing to the presence 
of the additional hydrogen atom the excitation 
energy passes readily into vibrational energy, so 
that the atom has the effect of stabilizing the 
polar ‘‘quenched complex’’ with respect to de- 
composition to give the initial reactants once 
again. 


QUENCHING BY HALOGEN MOLECULES 


The potential-energy surfaces for the quenching 
of sodium radiation by halogen molecules are 
shown schematically in Fig. 7. The curves on the 
left-hand face of the cube, for the molecules 
NaX, are all for '5+ states; these will not interact 
appreciably with the repulsive *2+ states which 
have therefore been omitted for simplicity. The 
curves for X2, on the right-hand face, may readily 
be constructed from the spectroscopic data, 
although this is not necessary here. The curve for 


TABLE II. Calculated quenching radii (A). J(Na) =118.2. 











Halogen E(X) D(X, Te 
Fluorine 95.3 64.6 5.97 
Chlorine 86.5 57.2 5.47 
Bromine 81.5 45.4 5.55 
Iodine 74.3 35.5 5.34 








Fic. 8. Potential- £ 
energy profiles in 
the plane ABCDE 
of Fig. 7, corre- 
sponding to the nor- 
mal X—X _ separa- 
tion. 














Na-X distance oo 


X,- has been constructed on the assumption* 
that the binding energy is one-half of that in Xo, 
and the equilibrium separation the same. The 
designations of the states are indicated in the 
figure. 

The quenching process involves a transition 
from the point B to the point C, and the only 
surface which intersects both of the valleys 
beginning at B and C is that corresponding to the 
polar state Na*+—X-~—X;; reaction will there- 
fore proceed initially in the plane including the 
points A, B and C and running parallel to the 
Na-—X axis. The energy profiles in this plane are 
shown in Fig. 8. The distance between the points 
A and B is readily seen to be I(Na)—E(X) 
+3D(X2) —48.3 kcal., where J(Na) is the ioniza- 
tion potential of sodium, E(X) the electron 
affinity of the atom X, and D is the energy of 
dissociation of the halogen molecule. Neglecting 
exchange forces, as is here permissible owing to 
the high value of r,, the value of r, is given by 


330 . 
= fi. 
I(Na) —E(X)+43D(X:) —48.3 





(18) 


The values of E and D (in kcal.) and of r, (in A) 
for the four halogens are given in Table II. 

Since the activation energy is zero and the 
transmission coefficient probably unity, the rate 
equation again reduces to Eq. (17) ; the values of 
r; may therefore be compared directly with those 
of a. There is here no rotational-energy barrier, 
since the surface is polar in the region in which it 
would appear (c. 5A) and the rapidly varying 
angular momentum (1/r?) cannot introduce a 
maximum into the polar state (1/r).% The acti- 


% Cf. J. L. Magee, J. Chem. Phys. 8, 687 (1940). 
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vated state should therefore be at the crossing 
point, i.e., 7y=e. 

The only data available for quenching by the 
halogens are those of Terenin and Prileshajewa,? 
whose values of o (cf. Table I) are 10 for Br2 and 
6.3 for Iz; the values are only rough but are 
appreciably higher than the calculated. However, 
in view of the approximations involved the 
agreement is not unsatisfactory. In particular the 
calculations explain the observed decrease in rate 
on passing from bromine to iodine; simple colli- 
sion theory would give an increase owing to the 
increased diameter. 

It is clear from Fig. 8 that, as in the case of 
quenching by hydrogen, the reaction is not 
completed in the profile ABCDE, but that ex- 
tension of the X— X bond has to take place until 
the polar and lower surfaces cross; this can be 
effected without expenditure of energy. 


QUENCHING BY HYDROCARBONS 


From the point of view of the efficiency with 
which they quench excited sodium atoms, the 
hydrocarbons appear to fall into two distinct 
classes, corresponding to whether or not an 
unsaturated group is present.2 The saturated 
hydrocarbons are only weakly effective, their 
quenching cross sections being only a fraction of a 
square angstrom, while the unsaturated hydro- 
carbons have quenching cross sections of 40-50A? 
(see Table I). The explanation of these differences 
can only be inferred from a consideration of the 
evidence with regard to the quenching of excited 
mercury and cadmium atoms, which is discussed 
in the last section of Part II (following paper). 
From the fact that triplet mercury and cadmium 
atoms remove hydrogen atoms from saturated 
hydrocarbons, one may conclude that a saturated 
hydrocarbon molecule interacts with an excited 
sodium atom in the same manner as does a 
hydrogen atom, the sodium atom approaching 
one of the hydrogen atoms of the hydrocarbon 
(H —R) to give a complex Na(?P) -H—R, witha 
slight extension of the carbon-hydrogen bond. 
The surface corresponding to this complex is cut 
by a polar surface for Nat — H~—R, and reaction 
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is completed by a transition, probably with a 
considerable extension of the carbon-hydrogen 
bond, from the polar to the lowest surface. The 
explanation of the low quenching efficiency is 
probably that a slight activation energy is 
involved, since there would seem to be little 
possibility of a restriction at either of the 
crossing points. 

The simple saturated hydrocarbons do not 
absorb except at very low wave-lengths; they 
therefore have no singlet levels of excitation 
energy less than the excitation energy of sodium. 
It is thus clear that electronic excitation of the 
hydrocarbons does not take place; probably 
most of the energy passes into vibrational energy, 
with a small amount going into translational and 
rotational degrees of freedom. 

As it does not appear that triplet mercury or 
cadmium atoms in general remove hydrogen 
atoms from unsaturated hydrocarbons it is 
inferred that the interaction, with these atoms 
and with sodium, is not between the metal atom 
and a hydrogen atom, but rather between the 
metal atom and the double bond. The interactions 
can hardly be worked out in detail at the present 
stage, but it would again seem likely that the 
first surface is cut by a polar surface, which is 
probably low-lying since the double-bond system 
readily takes up an additional electron; the 
crossing point of the two surfaces therefore corre- 
sponds to a high value of r,, so that the quenching 
efficiency is high. 

From the fact that ethylene does not absorb at 
wave-lengths greater than about 2100A* it can 
be concluded that there is no electronic energy 
state of ethylene which will combine with the 
ground state (i.e., no singlet state) having an 
energy less than about 135 kcal. Since only a 
singlet state could be produced by collision with 
an excited sodium atom, electronic excitation 
cannot take place, all of the energy going into 
vibrational, rotational and translational degrees 
of freedom, the majority probably going into 
vibrational. 


24 R. B. Mooney and E. B. Ludlam, Trans. Faraday Soc. 
25, 442 (1929). 
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The Mechanism of Processes Initiated by Excited Atoms 


II. Photosensitization by Excited Mercury and Cadmium* 
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The treatment of the preceding paper is extended to the reactions of excited mercury and 
cadmium in their triplet and singlet states. It is shown that with the hydrogen molecule 
quenching of a physical nature is prohibited for the triplet atoms by the necessity for con- 
servation of spin angular momentum, and that dissociation must take place with the initial 
formation of the metal hydrides. In the reactions of Hg *Po, Cd *P;, and Cd *Po these are also 
final products, but with Hg *P; and the singlet excited atoms, dissociation of the hydride into 
the metal and a hydrogen atom takes place at once. Quenching to the metastable state is found 
to be usually inefficient. The atoms react with saturated hydrocarbons by splitting off a 
hydrogen atom, but the unsaturated ones usually interact with the initial formation of an 
excited molecule which when activated by singlet atoms have vibrational energy only, but 
which with triplet atoms are excited to triplet electronic states. The subsequent reactions of the 
excited olefines are discussed with reference to the experimental material. 





INTRODUCTION 


F recent years a considerable amount of 

work has been done on the reactions of ex- 
cited mercury and cadmium.! Mercury possesses 
important resonance lines at 2537A and 1849A, 
but most of the experimental work has been 
done with the former, which corresponds to the 
transition 6'S»)—6*P;. The amount of energy 
liberated when a triplet mercury atom passes 
into its lowest state is 112.2 kcal. per mole, 
which is sufficient to cause dissociation of many 
molecules, particularly as a further amount (8.5 
kcal. per mole in the case of hydrogen) of energy 
may be liberated simultaneously if the quenched 
atom combines with one of the products of the 











TABLE I. 
Energy of Total available 
excited energy if hy- 
atom _ dride is formed 
Metal Line A Transition kcal. kcal. 
Mercury 2537 6'So— 6°P 1 112.2 120.7 
Cadmium 3261 51S9—58P 87.3 102.8 


Mercury 1849 
Cadmium 2288 


6'So—6'P; 153.9 162.4 
5'So—5'Pi 124.4 139.9 








*Contribution No. 1031 from the National Research 
Laboratories, Ottawa, Canada. 

'For reviews, see E. W. R. Steacie, Ann. N. Y. Acad. 
Sci. 41, 187 (1941); W. A. Noyes and P. A. Leighton, 
Photochemistry of Gases (Reinhold Publishing Corporation, 
1941), Chapter V. 
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dissociation. The corresponding energy for triplet 
excited cadmium is somewhat less, 87.3 kcal. 
per mole, but this again may be increased, by 
15.5 kcal., if cadmium hydride is formed. The 
energies liberated by the four excited atoms are 
summarized in Table I. 

Quenching experiments have so far only been 
carried out with mercury? and cadmium’ in 
their triplet states. The data are given in Table I 
of Part I (preceding paper). The chemical re- 
actions which are initiated by the excited atoms 
have been studied, using the usual kinetic 
methods, with mercury and cadmium in both 
the triplet and singlet states, although the 
majority are with the triplet atoms. Most of 
the work has been with hydrogen and the 
hydrocarbons. 

In the remainder of this paper the treatment 
of Part I is extended to the simpler primary 
reactions occurring between excited mercury 
and cadmium atoms and other substances. 


2H. A. Stuart, Zeits. f. Physik 32, 262 (1925); W. A. 
Noyes, J. Am. Chem. Soc. 49, 3100 (1927); A. C. G. 
Mitchell, Zeits. f. Physik 49, 228 (1928); J. R. Bates, Proc. 
Nat. Acad. Sci. 14, 849 (1928); J. Am. Chem. Soc. 52, 3825 
(1930) ; thid. 54, 569 (1932); M. W. Zemansky, Phys. Rev. 
36, 919 (1930); M. G. Evans, J. Chem. Phys. 2, 445 (1934); 
E. W. R. Steacie, Can. J. Research B18, 44 (1940). 

3H. C. Lipson and A. C. G. Mitchell, Phys. Rev. 48, 625 
(1935). 
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REACTIONS OF TRIPLET EXCITED MERCURY 
Reaction with Hydrogen 


The potential-energy surfaces for the reaction 
Hg(*P1)+H.2('2*+,) are shown in Fig. 1, the 
curves on the left-hand face of the cube being 
those for an infinite H—H separation, those on 
the right for an infinite Hg—H separation. The 
reaction Hg(*P1)+H may first be discussed in 
terms of the former curves. 

The curves for the homopolar states may be 
constructed from the spectroscopic data;‘ it 
should be noted that there is no repulsive state 
which dissociates into normal Hg and H. An 
Hg*(2S) ion interacts with an H~(!S) ion with 
the formation of a polar *=+ state, but the 
corresponding curve cannot be constructed owing 
to lack of information from crystal data; how- 
ever, neglecting the exchange interactions it is 
found that the polar curve cuts the *II curve, 
the lowest formed from Hg(*P;) and H(!2*), ata 
separation of about 3A. Since the curve for the 
polar state does not cut that for the lowest 
(?=*) state of the HgH molecule, which is known 
to be the case from the fact that the lowest state 
is homopolar, and since there is no other (re- 
pulsive) state, it is clear that quenching by the 
hydrogen atom can only arise from a jump from 
the polar to the *2* state; as these are probably 
well separated the quenching efficiency is low. 

The surfaces connecting the two sets of curves 
in the interior of the figure may now be dis- 
cussed. If a hydrogen atom (2S) approaches 
linearly an HgH molecule in its lowest state 














Fic. 1. Potential-energy surfaces for the linear system 
Hg—H-—H. The curves arising from the *P2 state of 
mercury are omitted for simplicity, and some of the 
repulsive curves are only continued a short distance. 


*Cf. R. S. Mulliken, Rev. Mod. Phys. 4, 1 (1932). 





LAIDLER 


(?=+) the resulting states are 'Z+ and *Z+, and 
of these, on account of the relationships 


Hg(1\S) +H.2('2+,) =HgH.2(!=+) (1) 
and Hg(}S) +H2(?2*.) =HgH2(32*), (2) 


the 'r+ state arises from the attractive state of 
hydrogen and the *Z+ from the repulsive. The 
curves for HgH(?=*+) and H.('Z*+,) are therefore 
connected by a 'Z*+ surface which has the same 
general form as the lowest surface in the Na+He 
reaction. A second, *=+, surface connects the 
curves for HgH(?2*) and H.(*2+,); this is a 
highly repulsive surface, its lowest point being 
the lowest point on the curve for HgH(?=+) when 
the second hydrogen atom is at infinity. 

The lowest state of HgH which arises from 
Hg(*P,) and a normal hydrogen atom is “II, and 
this interacts with a second hydrogen atom to 
give *II and 'II states. In view of the further 
relationships 


Hg(?P1)+H.('2+,) =HgH2(*l, *z+) = (3) 
and 


Hg(*P1) +H2(*2*.) 
= HgH,(‘Il, 37, TI, 5>+, ee 1y+) (4) 


the *II surface is clearly the lowest, and the only 
one of immediate interest. For the triangular 
configuration it may be assumed that the corre- 
sponding state is slightly attractive. 

There are two comparatively low-lying sur- 
faces representing ionic states of the HgH2 
complex. The lowest state of the H2~ molecule 
is 2E+, and this interacts with a normal (2S) Hgt 
ion to give *Z+ and !Z+ states; similarly the 
lowest ionic state of HgH (2+) interacts with a 
2S hydrogen atom to give *Z+ and '!Z* states. 
Two surfaces therefore connect the HgtH~(?2*) 
and H,~(?+) curves in Fig. 1, of which that for 
the !2+ is lower and corresponds roughly to the 
electron being on the hydrogen atom nearer to 
the mercury atom; the *Z+ state corresponds to 
the electron being on the farther atom. In the 
isosceles-triangular configuration the polar states 
have no such simple significance. 

The potential-energy profiles drawn parallel 
to the Hg—H axis for the linear configuration, 
and corresponding to the normal H—H separa- 
tion in hydrogen, are shown schematically in 
Fig. 2, in which only states corresponding to 
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Hy-H distance oo 


Fic. 2. Potential-energy profiles for HgH2 in the plane 
corresponding to the normal H—H distance. 


excited mercury in the *P,; state are shown. 
It is seen that owing to the high energy of the 
H2- ion at this H—H separation the polar 
curves do not cut the homopolar ones; extension 
of the bond must therefore take place for reaction 
to occur in some section such as that shown in 
Fig. 3. As in the case of the quenching of sodium 
the activated state may be either at the top of 
the rotational-energy barrier, or at a point 
where the homopolar *II surface is cut by the 
polar *2+ one, there being a restriction to crossing 
to the 'Z+ polar surface owing to the change of 
multiplicity. In either case there is a statistical 
factor of } in the rate expression due to the fact 
that the Hg+ and H.~ may approach one another 
initially either on the *Z+ or the *II surface of 
which only the latter leads to rapid reaction. 
The theoretical value of r;-must therefore be 
divided by v2 to be directly comparable with o. 
The value of o is about 2.5A (Table I, Part I) 
so that r; should be about 3.5A. This distance is 
a reasonable one for a crossing point but rather 
smaller than the value at the top of the rota- 
tional-energy barrier (4-5A). Agreement in the 
latter case could, however, be given by a value 
of « of less than one-half. 

The mechanism by which the reaction pro- 
ceeds to the final state may now be discussed, 
and this will be done in terms of the linear con- 
figuration of atoms. It has been seen that the 
polar surface on which the system finds itself is 





Hele +H 024 


Hy (i) Hl 











He-H distance co 


Fic. 3. Potential-energy profiles for HgH» corresponding 
to an extension of the H—H bond. Reaction initially 
proceeds along the #* homopolar surface, followed by a 
crossing to the *2* polar surface. The transition to the *=* 
repulsive surface takes place after an even greater ex- 
tension of the H—H bond has taken place. 


for a *Z+ state, so that a transition to the lowest, 
1>+, state can only take place with difficulty 
owing to the change of multiplicity. Reaction 
therefore occurs more readily by a transition 
from the *Z+ polar surface to the *=* repulsive 
homopolar surface connecting the HgH(?=~) and 
H.(*2+,) curves. At large H—H separations, 
which may readily take place on the polar 
surface, this transition may be effected at a 
point below that corresponding to the first 
crossing point. The subsequent path of the 
system corresponds to its finding its way to 
the lowest part of this repulsive, *=*, surface, 
that is to that part which corresponds to 
HgH(?2*) and H(2S). The quenching of excited 
mercury by hydrogen is therefore of a chemical 
rather than a physical nature, the H—H bond 
being broken.® 

The exact mechanism by which the process 
takes place is, of course, extremely complex and 
could not, at the present stage, be treated 
theoretically. One can be fairly certain, however, 
that a large fraction of the energy set free, 17.7 
kcal., passes into vibrational energy of the HgH 
molecule, the dissociation energy of which is 8.5 
kcal. per mole; the alternative suggestion that it 
becomes relative translational energy implies 


5G. Cario and J. Franck, Zeits. f. Physik 11, 161 (1922). 
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that the mass point proceeds along the lowest 
point of the valley, which would seem improbable 
in view of the complicated form of the surfaces. 
The HgH molecule will therefore decompose 
within the period of one vibration, so that the 
over-all reaction is indistinguishable from® 


Hg+H2=Hg+2H. (5) 


This conclusion is in accord with the findings of 
Olsen,? who showed spectroscopically that HgH 
is not produced in the reaction in the normal 
state, but only excited, presumably by secondary 
reactions. It is interesting to note that in the 
case of the reaction between hydrogen and 
mercury in the metastable (*Po) state the energy 
available for transfer into vibrational energy, 
12.7 kcal., may not be converted efficiently 
enough into vibrational energy to permit dis- 
sociation of the HgH molecule; it might, there- 
fore, be possible to detect HgH in its normal 
state in a system in which this reaction is 
occurring. 

In view of the application of these principles 
that is made below to the reactions of the 
hydrocarbons, the mechanism of the reaction of 
hydrogen may be considered from a slightly 
different point of view. The “‘physical’’ type of 
reaction, 


Hg(*P1)+H2('2+,) =Hg(*So) +H, — (6) 


analogous to that which takes place with sodium, 
is here restricted by the necessity for conserva- 
tion of spin angular momentum (the ‘Wigner 
spin conservation rule’’) unless the hydrogen on 
the right-hand side is in a triplet state; if this 
condition is satisfied, however, the molecule at 
once dissociates, so that the reaction is not of a 
physical nature. There is therefore little possi- 
bility of the hydrogen molecule taking up the 
excitation energy in the form of vibrational, 
rotational or relative translational energy, as it 
was able to do with excited sodium; reaction 
can only proceed effectively by a dissociation of 
the hydrogen. The detailed examination which 
has been made of the mechanism shows the most 


* Cf. J. Franck and H. Sponer, Nachrichten Gottingen 
241 (1928); H. Beutler and E. Rabinowitch, Zeits. f. 
physik. Chemie B8, 403 (1930). 
7L. O. Olsen, J. Chem. Phys. 6, 307 (1938). 
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probable reaction path to be 


Hg(?P;) +H2('S+,) -HgH2 (Il) 
—HgtH2 (@2+)—>HgH2(2*) (7) 
—HeH (2+) +H 2S) H¢('S) +2H(2S). 


It is to be particularly noted that the over-all 
change has taken place without there being, at 
any stage, a restriction due to change of multi- 
plicity, compensating changes in the hydrogen 
molecule taking place to preserve spin angular 
momentum. It will be seen in dealing with the 
reactions of hydrocarbons that arguments of 
this kind are frequently extremely useful in 
excluding mechanisms which might otherwise be 
considered. It may be mentioned that the 
selection rules imposed upon these molecular 
systems by the necessity of conserving spin 
angular momentum are obeyed almost rigorously. 


Quenching to the Metastable State 


Apart from quenching by hydrogen, oxygen 
and the hydrocarbons, the other gases given in 
Table I of Part I quench by transferring the 
excited mercury atom to the metastable state; 
this is known to be true as the presence of the 
metastable atoms can be demonstrated spectro- 
scopically. A potential-energy profile for quench- 
ing by a molecule M to the metastable state, 
which lies 5.03 kcal. below the 6*P; state, is 
shown in Fig. 4. The system first finds itself in 
the dip of the middle curve. If M is an atom the 
quenching process is clearly inefficient, but 
otherwise the dip is a section through a de- 
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Fic. 4. Potential-energy curves for the quenching of 
Hg(#P;) to the metastable state by an atom or molecule M. 
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pression in three or more dimensional space, and 
owing to the possibility of a redistribution of 
energy the complex will not necessarily redis- 
sociate into Hg(*P,:)+M, but may perform a 
series of vibrations until a transition to the lower 
surface is possible. The efficiency of quenching 
to the metastable state is lower than generally 
obtains for quenching to the normal state since, 
owing to the small difference between the energies 
of the *P; and *P» levels, the rate with which 
the complex redissociates into Hg(*P;) and M is 
by no means negligible. The apparent quenching 
efficiency is also lowered in reactions of this kind 
by the fact that *P) mercury atoms may become 
excited to the *P; state by collisions. 

Zemansky® has attempted to show a relation- 
ship between the efficiency of quenching to the 
metastable state and the nearness of a vibrational 
level of the quencher to 5.03 kcal., the energy 
of the transition *P,;—*P». The hydrocarbons 
must, however, be removed from his figure as 
they are now known to quench by a chemical 
mechanism (see below), and when they are 
removed the existence of the relationship can 
hardly be maintained. In any case it is doubtful 
whether such a connection has any theoretical 
justification, as the restriction to crossing caused 
by a failure of levels to fit is probably usually 
only slight. It would appear from the present 
treatment that the problem requires a detailed 
construction of the potential-energy surfaces 
relevant to each individual case, and a calcula- 
tion from them of the course of reaction. This 
could hardly be attempted at the present stage. 


REACTIONS OF TRIPLET EXCITED CADMIUM 


The potential-energy surfaces for the quench- 
ing of the cadmium 3261 line by hydrogen are 
similar in general form to those for the mercury 
reaction shown in Fig. 1. The chief difference 
arises from the fact that, owing to only 87.3 
kcal. of energy being liberated in the transition 
Cd®P,—'S», the process 


Cd(*P1) +H2=Cd('So) +2H (8) 


is endothermic to the extent of 15.6 kcal., and 





8M. W. Zemansky, Phys. Rev. 36, 1919 (1930). 
°Cf. J. L. Magee, W. Shand and H. Eyring, J. Am. 
Chem. Soc. 63, 677 (1941). 
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therefore slow since the activation energy is at 
least this amount. The reaction which forms 
CdH, 

Cd(@P,)+H2=CdH+H, (9) 


is much faster since it is only endothermic by 
c. 0.1 kcal. The quenching reaction is therefore 
largely reaction (9), which has a slight activation 
energy of at least 0.1 kcal. This activation energy 
accounts for the lower value of « as compared 
with the mercury reaction (see Table I, Part I). 

The chief difference between the mechanism 
of the quenching of cadmium and that of mercury 
is that in the former case very little energy can 
pass into vibrational energy of the hydride, any 
that does so having to be paid for as activation 
energy. The hydride CdH will therefore not 
dissociate and should be detectable in its normal 
state; this has been done spectroscopically by 
Bender.'® 

The quenching of cadmium atoms in the 
metastable state will also proceed by a mecha- 
nism analogous to (9), but the activation energy 
will be increased by 1.5 kcal., the difference 
between the *P; and *P, levels, and the process 
is consequently inefficient; this, however, has 
not been investigated, as the levels are separated 
so slightly as to be difficult to distinguish ex- 
perimentally. 

Owing to the comparative slowness of the 
quenching by this mechanism, complications 
may arise in the form of quenching of other 
types; thus the quenching may be partly due to 
a physical type of reaction which took place in 
the quenching of sodium by hydrogen, although 
as was seen for mercury such a reaction is re- 
stricted by a difference in multiplicity between 
the polar and the final state. A certain amount 
of the quenching may also be quenching to the 
metastable state, although this, as has been seen, 
is inefficient. In any case a considerable fraction 
of the total reaction must proceed by the 
chemical process since hydrogen atoms are known 
to be formed in quantity. 

Apart from the hydrocarbons most of the 
other quenching reactions proceed by trans- 
ferring the cadmium atom to the metastable 
state, and, as this lies only 1.5 kcal. below the 


10P, Bender, Phys. Rev. 36, 1535 (1930); cf. L. O. 
Olsen, reference 7. 
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8P, state, quenching is inefficient owing to the 
high probability of the system returning to its 
initial state (cf. Table I, Part I). 


REACTIONS OF THE SINGLET 
EXCITED ATOMS 


Owing to the non-existence of quenching data 
and the scantiness of the kinetic information, the 
reactions of singlet ('P;) excited mercury and 
cadmium will only be treated briefly; surfaces 
may, however, readily be constructed if needed. 
The term relationships required are, for mercury : 


Hg(?P:) +H2('2*+,)-HgH,('5*, 'IT), (10) 
Hg (So) +H2(!2*,)—-HgH,('2*), (11) 
Hg*(?S) +H2- ?2*)-Hg*H~2('5*, *E*), (12) 
and 
HgH ?2+) +H(??S)—HgH,(!5*, *=*), (13) 


with identical states for the cadmium system. 
The possibility of physical quenching, without 
dissociation of the hydrogen, may first be con- 
sidered. It is seen that this may occur via the 
two '=* states of HgH2 and the ‘=~ state of 
Hg*tH.°; reaction (10) first takes place to give 
HgH,(?*), followed by a crossing to the '2+ 
polar state and a second crossing to the lower 
1>+ state of HgHe, which then dissociates into 
Hg(!S) and H2(?Z*,) by the reverse of reaction 
(11). However in both the mercury and cadmium 
reactions the hydrogen so produced would prob- 
ably have sufficient vibrational energy for dis- 
sociation to take place at once, so that the over- 
all reaction will be 


Hg(?Pi) +H2(2*,)Hg(*So) +2H(7S). (14) 


Alternatively, the reaction may be envisaged 
as taking place via HgH(*2*) and H(2S), which 
may arise from the lower 'Z*+ state of HgHe 
(reverse of reaction (13)) ; however the HgH will 
probably have sufficient energy to dissociate. 
There thus seems to be no alternative to assum- 
ing that reaction (14) is the one which takes 
place. 

The quenching cross sections for the atoms in 
the singlet states should be higher than those 
for the triplet states, since the surfaces upon 
which the reaction initially moves are higher, 
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and the separation at the crossing-point to the 
polar state consequently larger; however the 
very small lifetimes of the singlet states make 
the actual magnitudes of quenching relatively 
small. 


PHOTOSENSITIZED REACTIONS 
OF HYDROCARBONS 


Reactions of Saturated Hydrocarbons 


The evidence from recent investigations on 
the reactions of ethane," propane’ and butane" 
photosensitized by triplet mercury and cadmium 
in all cases indicates an initial carbon-hydrogen 
bond split. From the theoretical point of view 
these reactions are very similar to the reactions 
of hydrogen, the normal state of the hydro- 
carbons being a singlet state. Quenching of the 
physical type, the energy being taken up as 
vibrational energy, is again excluded by the 
necessity for conservation of spin angular mo- 
mentum. The possibility of the hydrocarbon 
becoming excited to a triplet state has to be 
considered, but this would probably be unstable 
and decompose to give a hydrogen atom and a 
free radical. The mechanism is probably much 
the same as for hydrogen, reaction proceeding 
via a polar complex. 

The question of whether the hydrides HgH 
and CdH are formed is of interest in connection 
with the strengths of the carbon-hydrogen bonds 
in the hydrocarbons;! if they are, the upper 
limits to the bond strength are given by the 
figures in the last column of Table I, while 
otherwise the upper limits are not so high 
(fourth column). It has been seen that with 
hydrogen the hydrides are initially formed, but 
that if there is sufficient energy available, as in 
the case of triplet mercury, the hydride decom- 
poses at once. Presumably the situation is much 
the same with the saturated hydrocarbons; the 
hydride is formed initially but may decompose. 
In the absence of additional information, there- 
fore, the upper limit to the strength of the bond 


11 E. W. R. Steacie and N. W. F. Phillips, J. Chem. Phys. 
6, 179 (1938); Can. J. Research B16, 303 (1938); E. W. R. 
Steacie and R. Potvin, J. Chem. Phys. 7, 782 (1939). 

2 E,W. R. Steacie and D. J. Dewar, ibid. 8, 571 (1940); 
(i981 R. Steacie, D. J. LeRoy and R. Potvin, ibid. 9, 306 

1941). 
1 E. W. R. Steacie and R. Potvin, ibid. 7, 782 (1939). 
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is that which corresponds to the formation of 
the hydride. If, however, it could be satisfactorily 
proved that the hydride was not formed, the 
smaller figure might be accepted as a maximum. 

Work has not been carried out on the reactions 
of saturated hydrocarbons photosensitized by 
singlet excited atoms, but in the light of the 
treatment in the preceding section there would 
seem to be no reason to doubt that they proceed 
by an initial carbon-hydrogen bond split. 


Reactions of Ethylene 


The reaction of ethylene photosensitized by 
triplet mercury atoms has recently been studied 
by LeRoy and Steacie™ who obtain evidence for 
the following reaction scheme : 


Hg(®P1) +CoH4= Hg(4S0) +CoH«*, — (15) 


CoHy*+C2Hy=2C2H,, (16) 
CoH,* =C2H2+Hz, (17) 
Hg(*P1) +Ho=Hg(4S») +2H, (18) 
H+C.Hy=CoH,, (19) 
2C2Hs=CuHw. (20) 


The most significant facts in favor of this 
mechanism are the initial pressure increase, 
which seems to exclude a free radical mechanism, 
and the increase in acetylene concentration even 
when the total pressure is decreasing. It is of 
interest to consider the first three of these re- 
actions in the light of the present treatment. 

The lowest electronic state of ethylene is a 
singlet (1A,,) state, and since it is known from 
the high quenching cross section of ethylene 
(cf. Table I, Part I) that reaction (15) occurs 
without any restriction, it follows that for con- 
servation of spin angular momentum the excited 
ethylene produced in the reaction must be in a 
triplet state, so that reaction (15) should be 
written as 


He (@P1)+C2H,4(!A1,) 
=Hg(!S))+CsH,(triplet)*. (15a) 


Such a triplet state, of designation *B,,, has 


_'' TD). J. LeRoy and E. W. R. Steacie, J. Chem. Phys. 9, 
829 (1941); cf. H. S. Taylor and J. R. Bates, Proc. Nat. 
Acad. Sci. 12, 714 (1926); J. R. Bates and H. S. Taylor, 
|. Am. Chem. Soc. 49, 2438 (1927). 
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been shown theoretically" to lie about 20 kcal. 
above the normal state, the next higher being so 
high as to play no part here. Thus, of the 112 
kcal. set free when the *P, mercury is deactivated, 
about 20 kcal. passes into electronic energy, 
while the majority of the remainder, about 90 
kcal., goes into vibrational energy, a small 
amount being dissipated as translational energy. 
The excitation of the ethylene molecule to the 
triplet state is accompanied by a twisting of the 
double bond through 90°. 

The second reaction, the deactivation of ex- 
cited ethylene by collision with a normal mole- 
cule, takes place only with difficulty if the 
ethylene molecules produced are both in singlet 
states; however an alternative reaction, 


CoH, (*B,,.)*+CeH,('A1,) 
=CsH,4(*Bi,) +C2H4('Ai,), (16a) 


the removal of vibrational energy from the 
excited ethylene molecule, leaving the molecule 
in its triplet state, will take place with greater 
efficiency, although even here there will be a 
large negative entropy of activation'’® owing to 
the loss of rotation in the activated state.'’ 
Very little polymerization takes place so that 
the reaction 


C.H,(*B,,)*+C2H4('A1,) = CaHs (21) 


is apparently also slow. Since also radiation of 
triplet excited ethylene is restricted by the 
change of multiplicity, it appears that the 
molecule will react chiefly by splitting off hydro- 
gen (reaction (17)). 

Reaction (17) as written will be slow for the 
same reason, but the possibility now arises of 
an alternative reaction, the formation of acety- 
lene in a triplet state: 


C.H,4(B,,)* =CeH2(triplet) +He, (17a) 


which is not restricted by a low transition 
probability, but which has an activation energy 
higher than (17) by the excitation energy of the 
triplet acetylene molecule. The activation energy 
for the formation of unexcited acetylene and 


% J. L. Magee, W. Shand, and H. Eyring, J. Am. Chem 
Soc. 63, 677 (1941). 

16 R, N. Pease, J. Am. Chem. Soc. 53, 613 (1931); H. D 
Burnham and R. N. Pease, ibid. 62, 453 (1940). 
17F, P. Jahn, J. Am. Chem. Soc. 61, 798 (1939). 
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hydrogen from unexcited ethylene has been 
estimated!’ to be 85 keal., so that assuming the 
electronic excitation energy of the triplet acety- 
lene molecule to be the same as that of triplet 
ethylene, or 20 keal. (the actual magnitude is 
unknown), the activation energy for the forma- 
tion of hydrogen and triplet acetylene from 
normal ethylene is about 105 keal. It is thus 
clear that the excited ethylene has derived from 
the mercury atom sufficient energy for reaction 
17a) to proceed without difficulty. 

In the case of the reaction photosensitized by 
triplet cadmium the energy given to the ethylene, 
87 keal., is not sufficient for reaction (17a) to 
proceed without an activation energy. This is 
in accord with the results of Steacie, Potvin and 
LeRoy,!? who find that the acetylene formation 
is only a small fraction of the total reaction, 
which is chiefly polymerization taking place with 
a quantum vield of only about 0.01. 

The experimental facts with regard to the 
reaction photosensitized by singlet excited cad- 
mium atoms are rather more difficult to in- 
terpret.*° Polymerization takes place efficiently, 
but acetylene is a much less important product 
and there is no initial pressure rise. A possible 
explanation would seem to be that two initial 
mechanisms, vibrational excitation of the ethyl- 


8 E. W. R. Steacie, reference 1, p. 195. 

19 E,W. R. Steacie and R. Potvin, Can. J. Research B16, 
337 (1938); B18, 47 (1940); E. W. R. Steacie and D. J. 
LeRoy. unpublished. 

bys E. W. R. Steacie and D. J. LeRoy, J. Chem. Phys. 10, 
22 (1942). 
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ene, Which remains in its lowest clectronic state, 
and the breaking of a carbon-hydrogen bond, 
are here taking place simultaneously. 


Reactions of Higher Olefines 


Preliminary work"! on propylene, butene-1 and 
butene-2 with *P; cadmium shows that poly- 
merization is extremely slow. Quenching data 
are not available. 

Of the possible initial reactions, the splitting 
of a carbon-hydrogen bond and the excitation of 
the double bond to a vibrationally-excited triplet 
state are the most likely. Of the two the second 
probably proceeds more rapidly, as is known 
from the larger quenching cross sections for 
unsaturated as compared with saturated hydro- 
carbons ; however, subsequent reaction is difficult 
as polymerization proceeds only with a large 
negative entropy of activation,'®” and _ triplet 
cadmium does not provide enough energy for the 
splitting off of a hydrogen molecule. 

The author wishes to express his thanks to 
Professor G. Herzberg, Professor R. S. Mulliken 
and Dr. Anthony Turkevich for information 
regarding molecular electronic states, and to 
Professor Henry Eyring, Dr. Samuel Glasstone, 
Dr. D. J. LeRoy and Dr. E. W. R. Steacie for 
valuable discussions and suggestions in connec- 
tion with this and the preceding paper. 


21 E. W. R. Steacie and D. J. LeRoy, unpublished. 

2 W. E. Vaughan, J. Am. Chem. Soc. 56, 3863 (1932); 
G. B. Kistiakowsky and J. R. Lacher, ibid. 58, 123 (1936) ; 
J. B. Harkness, G. B. Kistiakowsky, and W. H. Mears, J. 
Chem. Phys. 5, 682 (1937). 
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" Thermodynamics of High Polymer Solutions' 
PauLt J. Fiory 
Isso Laboratories of the Standard Oil Development Company, Chemical Division, Linden, New Jersey 
(Received October 15, 1941) 
id A. statistical mechanical treatment of high polymer — classical values, which accounts for the large deviations of 
de solutions has been carried out on the basis of an idealized — high polymer solutions from ‘‘ideal’’ behavior. The entropy 
ta model, originally proposed by Meyer, which is analogous — of disorientation of a perfectly arranged linear polymer is 
to the one ordinarily assumed in the derivation of the found to be of the order of R cal. per chain segment. After 
1g “ideal” solution laws for molecules of equal size. There is introducing a suitable heat of mixing term, partial molal 
j obtained for the entropy of mixing of nm solvent and free energies are computed, and the calculations are com- 
of N linear polymer molecules (originally disoriented), pared with experimental data for all concentrations. Phase 
et AS= —k((n/B) In 01+ N In v2] where v; and v2 are volume — equilibria have been calculated in the region of partial 
id fractions and 6 is the number of solvent molecules re- miscibility. The theory predicts, in agreement with experi- 
mn placeable by a freely orienting segment of the polymer ment, that the critical composition for partial miscibility 
or chain. This expression is similar in form to the classical lies at a low concentration of polymer. Low intrinsic 
“i expression for equal-sized molecules, mole fractions having viscosities of polymers dissolved in poor solvents are 
been replaced by volume fractions. When the disparity attributed to the tendency for the molecule to assume a 
ult between the sizes of the two components is great, this more compact configuration in such an environment. 
ge expression gives entropies differing widely from the 
el 
he ~ 
INTRODUCTION Meyer and co-workers** have found that AS, 
* T is customary to correlate the thermodynamic for solutions of long chain compounds dissolved 
a properties of binary liquid systems with the ” simple solvents, e.8-, for carbon tetrachloride 
a so-called ‘‘ideal’’ solution laws. These rest funda- solutions of long chain esters, is much gene 
oo mentally on an entropy of mixing given by than Eq. (2) would predict. In other words, they 
‘8 have found strong negative deviations from 
wal ASiixing= —R(n1 In X1+n2 In X2), (1) Raoult’s law, even when AH,=0. Meyer and 
wi Liihdemann® suggested that this abnormal 
where m, and m2 are the numbers of moles of the entropy effect was due to the ability of sections of 
[WO COMpPONENtS, and X, and Xs are their mole 4 Jong chain solute molecule to act more or less 
a fractions in the mixture. Partial differentiation of independently of one another as kinetic units. 
6): (1) with respect to m vields for the partial molal Thus, the effective mole fraction should be 
J. entropy of dilution greater than the actual mole fraction. An equiva- 
AS,=—Rin Xi. (2) lent explanation has been advanced recently by 
Powell, Clark and Eyring.’ Apparently they 
Although (1) and (2) frequently have been ap- were unaware of the serious criticisms of this 
plied indiscriminately, at finite concentrations point of view advanced by Hiickel,? principally 
they can be justified theoretically only when on the grounds that it inherently assumes pri- 





molecules of the two species are equivalent in 
size and shape, i.e., when molecules of the two 
species are interchangeable in the liquid. The 
importance of this limitation has been empha- 
sized by a number of authors during recent 
vears.2—4 


* Presented before the Division of Physical and Inorganic 
Chemistry of the American Chemical Society, at the 
Atlantic City Meeting, September 11, 1941. 

> E. Hiickel, Zeits. f. Elektrochemie 42, 753 (1936). 

E. A. Guggenheim, Trans. Faraday Soc. 33, 151 (1937). 

*R. H. Fowler and G. S. Rushbrooke, Trans. Faraday 
Soc. 33, 1272 (1937), 


mary validity of the ideal solution laws when 
applied to mixtures of molecules differing greatly 
in size. Hiickel has also pointed out that the 
kinetic units cannot be truly independent, nor are 
they independently interchangeable in the liquid 
as applicability of (1) or (2) would require. 
Finally, in order to explain experimental results, 





5 KX. H. Meyer and R. Liihdemann, Helv. Chim. Acta 18, 
307 (1935). 

6 Ch. G. Boissonnas, Helv. Chim. Acta 20, 768 (1937). 

7R. E. Powell, C. R. Clark, and H. Eyring, J. Chem, 
Phys. 9, 268 (1941). 
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it is necessary to assume a size for the kinetic 
unit, which, though equal to the entire solute 
molecule at infinite dilution, decreases rapidly 
with increasing concentration, approaching a 
limiting size of some 15 to 30 chain atoms®’ at 
high concentrations. 

Recently Meyer’ has provided a more plausible 
explanation of the ‘‘abnormalities’’ of solutions of 
long chain, or polymeric, molecules. He has 
suggested an idealized model for such solutions, 
which consists of a quasi-solid lattice, not unlike 
that customarily used in the derivation of (1) for 
the case of a mixture of equal-sized molecules. A 
cell of the hypothetical lattice may be occupied 
either by a solvent molecule or by a segment of 
the polymer (solute). The segments of a given 
polymer chain must occupy a continuous se- 
quence of adjacent cells, which, however, may 
meander irregularly through the lattice. There 
will be many such sequences emanating from a 
given cell, and the polymer molecule may con- 
form to any one of them through suitable 
rotations about the valence bonds of the chain 
skeleton. According to Meyer, entropies of 
mixing in excess of ‘‘ideal’’ values given by (1) 
are the result of these numerous configurations 
available to the polymer molecules. 

Meyer has discussed his model in qualitative 
terms only. Recently quantitative statistical 
mechanical treatments*" of this model have 
been carried out independently by Huggins and 
the writer. On the basis of the preliminary 
accounts which have been published concerning 
these investigations, the results appear to be in 
substantial agreement. In this paper the writer’s 
treatment of the problem will be presented in 
full, with particular emphasis on an extension of 
the model and applications to hitherto unex- 
plained characteristics of polymer solutions. 


DERIVATION OF THE ENTROPY OF MIXING 


Throughout this paper polymer molecules 
which are linear, i.e., consist of a single linear 
sequence of structural units, are under con- 








8 K. H. Meyer, Zeits. f. physik. Chemie B44, 383 (1939); 
Helv. Chim. Acta 23, 1063 (1940). 

9M. L. Huggins, J. Chem. Phys. 9, 440 (1941). 

10 Paper presented by Dr. Huggins before the Wilder D. 
— Colloid Symposium, Ithaca, New York, June 20, 
1941. 

" P. J. Flory, J. Chem. Phys. 9, 660 (1941). 
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sideration. Somewhat arbitrarily, the symmetry 
number of these molecules is taken to be two, due 
to (assumed) indistinguishability of the two ends 
of the molecule. Concurrently with the derivation 
of the fundamental relationships, special promi- 
nence will be given to the assumptions involved. 

Assumption 1: the mixture of polymer and 
solvent molecules is assumed to conform to 
Meyer’s model. This involves assumption of (a) a 
quasi-solid lattice in the liquid, (b) interchange- 
ability of segments (not necessarily identical 
with the polymer structural units) of the polymer 
chain with solvent molecules in the lattice cells 
and (c) independence of lattice constants on 
composition. Although these assumptions obvi- 
ously are quite artificial, there is some consolation 
in the fact that analogous ones are inherent in the 
model from which Eq. (1) is derived for mixtures 
of small molecules of equal size. 

Assumption 2: all polymer molecules are as- 
sumed to be of the same size. 

Assumption 3: the average concentration of 
polymer segments in cells adjacent to cells 
unoccupied by the polymeric solute is taken to be 
equal to the over-all average concentration. 

Let m)= the total number of cells in the lattice. 

N=the number of polymer molecules, each 
composed of x segments, any one of which may 
replace a solvent molecule in the lattice. 

=a coordination number for the lattice, i.e., 
the number of “‘first neighbor”’ cells in the lattice 
available to the next consecutive segment of the 
polymer chain. 

In order to determine the total number of 
possible non-identical configurations, let us first 
consider the process of successive addition of 
polymer molecules to the empty lattice con- 
taining mp» cells. After N polymer molecules have 
been added, there will be 7» —xWN cells which may 
be occupied by the terminal segment of the next 
polymer molecule to enter the lattice. According 
to assumption 3, the expected number of avail- 
able unoccupied cells immediately adjacent to 
the cell occupied by the terminal segment of this 
polymer chain will be y(m)>—xN)/mo. For the 
third segment the expected number of available 
cells will be given by a=(y—1)(mo—xN)/no, 
since the cell occupied by the terminal segment 
eliminates one of the neighboring cells. 
Assumption 4: the expected number of avail- 











al 
er 
Is 


on 


on 
he 


es 
iS- 


of 


Ils 


ce. 
ch 
ay 


e., 
ice 


‘he 
of 


of 
on- 
ive 
lay 
ext 
ing 
ail- 

to 
this 
the 
ible 
/no, 
ent 


rail- 





THERMODYNAMICS OF HIGH POLYMER SOLUTIONS 53 


able positions for each successive segment is 
taken to be equal to a=(y—1)("»—xN)/no. This 
is of the nature of an approximation, as it fails to 
exclude impossible configurations in which two 
segments, separated by two or more intervening 
segments but belonging to the same chain, 
occupy the same cell. Use of this ‘‘assumption”’ 
obviously will lead to computation of too many 
configurations. 

For the expected number of configurations 
which the particular chain may assume without 
shifting the terminal group we have 


[y/(y-1) Jaz (3) 


and for the expected total number of configura- 
tions for the chain in the entire lattice 


yv4i=23(mo—*N)[y/(y—-1) Jar", — (4) 
23 (mo—xXN)*(y—1)7-!, mo'*, (4.1) 


where the factor } has been included in order to 
eliminate redundancy of configurations due to 
the indistinguishability of the ends of the polymer 
chain. 

The total number of possible configurations for 
the system consisting of N polymer and 1 solvent 
molecules in a lattice composed of mp=n+xN 
cells is given by 


W=(1/N!) TT vs (5) 


where the factor 1/N! eliminates redundant 
configurations which differ only by an interchange 
of one or more pairs of polymer molecules. 
Substituting (4.1) in (5) we obtain by suitable 
rearrangement , 


(2-1) N xy (ny/ x)! ! F 
C2) Gkceed 
no 2NN!/L(no/x—N)! 


Replacing factorials by Stirling’s approximation, 
N!=(N/e)* 


y-1 (2—-1)N (n +x N)*t 
w-(~ ) ie. 
n"N® 





For the entropy of mixing, according to the well- 
known Boltzmann relationship S=k In W, where 


k is Boltzmann’s constant, we obtain"! 


nN N 
AS mixing ™ -A|n In ( )+nin (- )| 
n+xN n+xN 


+k(x—1)N[In (y—1)-—1]—kRNIn 2. (7) 








This expression refers to the entropy change in 
the process of mixing ” pure solvent and N 
perfectly arranged (or oriented) polymer mole- 
cules, subject, of course, to the assumptions 
introduced above concerning the nature of the 
solution. When x=1, (7) reduces to the classical 
expression (1), except for the symmetry term. 

A few remarks should be made concerning the 
use in the above derivation of ‘‘expected values” 
for the numbers of configurations, since actual 
values in many cases may deviate enormously 
from expectancy. When a is greater than unity 
(i.e., when the concentration of polymer is not 
too great) the expected number of available 
configurations which extend from a given site in 
the lattice, is, according to Eq. (3), a very large 
number. (We consider x to be at least of the 
order of 100.) But the actual number of configu- 
rations may be zero; for example, there might be 
only one cell available for the second segment and 
none for the third. On the other hand, for a more 
favorably located initial cell, the number of 
available configurations may be an exceedingly 
large number. Since the total number of sites 
available to the terminal segment of the Nth 
polymer molecule to be added to the lattice is 
very large when a>1, the actual total number of 
configurations available to the particular mole- 
cule will deviate negligibly from the expected 
number calculated from (4). 

When a is appreciably less than unity, i.e., at 
high concentrations of polymer, the expected 
number of configurations (Eq. (3)) extending 
from a given site will be very much less than 
unity. If @ is sufficiently small, e.g., of the order 
of 0.5, then vy also may be extremely small (for 
large x). That is, at high polymer concentrations 
the expected number of configurations available 
for an additional polymer molecule is much less 
than unity. Of all the possible configurations for 
the N—1 polymer molecules in the lattice, only a 
small fraction (not exactly equal to vy) possess 
one Or more continuous sequences of x vacant 
cells any place in the entire lattice. According to 
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this interpretation, vy is the average ‘‘yield”’ of 
such sequences of unoccupied cells per configura- 
tion of the N—1 polymer molecules in the lattice. 
Thus, after a certain concentration is exceeded, 
addition of more polymer molecules to the lattice 
(possessing a fixed number of cells) decreases the 
total number of configurations for the system.'** 
In spite of the use of expected values from which 
actual deviations may be large, Eq. (6) for the 
total number of configurations should be a valid 
approximation at all concentrations, subject only 
to the original assumptions. 


THE ENTROPY OF DISORIENTATION 


The ‘‘configurational” entropy of a lattice filled 
with N polymer molecules in the absence of 
solvent, which is obtained by setting »=0 in (7), 
is!! 


AS*=RN In (x/2) 
+kN(x—1)[In (y—1)—1]. (8) 


When x is large 
AS*/ Nx 2k In [(y—1)/e]. (9) 


That is, the entropy change per mole of segments 
is of the order of R.”> These relationships express 
the entropy change for the transformation of NV 
polymer molecules from a_ state of perfect 
orientation and arrangement in the lattice to the 
state of random entanglement. This probably is 
the source of an appreciable portion of the 
entropy of fusion of polymer crystallites (e.g., in 
rubber, polyesters, etc.). Obviously, other terms 
arising from change in the lattice constants and 
in the arrangements and motions within the seg- 
ments when fusion occurs also contribute to the 
total entropy of fusion. 

It can be shown that the lattice “saturated” 
with polymer molecules will contain a few vacant 
cells, equivalent to the ‘‘holes’’” in simple liquids. 
From the point of view of the equations used 


1% Alternatively, one may consider that the polymer 
molecules are continually diffusing in the lattice, i.e., the 
system is continually changing from one configuration to 
another. Then vy is the time average number of available 
sequences of x units. 

2b H. Mark, J. App. Phys. 12, 41 (1941), using the 
statistical theory of the internal configurations of polymer 
molecules, has arrived at a similar value for the entropy of 
disorientation per chain bond about which there is free 
rotation. 

1H. Eyring, J. Chem. Phys. 4, 283 (1936); H. Eyring 
and J. O. Hirschfelder, J. Phys. Chem. 41, 249 (1937). 
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above, these vacant cells should be treated as 
molecules of a second component. However, their 
concentration is so small that no significant error 
has been committed in setting »=0 in the 
derivation of (8). 

The entropy change on mixing 7 solvent with 
N randomly entangled (liquid state) polymer 
molecules is 


Budtutes = ASn ixing — AS* 


n 
= =A In (—*.) 
n+xN 


smuin(—_)]|. a 
n+xN 


(10.1) 


Or 
AS mixing = —k[n In vw+N In v2 |, 


where 7; and v2 are volume fractions of solvent 
and polymer, respectively. The pure liquid 
polymer (disoriented) will be chosen as the 
standard state for the solute, as is customary for 
binary liquid mixtures. Equation (10) will serve 
as the starting point for the derivation of the 
thermodynamic relationships discussed below. 

It is interesting to note that (10.1) is so similar 
in form to (1); replacement of volume fractions 
by mole fractions in the logarithms of the former 
yields the latter. The entropy of mixing of 
polymer and solvent is greater than (1) would 
predict, but less than that for mixing the com- 
pletely dissociated polymer segments with sol- 
vent. Meyer’s® qualitative discussion of the 
problem leads to the same conclusion. 


PARTIAL MOLAL FUNCTIONS 


Differentiating (10) with respect to , there is 
obtained for the partial molal entropy of the 
solvent 


AS,=—R In (1—ve) —R(1—1/x)v2 (11) 
= Rvo(1/x+02/2+02/3+---). (11.1) 

Similarly, 
AS.=—R In (1—0;)+R(x—1)1. (12) 
Inasmuch as the forces acting between two 
polymer molecules, or between a polymer and a 
simple molecule, obviously are of the saine sort as 


those operating between pairs of simple mole- 
cules, the heat of interaction between a polymer 
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molecule and its environment should vary 
linearly with the composition. Consequently, the 
heat of mixing can be represented by a relation- 
ship used by Scatchard" for simple molecules 


AH=BV,VanN/(nVitNV2), (13) 


where B is a constant for a given pair of liquids 
and V,; and V2 are the molar volumes of solvent 
and polymeric solute, respectively. According to 
the model considered here, V2=xV. 
The partial molal heat contents are, according 
to (13) 
AH, = B V2" (14) 
and 
AH.=BV w= BxV 2. (15) 


Combining these equations with (11) and (12), 
there is obtained for the partial molal free energy 
of the solvent 


AF, = RT In (1 —Ve) 
+RT(1-1 x)vo+ BV we" (16) 
= — RT>?,[ 1, x+(1 — K)v» 2 
+02/3+---], (16.1) 
where 
K=2BV,/RT 
and for the partial molal free energy of the solute 


AF.=RT In (1-0) —RT(x—1)01 
+B Vix0/". (17) 
Equations (16.1) and (11.1) reduce to the 
“ideal” solution relationships at infinite dilution. 


This is virtually a necessary consequence® of the 
Gibbs-Duhem relationship 


JAF ,/d In X;=dAF,2/d In Xo. 


Since for any acceptable model for the system the 
activity a2 of the solute must be proportional to 
the concentration at infinite dilution 


8 In as/d In X2=(1/RT)dAF,2/d In X2=1 
dAF,/d In X,=RT. 
Or 
AF, = RT In Xi 


which at infinite dilution may be replaced by 
AF, = —RTX2= — RTv2/x, (16.2) 


which is identical with (16.1) as v2 approaches 
zero, 


4G. Scatchard, Chem. Rev. 8, 321 (1931). 
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THE EFFECT OF SEGMENT SIZE 


It has been assumed that a section of the 
polymer chain which occupies a volume equal to 
that of one solvent molecule is equally free to 
occupy any one of y—1 cells, entirely inde- 
pendent of the configuration of preceding seg- 
ments of the chain. If the solvent molecule is 
very small, or if there are limitations, perhaps 
steric in nature, on the tortuosity which the chain 
may comfortably assume, then the size of the 
section of the chain which is free to orient itself 
at random, irrespective of the orientation of the 
preceding sections, may be equivalent to the 
volume occupied by several solvent molecules. In 
this situation, we may construct a new lattice in 
which one cell will accommodate either 8 solvent 
molecules or one segment of the polymer chain. 

The entropy of mixing according to this revised 
model will be given by either of the Eqs. (7) or 
(10), provided that 7 is replaced by n’=n/8; x 
likewise has been decreased by the same factor, 
but a modified symbol is unnecessary. Differ- 
entiating the resulting expression with respect to 
n (not n’)"! 


AS,=-— (R/B)[In (1 —ve) + (1 —1/x)v2 | (18) 
= (Rv2/B)(1/x+v2/2+02?/3+---), (18.1) 


and 


AF ,=(RT/8)[In (1 —v) 
+(1—1/x)ve+Kv-?/2 | (19) 


= —(RTv2/8)(1/x+(1—K)v2/2 
+02/3+---)], (19.1) 


where K is re-defined as 28BV,/RT. Equation 
(19.1) reduces to AF, = —RTX; at infinite dilu- 
tion, as is required. Equation (12) for the partial 
molal entropy of the solute can be applied intact 
when 6>1, except that x must be taken equal to 
the number of re-defined segments per chain. 
Consequently, Eq. (17) can be written 


AF,=RT[In (1 —v}) 
—(x—1)o14+Kxv;/2]. (17.1) 
COMPARISONS WITH EXPERIMENTAL DATA 


Osmotic pressures of dilute polymer solutions 
have been investigated quite extensively,'*-'’ as a 


1 H. Staudinger and G. V. Schulz, Ber. 68, 2320, 2336 


(1935); G. V. Schulz, Zeits. f. physik. Chemie A176, 317 
(1936). 
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means for determination of molecular weight. 
Strong deviation from Van't Hoft's law, 1.e., from 
“ideality,”’ are invariably observed at low con- 
centrations, necessitating a frequently precarious 
extrapolation to infinite dilution. Experimental 
results for a number of dilute polymer solu- 
tions!S 16 1%20 support the simple relationship be- 
tween Osmotic pressure m and concentration c of 
polymer (in grams per cc) 


r=RTc M,+0c. (20) 


where VM, is the number average molecular 
weight of the polymer and 0 is a constant, inde- 
pendent of M,,)*!* representing the degree of 
departure from ideality. 

For dilute solutions Eq. (19.1) vields 


r= —AP, VERT? 2 V\8x+RT7(1-—K)v2" 218 


which reduces at once to (20) 3! overlooking the 
difference between volume and weight concen- 
tration in the second term 


b=RT\1-—K) 2148. (21) 


Huggins®'® has reported a similar deduction of 
(20) for the case K=0, 3=1. 

When various solvents or solvent mixtures are 
considered, according to (21) the greater the heat 
of mixing the smaller will be the value of } at a 
given temperature, since K depends directly on 
the heat of mixing constant B. As K approaches 
unity 5 approaches zero, and the solution ap- 
proaches classical “ideality.” When K becomes 
but slightly greater than unity, (x large), the 
region of partial miscibility is reached and the 
polymer will no longer dissolve to an appreciable 
concentration (cf. seg.). Thus, a dilute high 
polymer solution should be most nearly “‘ideal” 
when it is approximately on the verge of pre- 
cipitation. Osmotic pressure measurements on 
rubber dissolved in benzene-alcohol mixtures” 


%K. H. Meyer, E. Wolff and Ch. G. Boissonnas, Helv. 
Chim. Acta 23, 430 (1940). 

17 E, Wolff, Helv. Chim. Acta, 23, 439 (1940). 

%Q. Hagger and A. J. A. Van der Wijk, Helv. Chim. 
Acta, 23, 484 (1940). 

1*See also, H. Mark, Physical Chemistry of High 
Polymeric Systems (Interscience Publishers, Inc., New 
York, 1940), pp. 228-245. 

. G. Gee, Trans. Faraday Soc. 36, 1161 (1940). 

* Actually, the partial molal volume V; should be used 
in place of V:, but since the difference between them is so 
small in dilute solutions (see ref. (16)) they need not be 
distinguished here 
% G. Gee, Trans. Faraday Soc. 36, 1171 (1940). 
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and on solutions of polyvinyl! chloride in tetra- 
hydrofurane, a good solvent, and in dioxane, a 
poor solvent,” confirm this deduction.** However, 
to conclude, as Gee* has done, that a better 
extrapolation to infinite dilution is obtained in 
the presence of alcohol is not necessarily correct ; 
higher terms in the series expansion of (19) must 
be included if the extrapolation is to be improved. 

Although the statistical theory presented above 
predicts the correct form for the osmotic pressure- 
concentration relationship (20), quantitative 
correlation with the observed ‘deviation term”’ 
(i.e., bc?) requires a larger value of 8 than would 
be expected from steric considerations regarding 
chain flexibility. Thus, from osmotic pressure 
measurements on dilute solutions of rubber in 
toluene, Meyer, Wolff and Boissonnas found 
1,5=56 cal. at 298°K. The heat of mixing, 
deduced from the temperature coefficient of 
the osmotic pressure, was negligibly positive.*° 
Taking K=0, according to (21) bV:=300/8, or 
8225.3, which corresponds to about seven 
isoprene units, or 28 chain atoms, per segment. A 
similar value for )V, was obtained by Wolff" for 
gutta-percha in toluene. 

For the system toluene-polystyrene, according 
to osmotic pressure measurements by Schulz,!® 
bV,=39 cal. at ordinary temperature. Although 
reliable data concerning the heat of mixing do not 
appear to have been obtained,”® certainly the 
heat effect accompanying mixing of two such 
chemically similar materials is negligible, and it 
may be assumed that K=0. A value of B=8 is 








23H. Staudinger and J. Schneiders, Ann. d. Chemie 541, 
151 (1939). 

24 R. E. Powell, C. R. Clark and H. Eyring, in a paper 
presented at the Conference on Viscosity sponsored by the 
New York Academy of Sciences, February 14, 1941, have 
reached the fantastic conclusion from Gee’s results that a 
benzene-methanol mixture is a better solvent for rubber 
hydrocarbon than benzene alone. Apparently they have 
confused the literal meaning of the term ‘‘ideal” with its 
conventional meaning with reference to the so-called 
“ideal” solution laws. Addition of methanol to the benzene 
solution of rubber decreases —AF,, the free energy of 
dilution, thus decreasing the thermodynamic solvent 
power. Since —AF; is greater for pure benzene solution- 
than — AF; for an “ideal” solution, addition of methanol 
decreases the discrepancy between observed and ‘“‘ideal’’ 
values, but this is another matter. 

2 1.. Hock and H. Schmidt, Kautschuk 10, 33 (1934), 
have shown by direct calorimetric measurements that the 
heat of mixing of benzene and rubber is positive, but ver) 
small. 

% See Ch. G. Boissonnas and K. H. Meyer, Zeits. f. 
physik. Chemie B44, 392 (1939). 
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Fic. 1. Partial molal free energy as a function of the 
volume fraction of polymer for 8 = 1 and 8=5.3. e—Osmotic 
pressure results for toluene-rubber, reference 16; O—re- 
sults from vapor pressures of toluene-rubber mixtures, 
reference 16; (_}—results from vapor pressures of benzene- 
rubber mixtures, reference 28. 


required in order to bring the theoretically 
derived relationship into agreement with experi- 
ment; the corresponding segment contains about 
16 chain atoms. 

If rotation about each single bond in the chain 
were completely free, one or two single bonds 
should be sufficient to allow a new choice of 
direction of the chain independent of its preceding 
portions. Allowances for the limitations on the 
effective available configurations due to the 
potential restricting rotation about single carbon- 
carbon bonds and to the steric interferences of 
side chain substituents, presumably would raise 
the segment size somewhat. But it is difficult to 


account for segments of the size indicated by the 


above analysis on these grounds alone. Over- 
simplifications inherent in the model which has 
been postulated, and in particular assumption 4, 
probably are responsible for a considerable part 
of the discrepancy. The latter assumption allows 
too many configurations to be included, thus 
giving a configurational entropy which is too 
large. 

The segment sizes calculated above are similar 
to those deduced from the temperature coeffi- 





~~ 


cients of polymer viscosities,?” but the writer is 
not inclined to consider this as more than a 
coincidence. Vicsous flow isa kinetic phenomenon, 
whereas the thermodynamic properties of polymer 
solutions (taking the pure liquids as reference 
states) depend primarily upon configurational 
characteristics, according to the premises of this 
paper. Unless it can be shown that kinetic and 
configurational segments should be the same, 
correlation’ of these independently estimated 
segment sizes is likely to be illusory. 

In Fig. 1 partial molal free energies calculated 
from Eq. (19) are compared with experi- 
mental values for solutions of rubber in toluene 
(Meyer, Wolff and Boissonnas'*) and in benzene 
(Stamberger**). The dilute solution data, which 
are from those discussed previously, were ob- 
tained osmotically. Data at higher concentrations 
(v2>0.20) are based on vapor pressure measure- 
ments.”° Stamberger’’ found no perceptible differ- 
ences in the vapor pressure lowering for samples 
of rubber varying considerably in molecular 
weight. This is in agreement with Eq. (19) which 
predicts virtual independence of AF, on x when x 
is large and the solution is not too dilute. 

The upper curve in the figure has been 
calculated from (19) taking K=0, 8=1 and 
x= 2900, which corresponds to the average mo- 
lecular weight, 270,000, obtained by Meyer and 
co-workers'* from the limiting value of 2/c. For 
the lower curve a value of 8=5.3 has been chosen 
in order to obtain agreement with experiment at 
low concentrations (see above). The value of x 
has been lowered correspondingly to 550 in order 
to preserve consistency with the observed mo- 
lecular weight. With increasing concentration the 
experimental points depart from the 8=5.3 
curve, approaching the 8=1 curve at high 
concentrations. The comparison of observed and 
calculated values at high concentrations of 
polymer is shown more effectively in Table I by a 
tabulation of observed activities a;= ;/p1° and 
those calculated from (19) and the relation- 
ship In a,=AF, ‘RT, taking 8=1. 


27P. J. Flory, J. Am. Chem. Soc. 62, 1057 (1940); W. 
Kauzman and H. Eyring, ibid. 62, 3113 (1940). 

*8 P. Stamberger, J. Chem. Soc. 2318 (1929). 

29 Concentrations for the experimental data actually are 
in weight fractions, although volume fractions are indi- 
cated. The difference is unimportant, however, since the 
densities of the two components are similar. 
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The measure of agreement in this range is all that 
should be expected considering the crudity of the 
original assumptions. 

It would be possible, of course, to improve the 
agreement between theory and experiment by 
assuming an appropriate dependence of 8 on 22. 
However, this apparent variation of 8 may 
reflect imperfections in the theory, rather than 
an actual variation in segment size with con- 
centration. 


PHASE EQUILIBRIA 


When there is a large heat effect accompanying 
the mixing of two liquids composed of molecules 
of ordinary size, the observed entropy of mixing 
usually deviates considerably from the ideal 
entropy given by (1), i.e., the solution is not 
“regular.’”*° Due to differences in the energies of 
interaction between neighboring pairs of like and 
unlike molecules, distribution of the two species 
is not random, as is assumed in the derivation 
of (1). Furthermore, if AH <0 the comparatively 
large solvent solute interactions in the solution 
are likely to lead to a higher degree of local 
orientation of neighboring pairs than is found in 
the pure components, thus decreasing the 
entropy of mixing; if AH>0, the opposite is 
likely to be true. 

Similarly, in the case of polymer solutions for 
which AH is large in magnitude the entropy of 
mixing should not be expected to agree with that 
calculated assuming random distribution of the 
two species, random configurations for the 
polymer molecules, and packing and local orien- 
tations equivalent to the pure liquid com- 
ponents.!* ‘‘Deviations’’ should be expected to 
arise from the same sources discussed above for 
simple molecules. When A4H>0, moreover, com- 
pact configurations of the polymer chains should 











TABLE I. 
v2 ai (obs.)*6 a; (calc.) 
9.95 0.20 0.13 
0.90 0.36 0.25 
0.80 0.57 0.45 
0.70 0.73 0.60 
0.60 0.86 0.73 











80 J. H. Hildebrand, Solubility of Non-Electrolytes (Chem- 
ical Catalog Company, New York, 1936), second edition, 
p. 65. 
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be preferred, thus decreasing the configurational 
entropy (cf. seq.). 

In spite of these additional limitations imposed 
on the partial molal entropy and free energy 
Eqs. (18) and (19) when A/7>>0, it is of interest to 
examine the character of the phase-composition 
diagram which they predict when AI// is suffi- 
ciently large to cause separation into two phases. 

Equilibrium between two phases in a binary 
system requires that their compositions v2 and 
v’» shall be such that 


AF \(v2) = AF \(v'» ; AF .(v) = AF 2(v'2). (22) 
When A, =0, i.e., when K = 28BV,/RT=0, AP, 


decreases continuously as v2 increases from zero 
to unity. But when K is sufficiently positive, 
according to (19) the curve for AF, will possess 
a maximum and a minimum in this range. 
Differentiating (19) with respect to v2 and 
equating to zero 


Kxve — (Kx —x+1)v2+1=0, 
the solution of which is 


Kx—(x—1)+[(Kx—x+1)*—4Kx]? 
= . (23) 
2Kx 





The critical condition for incipient separation 
into two phases requires that the maximum and 
minimum in AF, coincide, i.e., that the two roots 
of (23) shall be equal. This occurs when™ 


K (critical) = (1++/x)?/x. (24) 
Substituting (24) in (23) 
Ve (critical) =1/(1+/x). (25)" 


When «x is large the critical value of K is 
slightly greater than unity, and it approaches 
unity as x approaches infinity ; at larger values of 
K two phases will co-exist in equilibrium. The 
critical composition, at which the two phases 
become identical, according to (25), occurs at a 
low concentration when x is large. The same 
relationships could have been derived by setting 
dAF./dv2 equal to zero, since 
Bxv, 


JAF ,/d02= -_ ( )eaP/av. (26) 


V2 


In Fig. 2, ~BAF\/RT calculated from (19) 
when x=1000 is ‘plotted against v2 for several 
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values of K in the neighborhood of K (critical) 
= 1.0644. Two scales have been used in order to 
cover in a single graph a broader concentration 
range. The character of the curves changes very 
rapidly with K. Since K depends inversely on the 
absolute temperature, these curves may be re- 
garded as representative of a series of tempera- 
tures in the vicinity of the critical temperature 
for total miscibility. The corresponding curves for 
—AF,/RT show related characteristics, as is 
required by (26); the occurrence of a maximum 
in AF, is accompanied by a minimum in AF, and 
vice versa. 

When K exceeds the critical value, there exists 
a pair of values v2 and v’2 which satisfies the 
equilibrium conditions (22). The smaller of these 
(v2) will lie between zero and the maximum in 
—AR, plotted against vo, as in Fig. 2; the larger 
(v2) will lie at a concentration greater than that 
at the minimum in —AF;. The value of AF;(v2) 
= AF ,(v'2) must necessarily be small, since AF, at 
the maximum in the curve is small. In fact, 
except when K is very near K(critical), it is 
permissible to set AF \(v’s) =0 and to compute v’2 
from K using Eq. (19); —dAF,/d22 is so large at 
v’2 that the error so introduced is insignificant. 

Since (19) can be solved explicitly for K but 
not for v’s, it is more convenient, having set 
AF,=0 to compute K for a given value of 2's. 
That is, we may take 


K=-—2[In (1—0'2) +(1—1/x)0’2]/(v'2)?. (27) 
Substituting (17.1) in the second part of (22) 


(In v2) /x—(1—1/x)(1—v2) +K(1 —v2)?/2 
= (In v’2)/x—(1—1/x)(1—v’2) 
+K(1—v'2)?/2, (28) 


which can be solved by trial for v2 using the above 
values of vw’, and K. If necessary, a second 
approximation can be computed by substituting 
the first approximation values of K and v2 into 
(19) to obtain AF, which in combination with v’2 
will yield a second approximation value for K. 
The second approximation value of v2 can then be 
obtained from (28). 

In the range where this procedure is applicable 
v2 is quite small (assuming a large value of x). 
Replacing (1—v2) and (1—v»)? by unity in (28), 
substitution therein of (27) for K, and expansion 




















Fic. 2. Partial molal free energies as functions of the 
volume fraction of polymer for various values of K in the 
critical region for partial miscibility. 


in series gives as an approximate relationship 
between v2 and v's 


(1/x) In (v'2/v2) =2/x+(v'2)?/2-3 
+ (v’o)3-2/3-4+(v'2)!-3/4-54+---. (29) 


These procedures are unsatisfactory very near 
the critical point, e.g., for K <1.10 when x = 1000. 
Here it is necessary to solve for K and ?-¢ by trial. 

The curve in Fig. 3, showing v2 and v’, as a 
function of 1/K when x=1000, has been calcu- 
lated as described above. Since the ordinate is 
proportional to the absolute temperature, the 
curve represents the calculated temperature- 
composition diagram. The concentration (v2) of 
polymer in the solvent phase rapidly approaches 
a negligible value as 1/K decreases. This is 
amplified by the following calculated composi- 
tions which extend beyond the range of Fig. 3. 
See Table II. 

Thus, the theory predicts an extremely unsym- 
metrical phase composition relationship. The 
solvent is readily soluble in the polymer even at 
large values of K (e.g., far below the critical 
temperature for complete miscibility), but the 
polymer will dissolve perceptibly in the solvent 
only very near the critical point, and then only to 
a small extent." 

These deductions are confirmed by the work of 
Brénsted and Volqvartz* on mixtures of various 
lauric acid esters with high molecular weight 


tJ. N. Brénsted and K. Volqvartz, Trans. Faraday Soc. 
35, 576 (1939) ; ibid. 36, 619 (1940). 
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Fic. 3. Phase composition curve calculated for a solution 
of a polymer composed of 1000 segments. @—z-propyl 
laurate-polystyrene, reference 31; O—chloroform-cellulose 
acetate, reference 32. 


(average ca. 300,000) polystyrene, which had 
been fractionated to remove low molecular 
weight constituents. Below the critical tempera- 
ture for complete miscibility, none of these were 
observed to dissolve the polymer to an appreci- 
able extent, although the polymer imbibed as 
much as several times its volume of solvent. 
Brénsted and Volqvartz concluded that their 
results were ‘“‘compatible with the existence of a 
critical point at which the solubility of solvent in 
the swelled phase is infinite.’’ Their results are 
also consistent with a critical value of v2 slightly 
greater than zero, which would be in accord with 
theoretical predictions. 

Brénsted and Volqvartz’s data for the n-propyl 
laurate-polystyrene system are shown in Fig. 3. 
The temperature scale, shown on the right, has 
been related to 1/K by an arbitrary propor- 
tionality factor (340) so that the absolute temper- 
ature is proportional to 1/K, as the definition of 
K requires. Results of Papkov, Rogovin, and 
Kargin® on the chloroform-cellulose acetate 
(average molecular weight about 175,000)* 
system are also shown using the same tempera- 
ture scale. 

Quantitative agreement between theory and 
experiment cannot be expected in view of the 

%S. Papkov, S. Rogovin, and V. Kargin, Acta Physi- 
cochim. U.R.S.S. 8, 647 (1938). 

33 This estimated molecular weight is based on a specific 
viscosity given by the above authors, and the intrinsic 


viscosity-molecular weight data of E. O. Kraemer, Ind. 
Eng. Chem. 30, 1200 (1938). 
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various approximations which have been intro- 
duced. But the ability of the theory to predict 
the two outstanding characteristics of such 
binary systems, namely, a very low critical 
composition, and a high swelling capacity ac- 
companied by negligible solubility, is a most 
gratifying confirmation of the validity of Meyer’s 
model and the statistical treatment of it which 
has been presented here. Many of the so-called 
gels composed of polymers saturated with solvent 
doubtless owe their existence merely to this 
peculiarity of the temperature-composition curve, 
rather than to any peculiar colloidal structure, as 
generally has been assumed. 


MOLECULAR CONFIGURATION IN 
POOR SOLVENTS 


Staudinger and Heuer** observed that the 
relative viscosity of a dilute solution of poly- 
styrene was less in a poor solvent than in a good 
one at the same concentration. They also found 
that the addition of a precipitant, e.g., ethanol, to 
a solution of polystyrene in a solvent such as 
toluene reduced the relative viscosity continu- 
ously to the point of precipitation. Staudinger 
and Schneiders* found lower viscosities for 
solutions of polyvinyl chloride in dioxane than in 
tetrahydrofurane, which their osmotic pressure 
results showed to be a better solvent. Similar 
results were obtained by Gee” on solutions of 
rubber in various solvents and solvent mixtures. 

The diminution of viscosities in poor solvents 
frequently has been attributed to decreased 
solvation. If by ‘‘solvation’’ some specific inter- 
action between solute and solvent is meant, it is 
difficult to conceive of solvation in a toluene 
solution of polystyrene, where the van der Waals 
forces between like and unlike components are 
both similar and small. Furthermore, if solvation 
does occur, how can it be so markedly decreased 











TABLE II. 

1/K v2 v’2 
0.862 1.4107 0.1915 
0.833 1.7X1077 0.2357 
0.770 1.3 107" 0.3300 
0.715 1 X<1072! 0.407 
0.625 5 x<10-* 0.527 
0.500 3 x<10-*% 0.683 








#4 H. Staudinger and W. Heuer, Zeits. f. physik. Chemie 
A171, 129 (1934). 
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by the addition of a few cc of alcohol which has no 
particular affinity for either component? 

When the heat of mixing of polymer and 
solvent is positive, i.e., in a poor solvent medium, 
the polymer molecules should tend to prefer more 
compact configurations where intramolecular 
contacts are more numerous. As is well known, 
globular or spherical polymer molecules enhance 
the viscosity of their solutions relatively little in 
contrast to the marked effects of uncoiled long 
polymer chains which assume more or less 
random configurations. Therefore, the low rela- 
tive viscosities in poor solvents very probably 
reflect changes in average configuration of the 
polymer molecules, dependent upon the solvent 
medium. Intrinsic viscosities at incipient pre- 
cipitation are of the order of one-half or less of 
the values for the same polymers in a good 
solvent medium.” This would lead to the con- 
clusion that the polymer molecule at the critical 
point pervades about the same volume as another 
molecule of half its chain length in a good solvent. 
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CONCLUSIONS 


Quantitative statistical treatment of Meyer’s 
model accounts for the following hitherto unex- 
plained peculiar properties of high polymer solu- 
tions: (a) the shape of the partial molal free 
energy-composition diagram throughout the com- 
position range (i.e., the large deviation from 
“ideal” entropies), (b) the virtual independence 
of partial pressure on molecular weight of the 
polymer except in dilute solutions, and (c) the 
extreme dissymmetry of the mutual solubility- 
temperature diagram and the limited high 
swelling capacity of linear polymers in certain 
media. These achievements confirm the correct- 
ness of the general point of view regarding these 
systems. On the other hand, quantitative agree- 
ment between calculations and observations is 
not particularly good, presumably due to im- 
perfections of the model itself and to limitations 
of the present treatment imposed by assump- 
tion 4. 
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This is a continuation of an effort to reduce to the solution of a characteristic value problem 
the rigorous calculation of thermodynamic properties of systems in which the intermolecular 
forces are sufficiently short ranged so that practically the entire potential energy of the system 
results from interactions between nearest neighbors. The partition function of such a system 
can be expressed in terms of the largest characteristic value of a linear operator equation and 
finally as the ratio of partition functions of systems with relatively few particles whose potential 
energy functions differ somewhat from those in the original system. A method of evaluating 
grand partition functions and one of introducing interactions between more distant neighbors 
are discussed. The general theory is applied to the calculation of magnetization, internal 
energy, and specific heat of two-dimensional ferromagnets on the basis of the Ising model. 
There seems to exist a \-point phase transition in the change from ferromagnetic to nonferro- 


magnetic states. 


I. INTRODUCTION 


HE intermolecular forces in many solids 
are sufficiently short ranged to permit the 


*Sterling Research Fellow in Chemistry, Yale Uni- 
tr Present address, Cornell University, Ithaca, New 
ork, 


total potential energy of the solid to be repre- 
sented by the sum of the interactions of each 
molecule with its nearest neighbors. With this 
approximation the evaluation of the potential 
energy contribution to the partition function 
and the thermodynamic properties of a solid can 
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be reduced to the solution of a characteristic 
value problem.'* Inasmuch as such problems 
have been much studied while the summation 
of highly multiple sums has been a rather 
neglected subject, the reduction is an advan- 
tageous one. 

Suppose a system can be divided into LX M 
XN cells so that there will be Z layers, each 
composed of NXM cells and each layer will 
have the same average behavior as any other 
layer. The molecules in each cell may have 
several or even a continuum of possible con- 
figurations. Since the same is true of each layer, 
let us define a set {a;} of all those configurations 
available to a physical coordinate a@ of a layer. 
For example, in a binary substitutional alloy a; 
might represent an ith arrangement of the two 
kinds of atoms in a layer. 

If v(am) is the total potential energy with 
respect to a coordinate a between all neighboring 
cells of a layer in the mth configuration, and if 








MONTROLL 


V(Qm, @n) is the total potential of interaction 
between cells of one layer in the mth configuration 
with their nearest neighbors in the mth con- 
figuration, the potential energy factor of the 
partition function with respect to the coordinate 
a is* 

Z=>-:- ¥ exp —OLv(ai)+---+2(az) 


{a1} fap} 


+0(a1, a2)+--++0(ars,a,)]. (1) 
Making the substitutions 


V(a:, aj) =v(a;)/2+0(a;, aj) +0(a;)/2); 
6=1/kT (2) 
(1) becomes 


L-1 
Z=>:--> II exp —O0V (ai, ais). (3) 
far} fag} 
If {&,} is the set of orthonormal characteristic 
vectors and {A,} the corresponding set of 
characteristic values of the matrix equation 


Ay =exp [ -6V]-¢ (4) 
Wa) | exp — 6V (a1, a) exp —6V(a, a2) ve (a) ) 
ny oe =i\cxp — 6V (ae, a) exp — 6V (ae, a2) See (az) 





LJ 


(when the exp —6V matrix is symmetrical), by 
expanding exp —@V(a;, a,) as a bilinear combi- 
nation of {¥,} it can be shown that! 


L-1 , 
Z= AmaxL 2 Yinex (ot) F (5) 
ay 
where Amaxis the largest characteristic value of (4). 


II. PARTITION FUNCTION BY METHOD OF 
SUCCESSIVE APPROXIMATIONS 


The secular equation of the matrix in (4) 
corresponding to a real crystal is of too great an 
order to be solved directly, so it seems necessary 
to seek a somewhat indirect method of solution 
that will lead to correct results. In studying 
long, narrow (three and four atoms wide), two- 
dimensional ferromagnetic strips! it was found 
that the method of successive approximations 
leads to asymptotic expressions for Amax which 
show excellent agreement with the exact values 

!E. Montrol!l, J. Chem. Phys. 9, 706 (1941). Henceforth 
this paper will be referred to as SM-I. 

2H. Kramers and G. Wannier, Phys. Rev. 60, 252 and 


Had (1941). These papers will be designated by KW-I and 
<W-IT. 





as calculated numerically. This suggests that 
the same approach might give good results when 
applied to real crystals. 

In the method of successive approximations one 
first guesses a reasonable vector to and calculates 


ti“=(exp -V6)-o, 
V1 “(ay) = y Yo(ae) exp — 6V (ay, a). 


Then some value of a, say a1=0, is chosen and 
one sets! 
¥i(0) =yo(0) = 1, 
V1 (a1) =Y"(ay)/~i"(0) 
= Yo(a2) exp —6V (a1, ax) / 
> Yo(a2) exp —OV(0, ay). 


1.€., 


3 Note the implicit assumption of independence of the 
various kinds of coordinates. For example, in a ferro- 
magnetic system this would be equivalent to the assumption 
that spin orientations are independent of lattice vibrations. 

* The superscript u is used to imply that ¥:"(0) has not 
yet been normalized to unity. For a detailed mathematica! 
discussion of the method of successive approximations sec 
Frazier, Duncan, and Collar, Elementary Matrices (Cam- 
bridge University Press, 1938), p. 140. 
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¥i(a1) is the second approximation of ¥(a,), and 
to find the (+1)st approximation one proceeds 
according to the general scheme 


¥i"(a1) = 2 Wi-1(a2) exp — OV (ax, a), 6) 
Wi(ar) =i" (ar) /pi"(0). 

The corresponding characteristic value is 

i= Vilar Wilae) exp —6V(a1, a2)/ 
- Evlavslas). (7) 


Using this procedure we shall now develop a 





power series in @ for Ainax ANd Wmax(a) of a real 
crystal. Since this development is the appropriate 
one for high temperatures, the logical first 
approximation wt) of ¢ should be that which 
corresponds to a system in a state of complete 
disorder; that is, all configurations of a given 
layer should be chosen as being equally probable. 
Such a & (unnormalized) is 


to=(1, - 1, “7% 1), 
1.€., Wo(a) =1 for all a. 


Then if there are v possible configurations per 
lattice point and v™ per layer 


¥"(a1) = Zz 1 [1 —6V(a, a2) + V?(ay, a2) /2— 7? | 


{a2} 


=v™N—6)) Var, a2) +6 2) V(r, a2)/2—--- 


{a2} 


and 
¥i(a1) =~"(a1) /"(0) 


{a2} 


=1-—0y-"" ¥ [V(a1, a2) — V(O, a2) J +0" * © [LV2(an, ae) — V2(0, ae) ]/2 





{a2} {a2} ~ VIO, as)[ Von, a») — V(0, a2) J} ose, 

Note that 
Yr(ax) —Yo(ar) = 0(6). 
Now 
. DX Vilar) {1—OV (ai, a2) +6 V?(a1, a2) /2— +--+} 
W2(a1) =Y2"(ar)/P2"(0) = 7 ‘ 
S DX vila2) {1—OV(0, a2) +6V2(0, a2)/2—-++} 
So 
¥(a) ~y2(a1) =1—0-¥% I [ Va, a) = V(0, a2) +0 y-N { pm [V?(a, a2) — V2(0, a2) |, 2 
{a2} {ao} 
+y?¥N bY [V(a1, a2) — V(0, ae) JL V (a2, 3) — V(O, a3) — Vas, a4) J} —-++. (8) 


{a2, ag, a4} 


Here we see that 


¥2(a1) —i(a1) = 0(6). 


The coefficients of higher powers of @ are rather long expressions so we will not record them here. 
However, by continuing the process started above one can show, using the same sort of reasoning 


as that used in the appendix, that 


¥i(a) — Pis(a) = 0(6'). (9) 
vi" (a) = »» ¥i-1(a2) exp —OV(a, ae) 
= ."(0) Pi" (ar) /Wi"(0) =Yi(ar)pr"(0) 


and ¥;(a) is equal to ¥:_1(a@) to within terms of order 6’, ¥,"(0) is equal to A;_1 (the /th approximation 


Since by definition 


tO Amax) to within terms of order 6’. Thus 


— yMN 
ho=v ’ 


h=v™N—y-NO SY V(ay,02), 


{aj, a2} 
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) ee ~dho= pMNn —_ p-N@ > V(ay, a2) + Py-¥N f } V?(a1, aa) /2 


{a1, ag} {a}, a} 


+o" p sing >» V(an, @2) V(az, a3) a aes | a V(a, a2) | : (10) 


{a1, a2, a3} {a}, a2} 


At this point let us search for a physical interpretation of the iteration process described above. 
If the successive approximations are made without bothering to expand in a power series in 6, 


Yo(a) =1 for all a. 


Then 
¥1(a) = y exp —6V(ai, a2)/2 exp —6V(0, a2), 
Y2(ay) = y Vilas) exp —O0V(ay, a2) / 2 Yil(az) exp — OV (0, az) 
= =. exp — OL V(a1, a2) + V(as, a3) | ‘as exp — 6 V(0, a2) + V(ase, a) ]. 
In general ; ; 


D+: Dexp —OLV (a1, a2) + Va, a3) +> +++ Var, ar41) ] 
DL: Lexp —OLV (a1, a2) + V(as, a3) +++ +V (arr, a) 
Wmax(a) — pila) = 0(6"*?). (12) 
The characteristic value corresponding to ¥(a) is 
Ar=t*(exp -OV)i/t.*- 
D-:+dD exp —OLV (a1, a2) +--+ + V (a2, ar41)] Ser 
"So gee wil Maa *<-4 Veal Sa 


Vi(a1) 





(11) 





but in the appendix it is shown that 
So141/ Se1= Si41/S1+0(8). 
Therefore, remembering the remarks made following Eq. (9), 
max = Sr41/Si+0(0"4). (13) 
S,;= a > exp — 6 V(a1, a2) +++-+V(ais, ai) |] (14) 


Now 


is the partition function of a system / layers thick in which the interaction between lattice points 
in the first (and in the /th) layer is one-half that between those of other layers of the system. Thus 
as Lo, the (L—1)st root of Z (see Eq. (5)), the partition function of a system L layers thick, 
is equal, to within terms of O[(1/k7)]'*+!, to the partition function of a system of (+1) layers 
divided by that of a system of / layers: 


ZULD = §),4/S,+0(0). (15) 


Since S; is the partition function of a system 1X MXN, 


Si=[u(2) max YOLEN ol, B) max, (16) 
where yu(2) max is the largest characteristic value of 
u(L)e(1, Bi) =X (1, Be) exp — Wi(B1, Bo)? (17a) 
B2 


with 


W1(Bs, Be) = wi(B1)/2+wi(B1, Be) +wi(Be)/2, 








a) 


b) 
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w,(8;) = potential energy of interaction between neighboring elements in the ith 1X WN layer in 
configuration £;, 


and 


w1(B:, 8;) = potential energy of interaction between elements in the ith /X N layer in configuration 
8;, and their nearest neighbors in the jth layer in configuration 8;. 


In the energy term W/,(8;, B;) the interaction between the N elements in the first row of the ith 
layer and their nearest neighbors in the first row of the jth layer (likewise for interactions between 
Ith rows of the two layers) is reduced by half. 
If we let 
Dm, 1= 2 ° oe exp — OL W1(81, Bo) ++ >» + Wi(Bn—1,m) J= Tims (18a) 
1 m 


by the same arguments as those used in the derivation of (13): 
M(L) max = Tm, 1/ Tm, 1+0(0"**) (18b) 
Try, 41/T 1, 41 +0(8") 
T 141, 1/T1,1.+0(6"") 
= Try1,41/T1,1+0(0"*"). 


and 





Z'L-D (M1) = 


(19) 


Now 77, is exactly the partition function of a system of 1X1 XN lattice points in which the potential 
energy of interaction between all pairs of lattice points is the same as that between pairs in the 
original system, except: (1) for pairs both of whose components lie on one of the four faces of width 
l and length N—the interaction between such pairs in one-half that between normal pairs; (2) for 
pairs both of whose components lie on one of the four edges (formed by the intersection of two of 
the 1X N planes) of length N—the interaction between such pairs is one-fourth that between normal 
pairs. This means that the 1X1 N system can be ‘‘compressed”’ in the N direction by the character- 
istic value method, and if w(/, 2) max is the largest characteristic value of the operator equation 


w(l,1)x(1, 131) =X x(I, 1; v2) exp — 0X1,i(M1, Y2), (20) 
Ti, ' aes [w(I, Dua) 2 x(L, ie "¥) mex} (21) 


Neglecting 1 as compared with L, M, N 


ZUEMN = w@(1+1, 1+1)max/o(l, 1) max +0(0"). (22) 
In Eq. (20) 
Xi, (V1, Y2) = X71, 1) /2+ x7, (M1, V2) +X7, (2) /2. 


x1,.(7:) =total potential energy of interaction between elements in the first layer while it is in 
configuration y;. The interaction between pairs of elements on the border of the /X/ square 
is only one-half that between other pairs of elements. 

x1,1(¥1, ¥2) =total potential energy of interaction between elements in the first layer in configura- 
tion 7; and their nearest neighbors in the second in configuration y2. The interaction between 
the four elements on the corners of the /X/ square in layer 1 and their nearest neighbors on 
the corners of layer 2 in one-fourth that between pairs of internal neighbors. The interaction 
between elements on the 1X2 faces is one-half that between pairs of internal elements. 


Thus we have shown that in order to calculate the partition function correct to terms of order 
(1/kT)!' of a system of LX MXN cells (as L, M, N=) it is sufficient to find the largest character- 
istic value of an operator equation whose kernel involves the potential energy of a system of 1X1X2 
cells. 
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The low temperature calculation of Amax can also be effected by the use of the method of successive 
approximations. Suppose there are two possible ordered configurations for a given layer® and that 
these are represented by 0; and 0». As a first approximation to ¢ we choose a vector which defines 
a completely ordered state : 


Yo(0e) = Wo(01) =1 ; Yo(a) =0 if 0:4aXF¥0o; Wi" (a) =exp [- 0 V(0;, a) ]+exp [- 6V (0c, a1) |. 


Here we shall keep 





Then ¥o(01) =yi(0;) = --- =1. 
betes _exp [—6V(0;, a1) ]+exp [—0V (02, a) ] 
exp [—6V(0;, 0:) ]+exp [—0@V(0y, 02) ] 
p> [exp —OV (a2, 0:)+exp —OV (a2, 02) |Lexp — OV (ai, a2) ] 
¥2(a1) = . 





> [exp —@V (ae, 0:) +exp — OV (ae, 02) |[exp —0V(01, a2) ] 


a2 


By continuing the approximation process we finally obtain 


Zz, nial Zz. }exp vient OV (an, 0:)++exp —60V (an, 02) } {exp — OL Va, a) + ated hy + V(an-1, Qn) }} 





— (23a) 
X- +L fexp —OV (an, 0:) exp —OV (an, 02)} fexp —O[V (01, a2) +++» + V(an—1, an) ]} 
and 
dn= E ¥n(ar)Yn(a) exp — OV (ar, a2)/ Xo Yuleu)pn (a2). (23b) 


It is apparent from Eqs. (23a) and (23b) that the low temperature expression for Ajax is the quotien t 
of the partition functions of two smaller systems. Proceeding in a manner analogous to that used 
for the high temperature expressions, one can finally evaluate Z as the ratio of characteristic values 
of operator equations whose kernels are functions of the potential energy between two small layers 
of lattice points. The details of this calculation will not be carried out here, but they will be apparent 
from the discussion of the special case of two-dimensional ferromagnetism at the end of this paper. 


III. GRAND PARTITION FUNCTION AND THE CHARACTERISTIC VALUE METHOD 


In systems of more than one component it is usually necessary or at least it is usually more con- 
venient to use the grand partition function for the calculation of thermodynamic quantities. So, we 
shall now investigate the applicability of the characteristic value method to the evaluation of grand 
partition functions. 

The grand partition function® of a k component system is defined as the sum 

P= 3 :°-Dd exp (uimitpomet+ : + + une — Ey.) 0 (24) 
N1, N2, +++, Nk, Mm 
over-all compositions 1, m2, --+, m, and over-all energy states E,,. The terms yp; are the chemical 
potentials of the various components and are so chosen that 


1 
N;=— > -- +> nn exp O(mim+-:+-+uin,.—EL,,), (25) 
P Alek 


equals the number of particles of the ith component of the system. 


_§ This would be the case if each lattice point could achieve either of two configurations. The method can without 
difficulty be extended to the case of more than two ordered configurations per layer. 

* For discussions of grand partition functions and their uses cf. Tolman, Principles of Statistical Mechanics (Oxford 
University Press, 1938), p. 619; Gibhs, Collected Works (Longmans, Green and Company, 1931), p. 187. 
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In a nearest neighbor system of the type discussed in Part I, retaining the notation introduced 
there, 


P= a > exp O[ uimi(ar) + + + + + xms(ar) —0(a1) + mi (ae) + + + + + em (a2) —v( a2) 


al 


+ pms (ar) ++ +++ pm (ar) —v(a,) —v(a1, a2) —0(a2, a3) — +++ —v(ar-1, a7) | 


with 1,(a;) being defined as the number of particles of kind 7 in the jth layer when it is in configu- 
ration a@;. 

If we make the substitution 
V(a;, aj) = —[pmim(ai) +++ + + ym(a;) —v(a;) ]/2+0(a;, a) —Lrmi(a;) + + + + um(a;)—v(a;) ]/2, (26) 


L— 


Pat... 


1 
al aL i=l 


exp — Ov(aiai+:) 
and if the exp — @V matrix is symmetrical the problem is formally the same as that of the evaluation 
of Z 
L-1 . . 
P = uaz! Do Winax(a) |. (2 7) 
The calculating devices introduced in studying Z are valid for studying P. 


As an example let us consider the grand partition function of a one-dimensional binary substitu- 
tional alloy. Let there be No particles of kind “‘0”’ and N; particles of kind ‘‘1,”’ and further let £;; 


sé 99 


be the energy of interaction of a particle of kind ‘“z’’ with a nearest neighbor of kind ‘‘j.”” Since 
any lattice point can be occupied by either kind of particle, 


No+N1-1 . 
P=Xmax oo” Wmax(a) P 
and the exp — 0V matrix is 
exp —6V(0,0) exp —é@V(0, ®) 


exp -ov=( 
exp —@V(1,0) exp —é@V(1, 1) 


- ( exp O(40 — Eno) exp O[ (wom) /2 ) 
exp OL (uo m1) /2— Ev ] exp 6(4i— En) ) 


By simple algebra the roots of the secular equation of exp —®V are 
= 3 {exp 6( uo — Eon) +exp 6( ui — Eu.) +(Lexp (40 — Eo) —exp 6(ui — Eu) P+4 exp 6(uit wo 2E12)) Y ’ 


ho = 5 {exp 0(4o — Eo) +exp O(u:1 — En) — (Lexp 6( 40 — Eno) —exp 6(u1—En1) ? +4 exp 0(ui+ Mo — 2Ew)) if 
Thus 


PUNotNi-) =), 
and the chemical potentials uo, 4; must be chosen to make simultaneously 
No=4 log di /0(u08), 
Ni =02 log \y/0(u18). 


Inasmuch as this problem was introduced merely as an example, we shall carry it no further. It 
might be mentioned however that there is no phase transition in this linear system. 
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IV. INTRODUCTION OF INTERACTIONS WITH 
MORE DISTANT NEIGHBORS AS APPLIED 
TO LINEAR FERROMAGNETS 


The characteristic value theory of the partition 
function developed in SM-I and in further 
studies of the theory as reported so far in this 
paper are adapted to ‘‘systems in which the 
molecular forces are sufficiently short ranged to 
permit the total energy of the system to be 
represented by the sum of the interactions of 
each molecule with its nearest neighbors.” One 
might now ask if it is possible to include direct 
interactions between more distant neighbors 
(as well as the correlation effects which are 
automatically taken into account in the present 
method) and still retain some remnants of the 
characteristic value theory. Formally it is made 
possible by placing more than one molecule in 
each of the cells mentioned in the introduction. 
The potential energy terms in (1) would then 
represent interactions between multiply popu- 
lated cells. Of course the exp —@V matrices 
would be of a higher order than those of the 
corresponding nearest neighbor systems. A 
further difficulty arises in that the exp —6V 
matrix of the extended problem may not be 
symmetrical even though that of the original 
problem is, so we cannot blindly use (5). 

In order to understand an approach that is 
satisfactory for evaluating some partition func- 
tions in which interactions between more distant 
neighbors have been introduced, and in order to 
see what effect the introduction has on thermo- 
dynamic properties let us consider in detail the 
statistics of a linear ferromagnetic chain. 


The seat of the magnetic properties of a ferro- © 


magnetic element is in the unpaired electrons of 
the d shell of the atoms at each lattice point. 
We shall let —e be the energy of interaction of 
two neighboring atoms whose unpaired spins are 
parallel and +e that when they are antiparallel. 
Furthermore, we shall let —f be the energy of 
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interaction of two ‘‘next nearest neighbors” with 
parallel excess spins, and f that with antiparallel 
spins. Also, we write J=e/kT and K=f/kT. 
A linear chain of N lattice points can be divided 
into N/2 cells containing two lattice points per 
cell 

5 2.2 8.3.6 N-1 N 
x xX xX|X._ xX |---| kK XI. 

I II Ill N/2 

By representing a lattice point with excess+spin 
by a “0” and one with excess—spin by a “‘1,”’ 
there exists a unique two digit number, a, in 
the binary system which represents a given cell 
configuration. If we define v(¢;, ¢;) as the spin 
interaction between the ith lattice point in the 
configuration ¢o; and the jth lattice point in 
configuration oj, the partition function of our 
linear system is 


Z=) ae p > exp — O[v(o1, a2) +v(01, o3) 
+v(02, ¢3)+0(02, 01) +++ +0(on, ov-r) ]. 
Then substituting 


x 


V(ai, Qi41) =v(02i-3, 02;-2)/2+(o2:-3, 2;-1) 
+0(02:~2, o2:-1) +0(o2i-1, o2:) 
+0(02:-1, o2;)/2, 
N/2-1 


Z=): =a > II exp —6V(ai, O41). 
a1 an/2 i=2 


Now let us examine the unsymmetrical matrix 
a2\a, 00 01 11 10 


OO { evt2k e! e72k ev) 

O1 et =e 2 +2K e! e7 2k 
exp —6V= 

11 e7 2k ed ert +2K e! 

10 { e/ e7 2k en J =e 2J+2K J 








To demonstrate the calculation of the matrix 
elements let us consider the element in the 01th 
row and the 11th column. This is the term 
exp —6V(01, 11) and corresponds to the con- 
figuration 


0 1|/1 1 
123 4 
Since 
v(o1,02)=€, (o1,03)=f, v(o2,03)=—€, v(o2,0s)=—f, (03, 04)=—€, 
V(01, 11) =¢«/2+(f—e—f) —e/2= —e= — JkT 
and 


exp —6V(01, 11) =e’. 
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1/J 


Fic. 1. Average internal energy of spin per lattice point 
(in units of e) for a nearest and next nearest neighbor 
system. 


The matrix 
(1 0 1 0) 








0 1 0 —-1 
transforms exp — 6V into 


( e +2K a e7 2k e! + e~ 0 0 
R exp (— 0V)R“'= 
0 0 erJ+2k stp e7 2k e/ iste e~/ 
q 1) 0 e~/ —e/ e724 42K _ 92K 








and the roots of the symmetrical upper left factor of the secular equation of R exp (—6V)R™- are 
Mi =e? +e?* cosh 2J+2 cosh J(e** sinh? J+1)}, 


ho =e"*X +¢?* cosh 2J—2 cosh J(e** sinh? J+1)?. 
The vectors ‘ 
i*=(A, By, 0, 0), $2* = (Ao, Bz, 0, 0), 
with 
Ne — en 2 +2K — e-2K 





[d.2— 20, (e-2/+2K 4 @-2K) 4 (el fe-J) 24 (e-24 42K 4 2K) 273 , 
B,=(1-A,*)!, 


are the pair of orthonormal characteristic vectors of R exp (— @V)R-, whose characteristic values 
are \; and Xo, respectively. The corresponding characteristic vectors of exp — 6V are 


ti* = o:*R=(Aj, By, Ai, Bi)/V2,  te*=o2*R=(Ao, Bo, Ao, Be), v2, 


also with characteristic values \; and )o. 
At finite temperatures >0, the roots \3 and , of the secular equation of the lower right sub-matrix 
of R exp (— 6V)R™ have absolute values less than \,. Since their parent matrix is not symmetrical, 
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Qa) ae An ay a2 An 

0 xX X x 0 BA 0 KX K vss K 1 yi 

0 XxX X x 0 B OX XK -+ K 1 2 
N 0 x x x al B's N 0 x x x 1 1's 

0 KX K «s+ K O By 0 K Moves K 1 yy 

O KX K «s+ K O By 0 KX K «s+ K 1 yy 

Fic. 3 


the characteristic vectors differ fundamentally from 6; and 2. If 63;*=(0,0, C3, Ds) is a right 
characteristic vector with the property 


[R exp — 0VR-! ]63” =d363" 


it does not follow that 63*[ ]=A363"; however, in our special case $3;”=(0, 0, C3, —D3) satisfies 
the equation $34[ ]= 363”. These vectors are biorthogonal to the other characteristic vectors 
(including 6,” and 64° which correspond to \,) in the manner 


6;"-6;=0; 6,*-6;"=0, | | 
44-93" =0; o34-6,"=0, (i=1, 2; 7=3, 4) 


and can be normalized so that 


o;%-9;"=1. 
The corresponding characteristic vectors of exp —6V are 
ts” = 634R=(C3, D3, —C3, —Ds3)/V2;  ts3*=(C3, —Ds, — C3, D3) /v2 
ty” = o4"R=(Cy, Ds, —Cs, —Ds)/V2; 3? =(Ci, —Ds, — Cs, Da) /v2. 


and 
By expanding 
exp —OV(a1, a2) = oi Ai” (an)Wi* (a2) 
(with yi“(a)=~i*(a)=Yla) and y2"(a) =y2"(a) =yY2(ax)) 
and applying the orthonormality of the characteristic vectors, 
Z= MN ye Wile) PAN? Dra Pola) PAN? LD! Yala) PMN?" DE Yala) P. 
A> Ae, As, Aa; ALV/2—HyDApN/A, AgN/2-1, Ny/2-H 
Z=MW*"C¥ ila)? and Z/¥=);} 


Since 


is the partition function per lattice point. Neglecting interactions between ‘‘next nearest’’ neighbors, 
f=0 and 
Ai =Amax = (e7 +e7")”, 
yielding 
Z"N =2 cosh J, 


which agrees with the known partition function of nearest neighbor linear chains.':? 
Since the average internal energy of spin interaction is 


E=0 log Z/d(—8), 
it can be found from our expression for Z. E/Ne is plotted in Fig. 1 for f=0 and f=e/2'. 
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V. APPLICATION TO TWO-DIMENSIONAL FERROMAGNETS 


As an application of Parts I and II, we shall calculate the high and low temperature partition 
function, internal energy, specific heat, and magnetization of two-dimensional NX M (as N, M—- =) 
ferromagnets. The definitions of interaction energies introduced in Part IV will be retained here. 
If our ferromagnet is inserted in a magnetic field H, the total energy of spin-field interaction is 


N 
—p(n,—n_)T= —pH > [n,(a;) —n_(a;) ], (28) 
i=1 


where 
#=magnetic moment of each spin, 
n4(a;)=number of ‘‘+”’ spins in the 7th row when it is in configuration a;, 
n_(a;)=number of “‘—”’ spins in the ith row when it is in configuration a;, 


n= >on,(a;) and n_=>on_(ai). 
We shall henceforth abbreviate C=yH/kT. 


1. High Temperature Development 


To derive a power series in 6=1/kT for Z that is valid at high temperatures, the two-dimensional 
analogue of (13) will be used. This necessitates the calculation of 


Sr=L-+-L exp —[V (ai, a2) + Var, a3) ++ +++ V(Qn-a, @n) |/RT (29a) 


and S,,4; with the value of x determined by the number of terms desired in the resulting power series. 

These S,’s correspond to partition functions of narrow strips of the sort exhibited in Fig. 2 in 
which the interactions between neighboring elements in chain a; (and likewise those between neigh- 
boring elements of chain @,) are only one-half those between other adjacent elements in the lattice. 
Instead of expressing S, as a sum over configurations of the vertical chains (i.e., a’s) in Fig. 3, 
one can equally well express it as a sum over configurations of horizontal chains (i.e., 6’s) 


Sr,=D---d exp —OLW(61, Bo) +--+ -+W(Bw-i, Bx) ], (29b) 
Bi BN 
n—l n 
W(Bi, B;) =} z v(oi, ky Gi, ra +3L0(o%.1, Gj, 1) +0(o;, ny Fj, n) J+  ¥ v(o;, ky Gj, k) 

k=1 k=1 
n—1 C n-1 

+3 DY (64%, 644) -— DL [ny(oi,n) +24(0;,4) —n_(o;,n) —n_(o;,x) ] 
k=1 Q k=2 


Cc 
~qlnete: 1) +n+(o;, 1) +n4(0;, a) +n4(o;, n) a n_(ai, 1) ie n(o;, 1) —n_(ai, n) = n_(o;, n) ] 


and o;, defines to the spin configuration of the /th element in the ith row. 


ni(o)=1 if o=0 (i.e., + spin), 
n4(o)=0 if o=1  (i.e., — spin), 
n(o)=0 if o=0, 
n(o)=1 if o=1. 


By expanding exp — W,,(8,8,) as a bilinear combination of the characteristic vectors of the equation’ 
u(n)o(n, Bi) = ~ exp — 0W,,(8:, B;)e(n, 8;), 
(29c) 
Su=u(0)menl La (0, 8). 


7 u(m) is a characteristic value corresponding to strips of width n. 
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If »=2, the exp —0@W matrix will be constructed so that the element in the ith row and jth column 
is the term 

exp —W(i—1, j7—1)8. 


The interaction energy in this term is that between a row in a configuration which corresponds to 
the two digit number (i—1) in the binary system and a neighboring row in a configuration corre- 
sponding to the two digit number (j—1). For example, 


W(0, 1) =[v(0, 0) +0(0, 0) +2(0, 1)+2(0, 1) ]/2— CL3n,(0) —n_(1) ]/4 


2 fam 
Q—1 Bo=1 
but 
v(0,1)=e, v(0,0)=—e, e=JkRT. 
Therefore W(0, 1) = —C/2. Similar calculations of the other matrix elements yield 
reste eCl2 Cle 1 
ecl2 1 ed 1 


exp — W0= 


ecl2 gi 





1 | 
| 
1 1 1 e2J-C J 


Expanding each matrix element as a power series in @ (i.e., in powers of J and C), choosing as a 
first approximation to the characteristic vector exp —W® the unnormalized vector of complete 
disorder (1, 1, 1,1) and finally proceeding by the method of successive approximations discussed 
in S.M.-I, Part-2, we obtain for the unnormalized characteristic vector 6(2) max of exp —Wé6 and 
for its characteristic value u(2) max 


[1 ] 
1—J—C/2+C/8+J?/2 
1—J—C/2+C?/8+ J?/2 
1—C+C?/2-—CJ J 

u(2)xam=4+3J?+C?+0(6). 


When n=3 the exp —W6 matrix is 


6(2) max = 




















e2C+ad e2J+3C/2 ef el e2F+3C/2 ef ef!2 1 ) 
e2J+3C/2 ected e724 +C/2 1 ef eC et eel2 

ec e72F+C/2 1 e~C/2 e724 +C/2 et e~Cl2 ec 

ef? 1 e~e/2 e2J—C e7 2d e2J—C/2 ec e24—3C/2 

exp — W6= | — 

e2J+3C/2 ef e724 +C/2 ed e2J+e ef!2 1 e~Cc/2 

ef eCl2 etd e-2J—C/2 ec? 1 e~2J—C/2 ec 

eC/2 e72d eel ee 1 e-2J—C/2 ey—ec e2J—3C/2 
1 eel? ec e2J—38C/2 e~ Cl ec e2J—3C/2 esJ—2€ 


The characteristic vector $(3)msx which corresponds to the largest characteristic value (3) max 
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(unnormalized) 
(1 


1—C/2-J+C/8+ 32/2 
1—2J—C+C?/2427 

| 1—3C/2—J+9C2/8—CI+J2/2 
1—C/2-J+C/8+F/2 
1—-2J—C+C?/24+ IC+2S2 
1—3J/2—J+9C?/8—CJ+ 32/2 
(1-2C+20?—3CJ 

(3) max =S+14.2+6C2+-++, 


6(3) max = 











and 


Then from Eqs. (13) and (29c) the first approximation to the partition function is (as M, N-« 
so that 17-—1-—M and N—1-—N) 


Z'WN-D ~ S3/ So +0(6*) 
and 


ZIMN ~ (3) max/p(2) max +0(6*) - 
30 
~24+27+C?+0(@). 
By the well-known relation between magnetization and the partition function 


mM/MN=p0 log Z/AC=uC+-:-:-. 


As the external field C--0, the magnetization per lattice point becomes zero and therefore in the 
temperature range where (30) is valid, a two-dimensional ‘‘ferromagnetic”’ system is not ferromagnetic. 
In the absence of a magnetic field we have calculated by the above method 


Zo BN ~ (5) max/ (4) max +0(6) _ 2(1 + J?+4J*/34+ Pee ), 


which agrees with the first three terms of the result of Kramers and Wannier* ® obtained by another 


method 


7 Zo MN =2(14+J2+4I4/34+77J°/454---). | (31) 
E=0 log Z/0(—0) and J=e/kT=€8, 

- E/MN~ — e(2J+10J*/3+64J°/15+---) (32a) 
C,/RMN~ J(2+10J2?+64J*/3+---). (32b) 


2. Low Temperature Development 


To derive the low temperature exponential development for Z we shall use Eq. (23) and therefore 
need 


P So.n=>D---d exp —LV(0, a1) + Ver, a2) +--+ + Van, 0) )/kT 
an a1 Gn 


Sin=L +++ DL exp —[V(0, a1) + Veer, a2) ++ +++ V(an, 1) J/RT. 


These S’s correspond to partition functions of narrow strips as shown in Fig. 3, in which the inter- 
actions between neighboring elements in the chains of 0’s and in the chain of 1’s is only one-half 
that between other adjacent elements in the lattice. Instead of expressing the S’s as sums over 
configurations of vertical chains (i.e., in terms of a’s) we shall express them as sums over horizontal 
(i.e., 8’ and y’) configurations. If we proceed in a manner analogous to that used in the low temper- 


* See also W. Opechowski, Physica 4, 181 (1937). 
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ature development 


. = N-1 1.0 N-1 ~—s 
So, at Si, nn, em & Pinax(B P+ n, re ® 2 Pmax(¥ )F, 


where », oM@max and », fmax are, respectively, the largest characteristic values of the matrices of inter- 


actions of two adjacent 8’ rows and between two adjacent y’ rows. Now 
n, max = exp [(2n+2)J+(n+1)C]+O(exp [(2n—6)J+nC)]), 
n, max =exp (2nJ+nC)+O(exp [(2n-—8)J+(n—1)C)). 
(n, 1@max/n, 0Mmax)* =exp —(2NJ+NC)+- 


So, 


as N-~«, and we can neglect S;,, as compared with So, ,,: 
So, nt Si, n= So, nL 1 +-O(e-?*I-*¢) ]. 


To find 1,o#max we need the exp —@W/(@';, B’2) matrix of the configuration 


. s.3 
0 x1 0 B’) 
0 Xe 0 B's 


with 
W(6, Bo) =v(0, x1) /2 +(x, 0) 2+0(0, Xe), ‘2 +(x, 0)/2+2(0, 0)/2 +v(0, 0)/2 +(x, Xe) —-C 


—[ (x1) —n_(x1) +04 (x2) — n_(x2) ]/2. 


This matrix is 


Xo\X1 0 1 
0 setvt2e ec 
exp —W0= ( ). 
1 ec 1 
By the method of successive approximations (choosing the ‘“‘ordered’”’ vector (1,0) as the first 
approximation 
PP ) ett +2e ef 
= etJ+20 
e-tJ—C 
gll+2C 9 
= (et +20 4 4) 
. ec 1 e-tJ-c e-t—C.4 98/30 F 
unt 


1, oMmax = C8 +204 eA 4 g-8J—204 


To continue further we must know 2, o4max and therefore develop the exp —@W/(8’;, B’2) matrix 


of the configuration 
QO xu Xp». O, 


0 Xo, Xo 0. 


Again all interactions except those between x’s of adjacent rows are one-half their normal value. 
The matrix is 


XoXo \ XX 00 A 3 7. 11 
01 e2J+2€ e2J+e ec 1 
=exp —V0. 
10 er +2C e-2s+e e2J+e 1 
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Starting with the “ordered” vector (1, 0,0, 0) as a first approximation the method of successive 
approximations yields the unnormalized characteristic vector 


(i 





e--C 4 e-8J-3C 4 Je-12J-5C _ 9-12-30 4... 
es —-C 4e-8J-3C 4 Jg-I2J-5C _ g-12J-3C 4... 


e 6/ “2 4 Qed iC 4 g—10J —6C 4 aie é 





and the corresponding characteristic value », ofmax is 


2, 0M@max = e6/ +8C 4 Je-2d+C 4 3e-6) aC W. 2e-10J-C 4 Ge-10d “IC Le 10/ ~5C 4 tee 


From Eq. (23) it is apparent that 


ZiU(M-) = So, n+1/ So, in 


and 


Z'M(M-1)(N-1) = 
asn—cx. 


Therefore as a first approximation 


n+1,0Mmax/ n, 0Mimax 


ZUM-DND Ao outmax /'1, oMmax = C2 +E(1 e826 4 Dee AC 4... ), 


To find the next approximation we must calculate 3, o4max from the matrix of the system 


O xn Xp 


QO X21 X22 


X13 «=O, 


X93 0, 


with all interactions except those between x’s of adjacent rows having one-half their normal values. 
It turns out by method of successive approximations that 


3,0éMmax > ESI FIC 4 36°C 4 Se-tC 4 eB ( ae 3 + 1 2g-*° 4-3e-*) + ihe 


Thus as a second approximation (neglecting the “ 


’ 


"in the terms N—1 and J/—1 since N, W—- =<) 


ZUNM ~ 3, oMlmax ‘2, omax = ¢2J+C/ 4 + e-8J—-2C 4 Je-t2J—4C 4 916-40 ( — 3+ 6e-2©+2e-1¢) + bwin ]. (33a) 


In the absence of an external magnetic field, 
C=0 and 


Zyl MN = @24 (1 + e-8/ 4 2e-12 + Se 4...) (33b) 


which agrees with the results of Wannier and 
Kramers, KW-II. 

For higher approximations it is necessary to 
find 4 oMémax, 5,0M@max, etc. by the same methods 
used for smaller values of m. In field free space 
4,0@max Was Of a value to make the next term in 
the Z/™ series 14e-!84. The low temperature 
magnetization per lattice point is 


m/MN =p log Z/aC 
= ul re Je-8J—2€ cet 8e—}247-4€ 


+ 2e-16J—-4€(7 — 18e 2 — Be-4C) 4+ --- J, (34) 





which shows that a two-dimensional ‘‘ferro- 
magnet” is actually ferromagnetic in the temper- 
ature range for which (34) is valid, in the sense 
that its magnetization is greater than zero in 
the absence of an external magnetic field. 

In zero external magnetic field and at low 
temperatures the average internal energy per 
lattice point is 


E/MN= —«(2—8e—™ 
~ De-12/ —72e/—.--) (35a) 
and the specific heat per lattice point is 
C/RMN=64J*e-* (14+-9e- /2 
+18e-*/+---). (35b) 
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Fic. 4. \max as a function of 1/J in the neighborhood of J-. 


3. Discussion 


Equations (32) and (33) provide us with 
expressions for the average energy and specific 
heat due to spin interactions at high and low 
temperatures. A reasonable problem now is to 
find the temperature ranges in which each 
expression is valid and to locate the temperature 
of a phase transition if one exists. In SM-I, 
p. 710, it is pointed out that a necessary condition 
for the existence of a phase transition at a 
temperature 7, is that the largest characteristic 
value of (4) be doubly degenerate. This condition 
of necessity results from supposing that the 
characteristic vector of (4) describing the most 
probable state of a layer of the system when 
T>T,. is te and when T<T, that it is qu. 
Then if \; is the characteristic value correspond- 
ing to wt, and Xz is the characteristic value 
corresponding to te 


Ar =Amax > Az if T< ln 
Ae =Amax > At if T> T~. 


Since at a transition point both phases (if there 
are two) are equally probable 


Ar =Ao=Amax if T= Ria 


Applying this criterion to the two-dimensional 
ferromagnet problem, if there is a phase transi- 
tion it must occur at the temperature at which 
the Zo!/"" of Eq. (33b) is equal to the Zo!" of 
Eq. (31). That is, at the temperature which 
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Fic. 5. Specific heat per lattice point in units of k. 


makes 


2(1+S2?+454/3+7758/45+---) 
=e? (1 +e-84 4 2-12 4 5e-16 4 24 e-20J 4. .), 


From Fig. 4 where the functions on both sides 
of this equation are plotted separately, the point 
of intersection of the two curves, J,.=0.44 
=e/kT., gives the temperature of a phase 
transition if one exists. At this point Z l/™% 
=2.52. Using a necessary condition that is a 
special property of a two-dimensional ferro- 
magnet for the existence of a sharp phase 
transition, Kramers and Wannier KW-I, II have 
calculated J.=0.4407 and Zo'!/"* =2.532 which 
agrees closely with our results. 

Unfortunately it is impossible from our 
equations to establish definitely the existence of 
a sharp phase transition at J,.=0.44. Since 
Zo!/™" as defined by (31) is equal to Amax when 
J>J., (32) gives the correct asymptotic expres- 
sion for the average internal energy and specific 
heat in this temperature range. When J<J,. 
(35) is a valid asymptotic expression. If the 
specific heat is plotted as a function of J (Fig. 5) 
the curves seem to be approaching a peak of the 
\-point variety at J=J,. Inasmuch as the series 
in (32) and (35) are probably asymptotic 
divergent series in the neighborhood of J, whose 
smallest terms (and therefore whose possible 
error’) are quite large, the calculated specific 
heat would probably be inaccurate around J.. 
As has been mentioned before, the magnetization 
at zero external field is zero in the region of 


*Knopp, Theory and Applications of Infinite Series 
(Blackie and Sons, 1928), p. 54. 
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validity of (30) while that in the region of 
validity of (34) is positive. Also, as J>J,, MM, 
as defined by (34) decreases; but (34) does not 
have enough terms to decide whether or not it 
decreases continuously to zero at J,. Approxi- 
mate calculations of higher order terms indicate 
that the SM series is probably an asymptotic 
divergent one around J,, so as in the case of the 
specific heat, this series gives there only a 
qualitative idea of the behavior of 9M. In spite 
of the inaccuracies near J,, our results as do 
those of Kramers and Wannier seem to suggest 
the existence of a A-point phase transition at J, 
that accompanies a change from a ferromagnetic 
state to an unferromagnetic state. 

In conclusion, it seems apparent that the 
characteristic value method of treating nearest 
neighbor problems is adequate to give asymptotic 
formulae for physical properties at high and 
low temperatures and to locate the temperature 
of a phase transition if one exists, but at its 
present state of development the method cannot 
provide accurate values of physical quantities 
which are derivatives of the partition function, 
in the immediate neighborhood of a transition 
point. 

The author wishes to thank Professor Lars 
Onsager for his enlightening discussions con- 
cerning this problem. 


APPENDIX” 
Proof that S2:,:/S2:=S141/S:+0(6') 
By definition 


t-1 
S:=2d-+-D Mexp—e@V (ai, ais). 
ay 


@,;i~1 


Therefore, if {As} is the set of characteristic values corre- 


0 The notation g(@) =O[f(@) ] as 6-0 is used here in the 
usual sense that there exists some positive finite constant 
A independent of 6 such that as 6-0. | g(@)/f(@)| finally 
becomes and remains less than A. 
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sponding to the set of characteristic vectors {tj} of the 
set of linear equations 


MN 


Av (ai) = 2 v(as)exp— 6V (a1, a) 


and if there are »v configurations available to each lattice 
point 
yun 
Si= 2 di! '[Lavi(a) F, 
i- 


Soi ‘ 1S7= Lidi*[Lavil(a) P+D’Ai?r,' “Savi (a) Fl<e¥; (a) } 
2 


(the prime on the 2 is introduced to dencte that the terms 
with i=j are omitted). Using a similar expression for 
SorSi41 


Sergi S1 — SoSi41 


MN 
= 2) Data tarts! Tavs (a) PLZ av (a) F. (36) 
t,97™= 

Suppose A; = 2a;;6/; A; =Amax- Then as 6-0 in the matrix 
exp—6V, \;—aiin. Since as 6-0 all the elements of exp—6V 
approach unity it is easily verified that \;=Amax—>v”%* 
and \;—0 if 7>1, or 4;=0(0). Also, when 7>1, either 
¥vi(a) =0(6) or for every a which makes yi(a@)=c there 
exists a unique 8 which makes ¥;(8) = —c; for if a single a 
yielded ¥i(a)=c and there existed no corresponding 8 to 

make yi(8) = —c, 


A= D Vilar) Vi(aslexp— OV (ara) 
bad ta | 
=c+0(6) 
contrary to our above conclusions. This implies 


vs 7.) sO) if i=1 
“adil@)=1 (9) if i>1. 
Since all 4; =0(6) except Ai which is 0(1), the dominant 
term in (36) as 6—0 is 
[ri24Aj2 4 — AA; Lavi (a) P[Zavi(a) FP (j¥1) 
= [0(1)0(6'*) —0(1)0(6") J[0(1) P[0(6) F 
=0(9'+1) —0(6"+2)-+0(6"+1), 


So we finally have 


Saigi Si — Sar St41 =0(0') 
or 
Sor41/S2i— St41/ Si =0(0') /S) S21 =0(0'), 


since S;S2—0(1) as 6—0, and our theorem is proven. 
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The Potential in a Solution Near a Charged Wall 


W. G. EVERSOLE 
Division of Physical Chemistry, The State University of lowa, 
Iowa City, Iowa 
December 1, 1941 
N recent papers in this Journal'? an equation was 
derived relating the potential, (y), in a binary salt 
solution near a charged wall to the surface potential, (Wo), 
of the wall and the distance, (X), from the Yo plane. The 
equation was derived without the usual simplifying as- 
sumption that se¥<CkT so it is valid for the relatively high 
potentials which may be encountered in electrokinetic 
measurements. 
This equation with numerical values substituted for 
water solutions at 25° may be written 


—log tanh (9.752’y) = —log tanh (9.752’yYo) 
+0.1422’(p’C)ix, (1) 


where y and yo are in volts, C in moles of salt per 1000 cc 
of solution, and X in angstrom units. 

In this form the equation is valid for binary salts and 
was shown to be a reasonable approximation for salts of 
other valence types if z’ and y’ are the valence and number 
of ions per salt molecule of the kind of ions having the 
opposite sign of charge to the solid wall. Equation (1) 
expresses y as an explicit function of X in which all of the 
quantities except Yo are known. 

Making the assumption that the electrokinetic poten- 
tial, (£), is the potential at the distance, ¢, from the yo 
plane, Eq. (1) may be written 


—log tanh (9.752’£) = —log tanh (9.752’yYo) 
+0.142t2’(v’C)*. 1’) 

It is obvious that if Yo and ¢ are independent of C, 
—log tanh (9.75z’£) is a linear function of C}. It is highly 
probable that the converse is also true since otherwise it 
would be necessary to assume that simultaneous changes 
in Yo and ¢ would exactly compensate. Thus if the experi- 
mental curve is linear both Yo and ¢ can be determined 
graphically. 

Further attempts to derive analogous equations for 
salts of other valence types have been successful in two 
other special cases: 

(1) s’=2s. This corresponds to a solution of a salt such 
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as BaCls or Th(SO,)2 in contact with a negative wall, or 
K2SO, or Ca2Fe(CN)¢ in contact with a positive wall. The 
equations analogous to Eqs. (1) and (1’) are 





log fe(y) =log fe(wo) +22’ CX (2) 
and 
log fe(é) =log fe(wo) +522’ Cit, (2’) 
where 
fo(y) = LEXP (s’ey/kT)+2 exp (s'ey/2kT) P4v3_, 2 
sit exp (2’ey/2kT)—1 v3’ 


e/kT =39 at 25° and with y in volts, 
v3 ( «Ne \* ; m 
b=, 55 (scoper) =0.123 for water solutions at 25°. 
(2) s’=}s. This case corresponds to a salt such as 
K2SO, with a negative wall or BaCls with a positive wall. 
The corresponding equations are 





log fs(w) =log fa(yo) +b32’C1X (3) 
and 
log f3(é) =log fs(Wo) +532’tC}, (3’) 
where 
f(y) = [2 exp (z’ey/RT)+1]}!+ v3 
J\0)* 02 exp (z’ey/kT) +1 — V3 
and 


b3 =2b2=0.247 for water solutions at 25°. 


Using data of Kruyt and van der Willigen* for a nega- 
tive glass wall, Yo and ¢ were evaluated graphically by 
Eq. (2’) for BaCl, solutions, and by Eq. (3’) for KsSO, 
solutions. Previously? these quantities were evaluated from 
the same data by Eq. (1’) which is an approximation for 
these salts. The results of the calculations were as follows: 

For BaCl. Eq. (2’) gave Yo=60 mv and ¢=31.4A, while 
Eq. (1’) gave yo=63 mv and t=31A. 

For K2SO, Eq. (3’) gave Yo=77 mv and t=29.4A, while 
Eq. (1’) gave Yo=77 mv and t=30A. 

In view of the excellent agreement of these calculations, 
it is doubtful if the accuracy of existing electrokinetic 
potential data will justify the additional inconvenience 
involved in making calculations by means of the more 
exact equations. 

1 W. G. Eversole and P. H. Lahr, J. Chem. Phys. 9, 530 (1941). 

2 W. G. Eversole and W. W. Boardman, J. Chem. Phys. 9, 798 (1941). 


3H. R. Kruyt and P. C. van der Willigen, Kolloid Zeits. 45, 307 
(1929). 





Ethylene-Ethane and Propylene-Propane 
Equilibria 
GEORGE B. KISTIAKOWSKY 
Department of Chemistry, Harvard University, Cambridge, Massachusett 
December 1, 1941 

E have determined the thermal equilibrium con- 

stants in the reactions C2Hs=C2H.+He2 and 

C;Hs=C;He+Hs at two temperatures each and give the 
results in Table I. 

The equilibrium was approached from both sides, using 
Cr.O3 as a catalyst.! It was found essential to have the 
catalyst aged, as otherwise less consistent data were ob- 
tained. The gases were recirculated over the catalyst and 
the position of equilibrium—once it was established—was 
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found to be independent of the number of passages over 
the catalyst and the initial gas composition.* 

A special analytical procedure was developed which 
enabled us to measure accurately the quantity of hydrogen 
and total olefine in the gas mixture, and also to determine 
the amount of ethylene or propylene in the olefinic fraction. 

Applying the van’t Hoff isochore to the data of Table I, 
the heats of the two reactions are calculated to be 34.0 x 10° 
calories per mole at 415°C, and 30.8 X 10° calories per mole 


TABLE I. 





Initial Number of Equilibrium 
excess experi- constant, Average, 
Temp. of: ments atmos. atmos. 
450°C CoHe 4 5.12 +0.09 K10~4 
5.16 X10 
CoHs+He 4 5.21 +0.04 1074 
380°C CoHe 3 3.94 +0.04 X107> 
4.04 X10 
CeHe+He 3 4.15 +0.13 K10™ 
375°C C3Hs 3 5.25 +0.08 X10~4 
5.25 X10~4 
CsHe+He 4 5.26 +0.01 X1074 
310°C C3Hs 1 3.69 X10~5 
3.68 X10~ 
CsHe+He 1 3.67 X1075 








at 343°C. These values differ from those measured directly 
by Kistiakowsky et al.2* at 82°C and extrapolated to the 
higher temperatures by means of the available heat 
capacity data (33.9108 and 31.1 10%, respectively) by 
only a few hundred calories, which can be caused by 
errors in the equilibrium constants of less than 2 per- 
cent. Errors of this magnitude we consider inherent in our 
measurements. 

We have calculated statistically equilibrium constants 
for the ethane-ethylene-hydrogen system, using the heat 
of reaction as given in the reference above, and molecular 
constants as given by Stitt;' Frost;5 Burcik, Eyster, and 
Yost;® Smith;? Giauque;* and Davis and Johnston.® 
The result is Kp=5.19X10-* at 450°C, for a potential 
barrier restricting internal rotation equal to 2750 calories 
per mole, as given by Kistiakowsky, Lacher, and Stitt.!° 
The equilibrium constant for .the propane-propylene- 
hydrogen system has also been calculated statistically, 
using the heat of reaction as given by Kistiakowsky et al.* 
The vibrational frequencies of propylene were taken from 
Wilson and Wells;" those of propane from Wu and 
Barker,” except that the four undetermined frequencies 
were taken as 2965, 1465, 1150, 1100. (This assignment 
gives heat capacities in satisfactory agreement with the 
experimental data reported by Kistiakowsky and Rice.!*) 
Taking as the potential barriers restricting internal rota- 
tion 3300 calories for propane and 2100 for propylene, we 
obtain Kp=5.6X 10 at 375°C, using Pitzer’s tables for 
the contributions from internal rotation. 

The excellent agreement of calculated and experimental 
equilibrium constants for the ethane reaction shows how 
trustworthy can be the statistical calculations although 
merely the simple rigid rotator and harmonic oscillator 
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molecular models are used, even at high temperatures, 
provided the molecular constants are accurately known, 
as is the case for this reaction. 

The good agreement of calculated and experimental 
equilibrium constants for the propane reaction is an addi- 
tional strong evidence, if such be needed, that propylene 
has a potential barrier of about 2100 calories and not, as 
Pitzer claims, 800. Calculation of the equilibrium con- 
stant by Pitzer’s assignment!® gives the value 10.0 10~4 
at 373°C. 

Detailed description of our work, as well as a discussion 
of the literature on the subject, will be published in this 
Journal with some delay because of pressure of national 
defense work. 


1 Frey and Huppke, Ind. Eng. Chem. 25, 54 (1933). 

* The slight effect of the initial gas composition on K evident from 
Table I is due to the inclusion in the averages of earlier runs made with 
an insufficiently aged catalyst. 

2 Kistiakowsky, Romeyn, Ruhoff, Smith, and Vaughan, J. Am. Chem. 
Soc. 57, 65 (1935). 

3 Kistiakowsky, Ruhoff, Smith, and Vaughan, J. Am. Chem. Soc. 
57, 876 (1935). 

4F. Stitt, J. Chem. Phys. 7, 297 (1939). 

5 A. V. Frost, Comptes, rendus, U.S. S. R. 161 (1933). 

6 Burcik, Eyster, and Yost, J. Chem. Phys. 9, 118 (1941). 

7L. G. Smith, J. Chem. Phys. 8, 798 (1940). 

8 W. F. Giauque, J. Am. Chem. Soc. 52, 4816 (1930). 

% Davis and Johnston, J. Am. Chem. Soc. 56, 1045 (1934). 

10 Kistiakowsky, Lacher, and Stitt, J. Chem. Phys. 7, 289 (1939). 

1 E. B. Wilson, Jr. and A. J. Wells, J. Chem. Phys. 9, 319 (1941). 

122V, L. Wu and E. F. Barker, J. Chem. Phys. 9, 487 (1941). 

13 G. B. Kistiakowsky and W. A. Rice, J. Chem. Phys. 8, 610 (1940). 

4K. S. Pitzer, J. Chem. Phys. 5, 469 (1937). 

8 K. S. Pitzer and W. D. Gwinn, J. Chem. Phys. 9, 485 (1941). 

i6 K. S. Pitzer, J. Chem. Phys. 5, 473 (1937). 





Equilibrium in the Gas Reaction, 
CHBr;+Br.—CBr,+HBr 


JOHN W. STRONG AND ROBERT N. PEASE 
Frick Chemical Laboratory, Princeton University, Princeton, New Jersey 
December 11, 1941 
ECENT calculations of partition functions for the 
bromomethanes! with other pertinent data indicate 
that such substitution reactions as 


CH.Br2+ Br2>—CHBr;+HBr (1) 
and 


CHBr; +Bro—CBr;+HBr (2) 


should be incomplete in the gas phase. This has been con- 
firmed by preliminary determinations of equilibrium in the 
latter reaction. 

Samples of liquid bromoform and bromine were weighed 
out in glass ampoules. These were introduced into a 
Pyrex tube which was then evacuated, sealed, and heated 
in a vapor bath (either naphthalene b.p. 218°C, or a-bromo- 
naphthalene, b.p. 281°C). Total pressures in the reaction 
tubes at the operating temperatures were about 1 atmos. 

After heating, the tubes were quenched in a water spray, 
opened under potassium iodide solution, and the solution 
titrated with thiosulphate and with alkali. Data on hydro- 
gen bromide formed were usually a littler lower (1-3 per- 
cent) than for bromine reacting, indicating possible side 
reactions other than substitution. No appreciable amounts 
of permanent gas were detected. 
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Principal sources of uncertainty follow from the facts 
that equilibrium was only approached from the one side, 
and that no allowance could be made for the formation 
of lower bromides (especially methylene dibromide). 
However, it was shown that variation of the heating time 
(0.5 to 4 hours) indicated a quite definite end-point short 
of complete reaction. Further, by increasing the ratio of 
bromine to bromoform up to as much as 5 to 1, the forma- 
tion of lower bromides could be minimized. 

Results calculated for the equilibrium constant, 
(CBr, ][HBr]/([CHBr;][Br2], gave values of 0.30 to 
0.32 at 218°C, and 0.26 to 0.30 at 283°C, (six experiments 
each) with no definite trends. 

An equilibrium constant was calculated from the sta- 
tistical data at 500°K (227°C) for purposes of compari- 
son. For the function, (Eo°— F°)/T, Stevenson and Beach! 
find 73.56 and 78.43 at 500°K, for CHBrs and CBr, 
respectively. Gordon and Barnes? find corresponding 
values for Bre and HBr to be 55.01 and 44.14, respectively. 
Thus, at 500°K (AEo°—AF*)/T = —6.00. 

The value of AEo® is uncertain. Judging from heats of 
formation at 18°C, as estimated by Bichowsky and 
Rossini,? AE29:°= —2650 cal. Calculated heat capacities 
indicate that this should be increased to about —3000 
cal. at O°K. 

Substituting this value for AEo®, one obtains AF50°=0, 
when Kso00=1, which is in fair agreement with the experi- 
mental data in the corresponding temperature range. A 
value of AEo® nearer — 1800 cal. would yield better agree- 
ment, and would be in better accord with the small ob- 
served temperature coefficient of the equilibrium constant. 

To this extent, the calculations of Stevenson and Beach 
are confirmed. 


1D. P. Stevenson and J. Y. Beach, J. Chem. Phys. 6, 25, 341 (1938). 

2 Gordon and Barnes, J. Chem. Phys. 1, 692 (1933), 

3 Bichowsky and Rossini, Thermochemistry (Reinhold Publishing 
Company, New York, 1936). 
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Erratum: The Partition Function of Molecules 
with Internal Torsion 


I, Single Asymmetric Top Attached to Rigid Framework 
[J. Chem. Phys. 9, 807-815 (1941) ] 


DONNA PRICE 
Mallinckrodt Chemical Laboratory, Harvard University, 
Cambridge, Massachusetts 
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(sub- 

script 

gives 
column) Line Corrected form 

810: _—_‘ Fine print, sixth equation. n?=C2+E?—22 edradyaEC =E2+(?. 

8101 Fine print, two lines above (14) Yada.sva?<0.058; 2M-'=0.013. 
~ 8102 1 (13) instead of (14). 

812) 4 and 6 h inserted before V2. 

8121 12 QzqC divided by ga. 

812; Eas. (29) and (30) R2, I2, S? instead of R2, I?, 82. 

812: Eq. (29), last term 3 instead of 5. 

812; Third line below (30) S=Cle4a4:0:0Pa- 

812+ Eq. (33), line 1 Do instead of Di. 

813; 2 Do instead of D,. 

8132 3 o =AC* Ze AessAs:0°. 

814> 8 


+ precedes sum. 








The last part of footnote 12, p. 814, should read: 
M+: and Mz are the same as above except for the factors 
of 4 and }4, respectively, instead of 3. 

Nii = 6D: — AAA; Vii, 

[gigs ¥+ 26 Yif; 1, i, jS3 
Cairgi VAL (qirdsy +9):03i') + Ceasiras;" 


7 +48 Ul], i,7>3 
[eign Y+Legiasy + 28 Yily J, i333, j>3 
where Y=(1—38+ 387), and L=C(1—467). 
(eigiE+68Yifil, i, jS3 
ng g J [gig E— 2B (Girdsj +g jai") C 
alin +B], i,7>3 
[gig7E— $B g:a37°C+B Yl; ], 433, j>3 


where E=(38-?—8"). F;; should also include the factor A 
as does E;j. 





